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Atmospheric Flight Dynamic Simulation Modeling 
Of Spin-Stabilized Projectiles 


Dimitrios N. Gkritzapis,* Elias E. Panagiotopoulos,” Dionissios P. Margaris” and Dimitrios 
G. Papanikas? 
* Athens, Greece, GR 11522 
> Fluid Mechanics Laboratory, Mechanical Engineering and Aeronautics Department, University of Patras, Patras, 
Greece, GR 26500 


Abstract: A full six degrees of freedom (6-DOF) simulation flight dynamics model is applied for the accurate 
prediction of short and long range trajectories of high and low spin-stabilized projectiles and small bullets 
via atmospheric flight to final impact point. The projectile is assumed to be both rigid (non-flexible), and 
rotationally symmetric about its spin axis launched at low and high pitch angles. The projectile maneuvering 
motion depends on the most significant forces and moment variations, in addition to wind and gravity. The 
computational flight analysis takes into consideration the Mach number and total angle of attack effects by 
means of the variable aerodynamic coefficients. For the purposes of the present work, linear interpolation 
has been applied taking data from official an tabulated database. The aforementioned variable flight 
model is compared with a trajectory atmospheric motion based on appropriate constant mean values of 
the aerodynamic projectile coefficients. Static stability, also called gyroscopic stability, is examined as a 
necessary condition for stable flight motion in order to determine the sufficient initial spinning projectile 
rotation. The efficiency of the method developed gives satisfactory results compared with published data of 
verified experiments and computational codes on atmospheric dynamics model flight analysis. 


Keywords: Constant and variable aerodynamic coefficients, low and high launched angles, Magnus effects, 
symmetric projectiles dynamics, stability criteria, trajectory dynamics prediction 


addition to wind and gravity forces. The variable 
aerodynamic coefficients are taken into account 
depending on the Mach number of flight and total 
angle of attack. 


Introduction 


Ballistics is the science that deals with the motion 
of projectiles. The word ballistics was derived 
from the Latin ‘ballista,’ which was an ancient 


machine designed to hurl a javelin. The modern 
science of exterior ballistics! has evolved as a 
specialized branch of the dynamics of rigid bodies, 
moving under the influence of gravitational and 
aerodynamic forces and moments. 


Pioneering English ballisticians Fowler, Gallop, 
Lock, and Richmond? constructed the first 
rigid six degrees of freedom projectile exterior 
ballistics model. Various authors have extended 
this projectile model for lateral force impulses** 
dual-spin projectiles*® etc. 


The present work addresses a full six degrees 
of freedom (6-DOF) flight dynamics analysis 
for accurate prediction of short and long-range 
trajectories of high spin-stabilized projectiles 
and small bullets. The proposed flight dynamic 
model takes into consideration the influence of the 
most significant forces and moment variations, in 
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Additionally, this flight trajectory analysis is 
compared with another flight model based 
on appropriate constant mean values of the 
aerodynamic coefficients. The efficiency of the 
method developed gives satisfactory results 
compared with published data of verified 
experiments and computational codes on dynamics 
model flight analysis of short and long-range 
trajectories of spin-stabilized projectiles. The 
present analysis considers the cartridge 105 mm 
HE M1 with a 105 mm howitzer such as the M103 
with M108 cannon, as a representative projectile 


type. 
Projectile model 


A typical formation of the Cartridge 105mm 
HE MI projectile is presented (see Figure 1), 
and is used with various 105 mm Howitzers (see 
Figure 2) such as US M49 with M52, M52A1 
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Figure 1. 105 mm HE M1 high explosive 
projectile artillery ammunition for howitzers. 


Figure 2. Howitzer of 105 mm projectile at firing 
site. 


cannons, M2A1 & M2A2 with M101, MIOIAI1 
cannons, M103 with M108 cannon, M137 with 
M102 cannon as well as NATO L14 MOD56 and 
L5 (see Table 1). 


Cartridge 105mm HE M1 is a semi-fixed type 
ammunition, using adjustable propelling charges 
in order to achieve desirable ranges. 


The projectile producing both fragmentation and 
blast effects can be used against personnel and 
materials targets. Technical characteristics and 
operational data of 105 mm HE M1 are illustrated 
(see Table 2). 


Physicomathematical flight 
simulation model and simplifications 


A six degrees of freedom rigid-projectile model 
79.11 has three rotations and three translations, the 
three translation components (x,y,z) describing the 
position of the projectile’s center of mass and the 
three Euler angles (¢,0,y) describing the orientation 
of the projectile with respect to translation from 
the body frame (no-roll-frame, NRF, ¢ = 0) to 
the plane fixed (inertial frame, if). For such flight 
bodies, the Xypp axis of the projectile in the no- 
roll-frame coordinate system usually lies along 
the projectile axis of symmetry and the Yxpp and 
Zyrp axes are then oriented so as to complete a 
right-hand orthogonal system (Figure 3). 


Table 1. Maximum range impact of 105 mm HE M1 depending on cannon charge and muzzle velocity at 


firing site. 
Cannon charge M52 MS52A1, M101, MIO1A1 M102 M103 
Cannon charge 1 7 1 7 
Muzzle velocity/m s ' 198.1 472.4 205.0 494.0 
Maximum range/m 3510 11270 3700 11500 


Table 2: Technical characteristics and operational data limits of 105 mm HE M1 howitzer projectile. 


Technical description 

The cartridge 105 mm HE M1 consists 
Projectile: Hollow steel forging 

High explosive charge: TNT (2.09 kg) 
Supplementary charge: TNT (0.14 kg) 
Cartridge case: M14 


Propelling charge: M67 (7 propelling charge increments) 


Percussion primer: M28B2 
Fuse: PD M557P1 or equivalent 
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Technical data 


Total length: 790 mm 
Total weight: 18.15 kg 


Temperature limits: 


Operation: —40 °C to +52 °C 
Storage: —62 °C to 71 °C 
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Figure 3. No-roll (moving) and fixed (inertial) 
coordinate systems for the projectile trajectory 
analysis. 


Newton’s laws of motion state that the rate of 
change of linear momentum must equal the sum 
of all the externally applied forces and the rate of 


change of angular momentum must equal the sum 
of all the externally applied moments, respectively. 
The force acting on the projectile comprises the 
weight W; , the aerodynamic force A, and the 
Magnus force M;. 


The moment acting on the projectile comprises 
the moment due to the standard aerodynamic 
force A,,, the Magnus aerodynamic force M,,, 
and the unsteady aerodynamic moment UA,,. 
All aerodynamic coefficients are based on Mach 
number and the aerodynamic angles of attack and 
sideslip. 


The independent variable is changed from time to 
dimensionless arc length and the equations that 
determine the pitching and yawing motion become 
independent of the size of the projectile which 
turns out to be very convenient in the analysis of 
free-flight range data. 


The transformations to the no-roll-frame (~) give 
us equations (1)-(12). 


x’ =Dd-—wWit+woe (1) 
2 D 

ee ea a (2) 

z'=—-D0+—wW (3) 
De 

o' =— p+— tan|6| 7 (4) 
7 
D 

'=—G 5 
yf (5) 

eee 7 (6) 
V cos 

2 D 1 7 w 7g vee 2D 
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The projectile dynamic model consists of 12 
highly first order non-linear ordinary differential 
equations, which are solved simultaneously by 
resorting to numerical integration using the fourth 
order Runge—Kutta method. 


In these equations, the following sets of 
simplifications are employed: the aerodynamic 
angles of attack a, 6 are small so 


arXwlV,BxvlV 


Also the yaw and pitch angles y, @ are small. A 
change of variables is to introduce velocity u 
replaced by total velocity V because the side 
velocities ¥ and W are also small. The Magnus 
force and moment are calculable quantities in 
comparison with the weight and aerodynamic 
forces so they are maintained in the forces and 
moments, respectively. The projectile is 
geometrically symmetrical so: 
Dy =Lyz =1yz, =9 and Ly =Iy 

The projectile is aerodynamically symmetric so 
the expressions of the distances from the center 
of mass to both the standard aerodynamic and 
Magnus centers of pressure are simplified. 


Static or gyroscopic stability 


Any spinning object will have gyroscopic 
properties. In a spin stabilized projectile, the 
center of pressure CP, the point at which the 
resultant air force is applied, is located in front of 
the center of gravity CG. Hence, as the projectile 
leaves the muzzle it experiences an overturning 
movement caused by air forces acting about the 
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center of mass. It must be kept in mind that the 
forces are attempting to raise the projectile’s axis 
of rotation. 


In Figure 4 two cases of static stability are 
demonstrated: in the top figure, CP lies behind CG 
so that a clockwise (restoring) moment is produced. 
This case tends to reduce the yaw angle and return 
the body to its trajectory, therefore statically stable. 
Conversely, the lower figure, with CP ahead of CG, 
produces an anti-clockwise (overturning) moment 
which increases « further and is therefore statically 
unstable. It is also possible to have a neutral case 
in which CP and CG are coincident whereby no 
moment is produced. 


There is clearly an important correspondence 


statically stable 


™, 
ei 


Figure 4. Static stability/instability conditions. 
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Oeil 


“trajectory 


yaw angle a 


Figure 5. Gust producing yaw angle on projectile 
in flight. 

between the distance from the centre of pressure to 
the centre of gravity and the static stability of the 
round. This distance is called the static margin. By 
definition, it is positive for positive static stability, 
zero for neutral stability and negative for negative 
stability. 


Firstly, we shall consider the case of a shell-like 
projectile in flight (see Figure 5). This is initially 
flying at zero yaw incidence along its flight 
trajectory and is then struck by a gust of wind so 
that the nose is deflected upwards, producing a 
yaw angle (a). The response of the projectile to 
this disturbed yaw angle determines its stability 
characteristics. In particular, the initial response 
determines whether it is statically stable or 
unstable. 


If the initial response is to move the nose back 
towards zero yaw (i.e. reducing the yaw angle) 
then it is statically stable. If the yaw angle initially 
increases as a response then it is statically unstable 
while if the disturbed yaw angle is retained the 
shell is neutral. It is therefore clear that it is the 
direction or sign of the resulting yawing moment 
generated that defines the static stability. This 
depends upon the aerodynamic force produced on 
the body due to the yaw angle and, in particular, 
upon the normal force (i.e. the aerodynamic force 
component perpendicular to the body axis) and the 
position at which it acts along the body’s axis. 


Classical exterior ballistics’ defines the gyroscopic 
stability factor S,, as: 
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Table 3. Physical and geometrical 
characteristics of the 105 mm projectile. 


105 mm, HE, M1 projectile 


Reference diameter 104.8 mm 
Projectile total length 49.47 cm 
Projectile weight 15.00 Kg 
Axial moment of inertia 0.02326 Kg m? 


Transverse moment of inertia 0.23118 Kg m* 


18.34 cm from the 


Center of gravity base 


Table 4. Trajectory aerodynamic parameters of 
atmospheric flight dynamic model. 


105 mm, HE, MI projectile with constant 
aerodynamic coefficients 


Cy = 0.243 C,=1.76 Cz» =—0.0108 
Cig =—9:300 Cy = 3.76 Cypy = —0.381 
Cyp, = 0.215 
pi 
x DP (13) 


®~ 2pl,SdV°C,,, 


The static factor is proportional to the product 
of four terms, which are the geometric technical 
characteristics of the projectile shape model, the 
square axial spin to velocity ratio, the aerodynamic 
Magnus effect variation and the proposed density 
atmospheric model. 


This may be rearranged into: 


21 Noe 1 yl 
| Fes a 


3 
Ind 

The estimation of the projectile stability factor is 

effected from the last three terms separately, due to 

the fact that the first term is a constant parameter 

for a given projectile body motion. 


& 


Computational simulation 


The flight dynamic model of the 105 mm HE, M1 
projectile involves the solution of the set of the 
twelve non-linear first order ordinary differential 
equations(1)—(12) whichare solved simultaneously 
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——yvariable model [70° intial pitch 
constant model (70° initial pitch 


—x—McCoy Model 


O° initial pitch 
— MeCes Model 


——variable model (45° initial pitch 
——constant model (45° initial pitch 
al initial pitch 


altitude, m 


0 2000 4000 += 6000» 8000» 10000 = 12000 
range, m 


Figure 6. Impact points and flight path 
trajectories with constant and variable 
aerodynamic coefficients for a 105 mm projectile 
compared with McCoy 8 trajectory data. 


by resorting to numerical integration using a 
fourth order Runge-Kutta method. Physical 
characteristics of the projectile 105 mm, HE, M1 
are listed (see Table 3). 


The constant dynamic flight model uses mean 
values of the experimental aerodynamic coefficient 
variations! (see Table 4). 


For the trajectory parameters of the projectile with 
variable coefficients, linear interpolation has been 
applied taking data from the official tabulated 
database.’ Initial data for both dynamic trajectory 
models with constant and variable aerodynamic 
coefficients are: 


x=0.0 m, y= 0.0 m, z= 0.0 m, ¢ = 0.0°, = 45.0° 
and 70.0°, y=0.0°, 7 =494ms', > =0.0ms1, 

=0.0ms-l, P = 1644 rads 1q = 0.0 rads! 
ue r =0.0 rad s |. Giving the initial rifling twist 


rate 7 at the gun muzzle (calibers/turn) the axial 
spin rate is calculated from the expression: 


p=2nV/nD (rads') (15) 


The density and pressure are calculated as a 
function of altitude from the simple exponent 
model atmosphere, and gravity acceleration taken 
with the constant value g= 9.80665 ms * 


Results and discussion 


The flight path trajectories of the present dynamic 
model with initial firing velocity of 494 ms ' and 
rifling twist rate 1 turn in 18 calibers (1/18) of 
the M103 Howitzer at initial pitch angles of 45° 
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—45deg 9 ——7Odeg 


Square axial to velocity ratio 


0 2000 4000 6000 8000 10000 =12000 
range, m 


Figure 7. Square axial spin to velocity ratio 
versus range at quadrant elevation angles of 45 
and 70 degrees. 


— 45deg ——7Odeg 


+ + — — + | 
0 2000 4000 6000 8000 10000 12000 


range, m 


Figure 8. Inverse of density versus range at 
quadrant elevation angles of 45 and 70 degrees 
for 105 mm projectile. 


—45deg + 7Odeg 


Inverse Pitching Moment 


0,2 04 0,6 0,8 1 1,2 14 1,6 
Mach Number Flight 


Figure 9. Inverse pitching moment versus Mach 
number flight at quadrant angles of 45 and 70 
degrees. 
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and 70° are indicated in Figure 6 for two cases: 
constant and variable aerodynamic coefficients. 
The impact points of the above trajectories 
with the present variable coefficients method 
compared with accurate estimations of McCoy’s 
flight atmospheric model’ provide satisfactory 
agreement for the same conditions. The diagram 
shows that for the 105 mm MI projectile, fired 
at sea-level with an angle of 45° (dark blue solid 
line) and no wind, the predicted range to impact 
is approximately 11520 m, the time of flight is 
slightly greater than 53 s, and the maximum height 
is 3465 m. At 70° (green solid line), the predicted 
level-ground range is 7340 m, the time of flight to 
impact is about 72.5 s, and the maximum height is 
slight over 6030 m. In the same diagram, the basic 
differences from McCoy’s flight data are remarked 
in trajectory models with constant aerodynamic 
coefficients. 


The three final bracketed terms of the gyroscopic 
stability factor (equation (14)) contain trajectory 
variables which vary significantly during a 
projectile’s flight and each will be considered 
independently below: 


(a) Square axial spin to velocity ratio effect 


The ratio of square axial spin to velocity (Figure 7) 
against range at 45°and 70° increases from a value 
of 11 at the muzzle, to values of 49 and 204 at the 
apogee, and then decreases to values of 22 and 19 
at impact, respectively. 


For the majority of projectiles, the velocity V 
will decay at a faster rate than the spin rate p. 
Consequently the launch conditions will be critical 
— S, will increase (i.e. improve) as the projectile 
moves down-range. 


It can be stated that if'a projectile is statically stable 
at the muzzle, it will be statically stable for the 
rest of its flight. This we can understand because 
the static factor is proportional to the ratio of the 
projectile’s rotational and transverse velocity. As 
the rotational velocity is much less damped than 
the transverse velocity (which is damped due to 
the action of the drag), the static factor increases, 
at least for the major part of the trajectory. 


(b) Atmospheric density effect 


The inverse of density against range at 45° and 70° 
(Figure 8) rises from 0.8 at the muzzle, to values 
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0 + + — + + 


0 2000 4000 6000 8000 10000 12000 
range, m 


Figure 10. Comparative static stability on the 
range with constant and variable coefficients at 
quadrant angle of 45 degrees. 


0 2000 4000 6000 8000 
range, m 


Figure 11. Gyroscopic stability factor variations 
of constant and variable coefficients at 
70 degrees initial pitch angle. 


of 1.15 and 1.52 at the apogee, and then decreases 
to a value of 0.8 at the point of impact. 


Also S, will be at its lowest when the density is 
at its highest. This will occur at low altitudes and 
at low temperature, high pressure atmospheric 
conditions. It is for this reason that firing trials 
for projectile impacts are often carried out in cold 
countries. 


(c) Magnus aerodynamic effect 


On the other hand the aerodynamic Magnus effect 
coefficient (Figure 9) clearly varies with anything 
which affects static margin. As a consequence, any 
variations in the positions of either CP or CG will 
directly affect the value of Cy,,. The position of the 
CP will certainly vary during flight if the projectile 
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passes through the transsonic flow regime, due to 
complicated shock and expansion fan movements. 
In particular, it rises abruptly as the inverse Cy, 
term is highly reduced. This therefore means that 
S, reduces under transsonic conditions. 


The gyroscopic stability factor with constant 
and variable aerodynamic coefficients at 45° 
(Figure 10), which was 3.1 at the muzzle, increases 
to 23 and 19 at the summit of the trajectory and 
then decreases to values of 8 and 5.8 at impact, 
respectively. Also estimates of S, have been 
made for aerodynamic coefficients with constant 
and variable values at 70° (Figure 11). It begins 
with an initial 3.1 at muzzle (the same value as at 
45°), rising to 123 and 117 at the summit of the 
trajectory and then decreasing to values of 6.6 and 
4.19 at the point of impact, respectively. 


Conclusion 


The six degrees of freedom (6-DOF) simulation 
flight dynamics model is applied for the accurate 
prediction of short and long range trajectories 
of high and low spin-stabilized projectiles and 
small bullets. It takes into consideration the Mach 
number and total angle of attack variation effects 
by means of the variable and constant aerodynamic 
coefficients. The criteria and analysis of gyroscopic 
stability are also examined. The computational 
results are in good agreement compared with other 
technical data and recognized projectile flight 
dynamic models. 
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Appendix — Equations of motion, forces and moment and symbols 


The dynamic equations of motion are derived in the non-rolling frame and provided in equations (1) 
through (4): 


Xif cosAcosy —siny sin cos] | typr 
Vy-p=lcosPsiny cosp = sindsinwW}} Vy, (1) 
Zi —sin0 0 cos 0 Wyer 


For position of the projectile and for orientation with Euler angles help 
0 ly Pyer 


1 
=|0 1 0 4 ver (2) 
0 0 1/cos@)| Tyee 


| DSi 


From the second law of Newton we have the force and the moment which act on the projectile in equations 
(3) and (4), respectively. 


Une F- XTorAL /m 0 r; NRF —4 ver Uner 
Vue f = 7 Vrora (tt +|—Fyee 0 —Typr |) Vine (3) 
Wyer FZror4, 1m Quer ver ty 0 Wyer 
Leora 0 =, NRF Auer \|Lxx xy 4¢yz|| Pryer 
-1 fe Bs se a 
[Z M toTaL {| “ner 0 Furr to\|Lix Ley yz || wre (4) 
N. TOTAL — ner =F; wer £9 0 Lox ZY ZZ i; NRF 


The force acting on the projectile in equation (3) comprises the weight W;, the aerodynamic force A,and 
Magnus force M;: 


Prior, | |Xw,| |X 


4 M 
FY tora, {| = Yw, a Yy +5 Vy (5) 
E osiae Zw, Zf Ly 


The moment acting on the projectile in equation (4) comprises the moment due to the standard aerodynamic 


force A,,, due to the Magnus aerodynamic force M,, and the unsteady aerodynamic moment UA,,: 
Lrorat Ly Lig Lis 
Mora | =)! +) Miag( +) Mw (6) 


All aerodynamic coefficients are based on the Mach number and the aerodynamic angles of attack and 
sideslip: 


a=tan! 


Wyre 
a 1) 


Uner 
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G=tan"' 


Vivre 
oe | (8) 
u NRF 


The total aerodynamic is velocity given in equation (9): 


V, = Ties + ees a Wrage (9) 


The weight force in the no-roll system is: 


Xw, —sind 
Yw, =mg; 0 (10) 
Zw, cos0 


The aerodynamic force, which acts on the projectile at the aerodynamic center of pressure, is: 


~ 2 ~ 2 
Ww Vv 
Cy = Ce A =Cys Ze 
X Ee T 
a i 
Vy t= = OV 7 Srey =O, (11) 
Z 5 V, 
Zz. - 
, Ww 
—Cy, NRF 
V, 


The Magnus force, which acts on the projectile at the Magnus force center of pressure, is: 


Xy 1 D i 
yy. = — p78 Purr? UnpaWner 12 
iM 3? T ref 2V,? (12) 
M e by 
—PyperD Cypser 
2,7 
The moment due to the aerodynamic force is: 
L,| [0 0 0 (xX, 
Myt=|0 0  —Rowc|| Vy (13) 
Ne WO Res O Wize 
The moment due to the Magnus force is: 
i,| [0 0 0 |[X,| 
Myt=|0 0 =Reyaz||Yy (14) 
Nel 10> Riss 0 llZe 


In addition, for the unsteady moment UA,, is: 
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PwrrDCip 


i 2V 
WA ~ 
~ I Suey Gurr PCy Cc 
M wut =~ PV,D' | -—+ + 4 15 
. 2°? 2 V ee 
uA Furr DCyo 7 Cy 
2V V 
Forces and moments 
Dd 
a Pp e Siop&y 
Axial force = 2 C,= axial force coefficient 
: VS oC 
a ee 2 ie ref NA Cy4= normal force coefficient 
2 pd 
5 OV Srey [22 Cups Cyp4= Magnus force coefficient 


Magnus force = 
Pitch damping moment = 


1 


4 Cy: Coefficient dueto q, = q +?’ total transverse angular velocity 
4 
sev" S,,d | ; \(Cu +Cy;) 


d 
Cy :Coefficient due toa, = a rate of changeof total yawangle 


Roll damping moment = 


1 7 ; : 
5 pV? S., d ee C,, C,p= roll damping moment coefficient 


1 d : 
oye Si A [et Cyp,8ina, Cypy= Magnus moment coefficient 
Magnus moment = V 


Pitch or overturning moment = 


: 0 Vy? S. dC,,, sina, Cv = pitch or overturning moment coefficient 
List of symbols 

Cy zero-yaw axial aerodynamic coefficient 

Cy square yaw axial aerodynamic coefficient 

C, lift aerodynamic coefficient 

Crp roll damping aerodynamic coefficient 

Cuo pitch damping aerodynamic coefficient 

Cv pitch moment coefficient 
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Cypy Magnus moment coefficient 


Cyp4 Magnus force aerodynamic coefficient 

= ss rojectile position, in inertial space 

Kips Vigo Zig ee . 

Se see projectile velocity components, expressed in no-roll-frame 
NRF ? © NRF ? "" NRF 

fe projectile roll rate, expressed in no-roll-frame 

P yrr 

ee - projectile pitch and yaw rates expressed in no-roll-frame 

rr >" rr 

ry wy projectile pitch and yaw angles 
s 

é projectile roll angle 

I projectile inertia matrix 

Lyx Lyy, Izz diagonal components of the inertia matrix 

Lyy, Lyz, Lyz off-diagonal components of the inertia matrix 

V total aerodynamic velocity 

p atmospheric density 

Sref projectile reference area (zd2/4) 

m mass of projectile 

t time 

a, aerodynamic angles of attack and sideslip 

R vector from the projectile center of mass to the center of pressure 
®MAC 

R vector from the projectile center of mass to the Magnus center of pressure 
MAX 

D projectile reference diameter 

g gravity 
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Atmospheric Flight Dynamic Simulation Modeling 
Of Spin-Stabilized Projectiles And Small Bullets With Constant Aerodynamic 
Coefficients 


Dimitrios N. Gkritzapis,* EliasE. Panagiotopoulos,” Dionissios P.Margaris® and Dimitrios G.Papanikas” 


* Laboratory of Firearms and Tool Marks Section, Criminal Investigation Division, 
First Lieutenant of Hellenic Police, Postgraduate Student of Fluid Mechanics Laboratory (FML), 
University of Patras, GR-26500, Patras, Greece 
Email: gritzap@yahoo.gr, web page: http://www.mech.upatras.gr/~fml 
> Fluid Mechanics Laboratory (FML), Mechanical Engineering and Aeronautics Department 
University of Patras, GR-26500, Patras, Greece 


Abstract. 4 full six degrees of freedom (6-DOF) simulation flight dynamics model is applied for the 
accurate prediction of short and long range trajectories of high and low spin-stabilized projectiles and 
small bullets via atmospheric flight to final impact point. The projectile is assumed to be both rigid (non- 
flexible) and rotationally symmetric about its spin axis launched at low and high pitch angles. The projectile 
maneuvering motion depends on the most significant forces and moments, in addition to wind, gravity and 
Magnus effects. The computational flight analysis is based on appropriate constant mean values of the 
aerodynamic projectile coefficients taken from an official tabulated database. The newer ICAO atmospheric 
model simulates the height distributions of density, pressure and temperature properties of the air. Static 
stability, also called gyroscopic stability, is examined as a necessary condition for stable flight motion in 
order to determine the sufficient initial spinning projectile rotation. The efficiency of the method developed 
gives satisfactory results compared with published data of verified experiments and computational codes 
on atmospheric dynamics model flight analysis. 


Keywords: Trajectory dynamics simulation, constant aerodynamic analysis, high and low pitch angles, 


Magnus effects, symmetric projectiles, static stability criteria 


Introduction 


Ballistics is the science that deals with the motion 
of projectiles. The word ballistics was derived 
from the Latin ‘ballista’, which was an ancient 
machine designed to hurl a javelin. The modern 
science of exterior ballistics! has evolved as a 
specialized branch of the dynamics of rigid bodies, 
moving under the influence of gravitational and 
aerodynamic forces and moments. 


Pioneering English ballisticians Fowler, Gallop, 
Lock, and Richmond?’ constructed the first 
rigid six degrees of freedom projectile exterior 
ballistics model. Various authors have extended 
this projectile model for lateral force impulses,*“ 
dual-spin projectiles,” ete. 


The present work addresses a full six degrees 
of freedom (6-DOF) projectile flight dynamics 
analysis for the accurate prediction of short and 
long trajectories of high spin-stabilized projectiles 
and bullets. The proposed flight model takes into 
consideration the influence of the most significant 
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forces and moments, based on appropriate constant 
mean values of the aerodynamic coefficients, 
in addition to wind and Magnus effects. The 
efficiency of the computational method developed 
gives satisfactory results compared with published 
data of verified experiments and computational 
codes on projectile trajectory analysis with various 
initial flight conditions at the firing site. 


Projectile model 


The present analysis considers two different types 
of representative flight projectile vehicles. A 
typical formation of the cartridge 105 mm HE M1 
projectile is presented in Figure | and is used with 
various 105 mm Howitzers such as M103 with 
M108 cannon, M137 with M102 cannon as well 
as NATO L14 MODS6 LS. Cartridge 105 mm 
HE M1 is of semi-fixed type ammunition, using 
adjustable propelling charges in order to achieve 
desirable ranges. The projectile producing both 
fragmentation and blast effects can be used against 
personnel and materials targets. 
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Figure 1. 105 mm HE M1 high explosive 
projectile artillery ammunition for howitzers. 


2,830 IN, 
(71.9 MM) 


Figure 2. 7.62 mm M852 match ammunition with 
a diameter of .30 caliber representative bullet. 


Also a 0.30 caliber (0.308 inch diameter), 
168 grain (=10.9 g) Sierra International bullet used 
in a National Match M14 rifle is loaded into 7.62 
mm M852 match ammunition for high power rifle 
competition shooting, as shown in Figure 2. The 
cartridge is intended and specifically prepared for 
use in those weapons designed as competitive rifles 
and also for marksmanship training. This bullet is 
not for combat use. The cartridge is identified by 
the cartridge case head stamping of MATCH. It 
also has a knurl at the base of the cartridge case 
and a hollow point boat-tail bullet. 


Table 1. Physical and geometrical data of 105 mm 
big projectile and 7.62 mm small bullet. 


Characteristics eee on 
M1 projectile bullet 

Reference diameter, 104.8 762 

mm 

Total length/nm 494.7 71.88 

Weight/kg 15.00 0.385 

Axial moment of 2326-102 7.2282-10° 

inertia/kg m 

Veversemomemt sailed’ ~S3787-407 


of inertia/kg m7 


Center of gravity 
from the base/mm 183.4 12.03 
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Physical and geometric characteristics data of the 
above mentioned 105 mm HE M1 projectile and 
7.62 mm bullet are illustrated in Table 1. 


Trajectory flight simulation model 


A six degrees of freedom rigid-projectile and 
bullet model’'° has three rotations and three 
translations. The three translation components 
(x, y, z) describing the position of the projectile’s 
center of mass and the three Euler angles (¢, 9, 


Figure 3. No-roll (moving) and fixed (inertial) 
coordinate systems for the projectile trajectory 
analysis. 


y) describing the orientation of the projectile body 
with respect to Figure 3. 


Two mean coordinate systems are used for the 
computational approach of an atmospheric flight 
motion. The one is a plane fixed (inertial frame, 
if) at the firing site. The other is a no-roll rotating 
coordinate system on the projectile body (no-roll- 
frame, NRF, ¢ = 0) with the Xypp axis along the 
projectile axis of symmetry and Yyrp, Zypp AXeS 
oriented so as to complete a right-hand orthogonal 
system. 


Newton’s laws ofthe motion state that rate of change 
of linear momentum must equal the sum of all the 
externally applied forces and the rate of change 
of angular momentum must equal the sum of all 
the externally applied moments, respectively. The 
force acting on the projectile comprises the weight, 
the aerodynamic force and the Magnus force. 
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The moment acting on the projectile comprises 
the moment due to the standard aerodynamic 
force, the Magnus aerodynamic moment and the 
unsteady aerodynamic moment. 


Therefore, the twelve state variables x, y, z, ¢, 0, 
y, U, V, W, P, g and r are necessary to describe the 
position, flight direction and velocity at every point 


of the projectile’s atmospheric flight trajectory. 
Introducing the components of the acting forces 
and moments with regard to the no-roll-frame (~) 
rotating coordinate system with the dimensionless 
arc length as an independent variable, we derive the 
following equations of motion for six-dimensional 
flight: 


; DD. _... DD 
x = Dcosw cos 6 — —siny) ¥ + weos w sin d — (1) 
V V 
' 2 TD es acti 
y = ee eee ee (2) 
D 
z' =—Dsin@ + — weos@ (3) 
V 
D D 
gy! = — pt+—tan|6]7 (4) 
V V 
D. 
o=— q (5) 
V 
Ds. 
= 7 (6) 
V cosé 
D D D 
Pao en" AP eb ave oD pve oa 4a ee (7) 
V 8m x 8m m x V V 
Tv 
:= DY © Cua? iw) +D —— PeCy py pwtan|6| — DF (8) 
m 
j Tv a. 4 7 | 2 D_. 
w =—gcosd—D — P Cyg (WW) —D Sag NE FOTO Ry (9) 
. 5 
p'=D beCyp (10) 
~p 3 Tv a és 
D p Cy (Wy )LEqicp + D an yea? | Lec + 
(11) 
z D Aye . 2 
+ D° —— pCyg 4+ D’ —— pCyy -—F *~ tan |o| 
l6lyy V I 
~/ 
= —D° —— p Cy (Fw) LEycp + D Pp a \LEucu + 
(12) 
i pCygF-D pCypg +— ba + — Gr tanp| 
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The projectile dynamics trajectory model 
consists of twelve highly first order non-linear 
ordinary differential equations, which are 
solved simultaneously by resorting to numerical 
integration using the fourth order Runge-Kutta 
method. In these equations, the following sets of 
simplifications are employed: the aerodynamic 
angles of attack a and sideslip £ are small a= w/V, 
B= V/V, the projectile is geometrically symmetrical 
Ly = Tyz = Lyz = 0, [py = Izz and aerodynamically 
symmetric. With the afore-mentioned assumptions 
the expressions of the distance from the center 
of mass to both the standard aerodynamic and 
Magnus centers of pressure are simplified. 


Aerodynamic model 


For the projectile trajectory analysis, a constant 
flight dynamic model is proposed for the 
examined test cases. The above calculations are 
based on appropriate constant mean values of the 
experimental average aerodynamic coefficients 
variations taken from an official tabulated 
database’! as shown in Table 2. 


Initial spin rate estimation 
In order to have a statically stable flight projectile 


trajectory motion, the initial spinrate P ) prediction 
at the gun muzzle in the firing site is very important. 
According to McCoy’s definitions,’ the following 
form is used: 


PD, = 2nV,,/7 D (rad/s) (13) 


where V, is the initial firing velocity (m s'), 7 
the rifling twist rate at the gun muzzle (calibers 
per turn), and D the reference diameter of the 
projectile type (m). Typical values of rifling 7 
are 1/18 calibers per turn for big projectiles and 


12 inches per turn for small bullets, respectively. 


Static or gyroscopic stability 


Any spinning object will have gyroscopic 
properties. In a spin stabilized projectile, the center 
of pressure, the point at which the resultant air 
force is applied, is located in front of the center of 
gravity. Hence, as the projectile leaves the muzzle 
it experiences an overturning movement caused 
by air forces acting about the center of mass. It 
must be kept in mind that the forces are attempting 
to raise the projectile’s axis of rotation. 


In Figure 4 two cases of static stability are 
demonstrated: In the top figure, CP lies behind 
the CG so that a clockwise (restoring) moment is 
produced. This case tends to reduce the yaw angle 
and return the body to its trajectory, therefore is 
statically stable. Conversely, the lower figure, 
with CP ahead of CG, produces an anti-clockwise 
(overturning) moment which increases further and 
is therefore statically unstable. It also possible 
to have a neutral case in which CP and CG are 
coincident whereby no moment is produced. 


There is clearly an important correspondence in 
the distance between the center of pressure and 
the center of gravity and the static stability of the 
round. This distance is called the static margin. By 
definition, it is positive for positive static stability, 
zero for neutral stability and negative for negative 
stability. 


Classical exterior ballistics'' defines the gyroscopic 
stability factor S, in the following generalized 
form: 


T 2>~2 
55 ee (14) 
2 pI, SDVC,, 


Table 2. Aerodynamic parameters of atmospheric flight dynamic model. 


Aerodynamic coefficients 


105 mm HE M1 projectile 


7.62 mm bullet 


Drag Cp 0.243 0.235 
Lift C, 1.76 2.205 
Roll damping Crp —0.0108 —0.01 
Pitch damping Cw —9.300 —4.7 
Overturning moment Gus 3.76 2.92 
Magnus moment Cypy —0.381 —1.26 
Magnus force C 0.215 0 


Fa — a) 7 ee 
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Figure 4. Static stability/instability conditions. 


This may be rearranged into: 


Fltelte] 


Equation (15) above shows that the static factor is 
proportional to the product of four terms, depending 
on the geometric technical characteristics of 
projectile shape model, the square axial spin 
to velocity ratio, the aerodynamic overturning 
moment coefficient and the proposed atmospheric 


21. 
IrD 


Ss = 


& 


density model. 


Computational simulation 


The flight dynamic model of 105 mm HE M1 and 
7.62 mm projectile types involves the solution of 
the set of the twelve non-linear first order ordinary 
differential, equations (1)-(12), which are 
solved simultaneously by resorting to numerical 
integration using a fourth order Runge-Kutta 
method. 


Initial flight conditions for both dynamic trajectory 
projectile models with constant aerodynamic 
coefficients are illustrated in Table 3. 


Results and discussion 


The flight path trajectory motions with constant 
aerodynamic coefficients of the big 105 mm 
projectile with initial firing velocity of 494 ms", 
initial yaw angle 3°, rifling twist 1/18 caliber 
per turn and initial yaw rates 3.61 rad s' and 
3.64 rad s' at 45° and 70°, respectively, are 
indicated in Figures 5 and 6. The calculated impact 
points of the above no-wind trajectories with 
the proposed constant aerodynamic coefficients 
compared with accurately estimations of McCoy 
flight trajectory analysis' provide basic differences 
for the main part of the atmospheric flight motion 
for the same initial flight conditions. Comparative 
computed trajectories of the 105 mm projectile 
with a 5.0 ms | mean crosswind blowing are also 
indicated in the above diagrams. 


For the 105 mm M1 projectile, fired at 45° 
at sea-level neglecting wind conditions, the 


Table 3. Initial flight parameters of the projectile examined test cases. 


Initial flight data 105 mm HE M1 projectile 
x/m 0.0 

y/m 0.0 

z/m 0.0 

¢/deg 0.0 

A/deg 45 and 70 
y/deg 3.0 

u/ms ' 494.0 

vim s! 0.0 

wim s' 0.0 

pirad s ' 1644.0 

qirad s '' 3.61 and 3.64 
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7.62 mm bullet 


0.0 

0.0 

0.0 

0.0 

0.84 and 32 
2.0 
792.48 
0.0 

0.0 
16335.0 
25.0 
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—— Present Analysis Pitch Angle 45 deg 
—— = McCoy Model Pitch Angle 45 deg 
4000 -———— —_ Wind Model Pitch Angle 45 deg 


altitude, m 


+ 1 : T : T : T + T : 
0 2000 4000 6000 8000 10000 12000 


range, m 


Figure 5. /mpact points with constant 
aerodynamic coefficients for 105 mm projectile at 
quadrant elevation angle of 45 degrees. 


predicted range to impact of the McCoy model is 
approximately 11500 m and the maximum height 
is 3490 m, as shown in Figure 5. The corresponding 
maximum range of the wind trajectory (green 
solid line) is approximately 11510 m with 3780 m 
maximum height. From the computational results 
of our presented analysis at 45°, the predicted range 
to impact is almost 11550 m and the maximum 
height is slightly over 3795 m. 


Also in Figure 6 the predicted level-ground range 
of the McCoy model is 7300 m at 70° with a 
maximum height at about 6 km. For the crosswind 
trajectory estimation the corresponding values are 
5400 m and 6350 m, respectively. The presented 
method gives overestimated values of range 
impact and flight height calculations. 


— Present Analysis Pitch Angle 70deg 
—— McCoy Model Pitch Angle 7Odeg 
— Wind Model Pitch Angle 70deg 


altitude 


to 0 
i) 1000 2000 3000 4000 5000 6000 7000 8000 


range, m 


Figure 6. Flight path trajectories with constant 
aerodynamic coefficients at quadrant elevation 
angle of 70 degrees for 105 mm projectile. 
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—— Present Analysis Pitch Angle 0.84 deg 
| — Wind Trajectory Pitch Angle 0.84 deg 


altitude, m 


r 7 T T T ra r T 
tt) 200 400 600 800 1000 1200 


range, m 


Figure 7. Constant aerodynamic coefficients 
impact points for 7.62 mm bullet. 


A small bullet of 7.62 mm diameter is also 
examined for its atmospheric constant flight 
trajectory predictions in Figures 7 and 8 at low 
and high pitch angles 0.84° and 32°, respectively, 
with initial firing velocity of 793 ms‘, initial 
yaw angle 2°, yaw rate 25 rads ' and rifling twist 
12 inches per turn. The impact points of the above 
trajectories are compared with a corresponding 
5.0ms ' mean crosswind blowing motion and an 
accurate flight path prediction with Nennstiel’s 
trajectory analysis'! for a cartridge 7.62 mm ball 
M80 bullet type with initial firing velocity of 
838 ms |. The main differences are presented at 
high altitudes from the firing site sea-level. 


At 0.84° the 7.62 mm bullet fired at sea-level has a 
range with the wind model of almost 1090 m (red 


1400 —— Present Analysis 7.62mm Pitch Angle 32deg 
| ——_Nennstiel Model 7.62mm Pitch Angle 32deq 


altitude, m 


T T r T 
tt) 1000 2000 3000 4000 


range, m 


Figure 8. Constant flight atmospheric trajectory 
analysis for small bullet compared to Nennstiel 
prediction computational algorithm. 
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— Present Analysis Pitch Angle 45 deg 
—— Present Analysis Pitch Angle 70 deg 


a 
n 


angle of attack, deg 
3 
1 


i 
8000 10000 12000 


r + T ~ 
0 2000 4000 6000 


range, m 


Figure 9. Comparative angle of attack on 
the range of 105 mm projectile with constant 
coefficients at 45° and 70°. 


solid line). The presented method for the same 
initial pitch angle gives a corresponding value 
of 1105 m, and the height is slightly over 5.1 m. 
At 32 degrees initial pitch angle, the predicted 
range of the 7.62 mm bullet to the impact point 
of Nennstiel’s analysis’! is approximately 4100 m, 
and the maximum height is 1170 m. On the other 
hand, the proposed computational no-wind flight 
analysis gives an impact point at almost 4 km and 
a maximum height at about 1.3 km. 


The strong characteristics alterations of the total 
angle of attack influence the distributions of the 
most basic projectile trajectory phenomena taking 
into account constant aerodynamic coefficients 
during the whole atmospheric flight motion. Its 
effects are indicated in Figure 9 for the 105 mm 


— Present Analysis Pitch Angle 45deg 
—— Present Analysis Pitch Angle 70deq 


r t + T > T —+ r + r 
tt) 2000 4000 6000 8000 10000 12000 


range, m 


Figure 10. Comparative static stability variation 
with constant aerodynamic model at low and high 
quadrant angles for the 105 mm projectile 
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— Present Analysis Pitch Angle 32 deg 
Present Analysis Pitch Angle 0.84 deg 


angle of attack, deg 
3 


+ : T T + 
0 1000 2000 3000 4000 


range, m 


Figure 11. Angle of attack versus range at pitch 
angles of 0.84° and 32° for 7.62 mm bullet. 


M1 projectile type fired from a 1/18 twist cannon 
with a muzzle velocity of 494 m s' at quadrant 
elevation angles of 45° and 70°, respectively. 


After the damping of the initial transient motion, at 
apogee, the stability factor for 45° increases from 
3.1 at the muzzle to 23, and then decreases to a 
value of 8 at the final impact point, as presented in 
Figure 10. The corresponding flight behavior at 70° 
initial pitch angle shows that the transient motion 
damps out quickly and the yaw of repose grows 
nearly 14° at the apogee, where the gyroscopic 
static stability factor has increased from 3.1 to 121 
and then decreased to a value of almost 6.6 at the 
impact area. 


Figure 11 shows the angle of attack variations 
with range for the 7.62 mm bullet, fired from a 


—— Present Analysis 
Pitch Angle 32 deg |” 


+ r ~ T r T + 
tt) 1000 2000 3000 4000 


range, m 


Figure 12. Static stability versus range at pitch 
angle of 32°, for 7.62mm bullet. 
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12 inch twist per turn with a muzzle velocity of 
793 ms‘ at quadrant elevation angles of 0.84° 
and 32°, respectively. The transient motion damps 
out very quickly for the two cases. The gyroscopic 
stability factor for 32° was 1.5 at the muzzle, 
increased to 54 at the summit of the trajectory, and 
then decreased to a value of 24 at the impact point, 
as indicated in Figure 12. 


The estimation of static stability factors, 
equation (14), for the examined projectile test cases 
is based only on the square axial spin to velocity 
ratio due to the constant overturning aerodynamic 
moment. The magnitude of this suggests that 
if a projectile or bullet is statically stable at the 
muzzle, it will be statically stable for the rest of 
its atmospheric flight motion. As the rotational 
velocity is much less damped than the transversal 
velocity (which is damped due to the action of 
the drag force), the static factor S, increases at 
least for the major part of the predicted projectile 
trajectory. 


Conclusions 


The six degrees of freedom (6-DOF) simulation 
flight dynamics model is applied for the accurate 
prediction of short and long range trajectories 
of high and low spin-stabilized projectiles and 
small bullets. The parallel computational analysis 
takes into consideration the effects of constant 
aerodynamic coefficients. The criteria and 
analysis of gyroscopic stability are also examined. 
The computational results of the dynamic flight 
trajectories with constant aerodynamic force 
and moment coefficients are in good agreement 
compared with other technical data and recognized 
projectile flight dynamic models. 
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5th Workshop on Pyrotechnic Combustion Mechanisms, 
held October 6th 2007 in Beane, France 


Ernst-Christian Koch* and Rutger Webb” 
a NATO Munitions Safety Information Analysis Center, Brussels, Belgium 
b TNO, Defence Safety & Security, Rijswijk, The Netherlands 


For the fifth time the “Workshop on Pyrotechnic Combustion Mechanisms” was held and organized in 
conjunction with the International Pyrotechnics Seminar. This time it was the 34th IPS Seminar in conjunction 
with EUROPYRO 2007, in Beaune, France. 


This workshop had a total of 21 participants, of which 4 were from France, 4 from Germany, 4 from USA, 
3 from the UK, 2 from The Netherlands, I from South Africa, 1 from Japan, 1 from Russia, and 1 from 
Finland. 


This event received very positive feedback from both participants and presenters. 
Papers based on the following presentations have already, or are expected to be published elsewhere. 


¢ Volker Weiser, Fraunhofer ICT, Germany, Modeling Spectral Emission and Radiation Intensity of 
Pyrotechnic Reactions 


¢ Alexander Dologoborodov, Semenov Institute Moscow, Russia, Mechanoactivated Energetic Composites 
on the Base of Metal - Oxidizer Mixtures 


¢ Selena Burn, BAE Systems, UK, Personal Protective Equipment 


The next workshop on pyrotechnic combustion mechanisms will be held on July 13 2008 in Fort Collins 
Colorado in conjunction to the 35th International Pyrotechnics Seminar. The topics then will focus on safety 
and sensitivity of pyrotechnics. For further information on that event please refer to www.pyroworkshop.net 
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Ignition and Combustion of Aluminium in High Explosives 


Arno Hahma 
University of Jyvaskyla, Department of Chemistry, Finland 


Abstract: /gnition and combustion properties of aluminium powder up to 250 kbar pressure were successfully 
determined. Hydrodynamic modelling combined with parametrized reaction profiles was used to match the 
experimental metal plate profiles to calculated values and the metal reactivity was determined from the 
parameters thus obtained. The reaction profiles over a wide range of specific volumes and pressures were 
obtained. 


Aluminium ignited only if the detonation temperature was higher than 2300 K and the aluminium oxide was 
decomposed, even a strong shock wave was not able to clean the aluminium surface to promote ignition 
at lower temperatures. Aluminium reactivity at the Chapman—Jouguet plane (C-J plane) was concluded 
to be insignificant in all cases studied. A maximum of only 17 % reactivity by an expansion ratio of 1 : 50 
was reached. The aluminium combustion could be best described with a power law Ap" and the pressure 
exponent of aluminium was found to be slightly positive: +0.33. Nitrogen rich explosives were concluded 
to be necessary to convert Al energy effectively to mechanical work and to unambiguously measure the Al 


reactivity. 


Keywords: Aluminium, Chapman—Jouguet plane, ignition, combustion, hydrodynamic modelling 


Introduction 


In order to measure the combustion characteristics 
of aluminium and other metal powders over the 
entire pressure range of interest, at least three 
different methods need to be used. Combustion 
at pressures ranging from subatmospheric up to 
200-300 bar can be studied with a windowed 
pressure vessel, where the combustion can be 
observed with a fast camera.' Optionally, the 
powder may be dropped into a preheated gas”? or 
electrostatically levitated* to ensure the observed 
particles are not agglomerated. 


The pressure range from approximately 200 bar 
up to 4000 bar can be studied with a gun 
experiment, where the metal powder is mixed with 
a propellant and the pressure—time history of the 
gun is analyzed. This method is accurate and can 
produce combustion rate data over a wide pressure 
range yielding a two parameter burn equation for 
the metal powder. Details of this method are not 
within the scope of this paper, but will be published 
elsewhere. 


For pressures higher than 4000 bar explosives 
are needed to generate the pressure, temperature 
(ignition source) and the gaseous environment, 
where the metal powder burns. Ideally, the 
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explosive would use the aluminium powder to 
be studied as the only source of energy to the 
detonation, but this has been shown not to work.** 
Aluminium cannot support a true detonation even 
with the finest available aluminium powders due 
to strong gas volume reduction by the reaction 
Al + oxidizing gas — AlI,O3. In addition, the 
covolumes of gaseous Al,O,; and other Al,O,H, 
species are much smaller than that of any gaseous 
starting materials separately leading to further 
reduction of the product specific volume even to 
less than that of the starting materials. Both effects 
prevent a detonation from occurring, even though 
there is plenty of energy available to drive it. 


Direct pressure measurement from flash powder 
mixtures is too unreliable to produce useful dataand 
cannot be used to track the reactions. In addition, 
such mixtures do not reach as high pressures as 
explosives do. As a result, flash powder or similar 
mixtures are not an option to study aluminium 
reactions at high pressures. 


This study also addresses the role of Al under 
detonation, which has been a_ subject of 
controversy. It has been claimed Al reacts fully and 
contributes to the shock energy of a detonation” 
while claims about its inertness® have also been 
made. Measurements of the detonation velocities 
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of fuel-oxidizer solutions with and without 
Al give rise to different interpretations, since 
perchlorate solutions (ethylene glycol/lithium 
perchlorate) show an increase in their detonation 
velocity as a function of aluminium content up to 
the stoichiometric ratio.” This increase is due to 
the increased density of the mixture rather than 
aluminium combustion at the Chapman—Jouguet 
(C-J) plane easily leading to false conclusions about 
Al reactivity. Similarly, ammonium dinitramide 
(ADN) increases its detonation velocity if nano- 
sized aluminium is mixed in.'? ADN is a highly 
non-ideal explosive that does not burn completely 
at the C-J plane and thus its detonation velocity 
is far less than theoretically calculated. Any 
augmentation of the intrinsic reactions of ADN at 
the C-J plane may therefore enhance the detonation 
velocity, such as surface catalysis on the nano- 
sized aluminium particles or slight reactions of 
the aluminium increasing local temperature, such 
as the release of the intrinsic lattice energy of the 
electroexploded submicron aluminium used in the 
study. '!"! 


Morerecently ithas been shown!*'’ that aluminium 
does not react significantly at the C-J plane or 
in its vicinity. As a result, Al cannot contribute 
significantly to the shock wave energy due to 
too slow reactions, but can ignite and burn while 
the product gases expand.'**! Aluminium can, 
however, greatly enhance the metal acceleration 
performance of an explosive'®*' as shown even 
in this study, provided that the metal burns fast 
enough before the gases have expanded beyond 
approximately 20-30 times their original volume, 
and that there is enough working fluid available to 
convert thermal to mechanical energy. 


Prerequisites 


To study the reactivity of aluminium in explosives, 
the explosive has to be selected carefully and 
not just any high explosive will work. Many 
common explosives yield almost no difference 
to both detonation velocity and the mechanical 
energy as a function of aluminium reactivity by 
the expansion of approximately | : 50. The metal 
content does contribute to the blast energy of the 
explosive, if the metal powder ignites and burns 
within the expanding gas cloud generated or when 
mixed with surrounding materials. This energy, 
however, is released at very low pressure and is 
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not available to study the reactivity of aluminium 
at high pressures nor to accelerate the metal. 


Consequently, it is difficult to determine the 
reactivity of aluminium by using explosives, 
where there are two or more possible calculated 
reactivities at the C-J plane and at the following 
reaction zone, all fitting experimental data. It 
is even more difficult to judge the reactivity by 
measuring detonation velocity only, since it may 
or may not change as a function of aluminium 
reactivity and there may be several reactivities 
fitting the same detonation velocity even for 
nitrogen rich explosives (Figure 1). 


If the explosive produces enough gas relatively 
inert towards burning aluminium, e.g. nitrogen, 
differences in mechanical energy and even in 
detonation velocity will be seen, since there is 
enough working fluid to compensate for the gas 
volume reduction by Al combustion. 


In order to study aluminium reactivity, the 
explosive has to be rich in nitrogen and have as 
low a carbon content as possible to avoid carbon 
coagulation on the aluminium particles quenching 
combustion. Low carbon content also makes the 
explosive gases less luminous and makes it more 
possible to employ photographic techniques for 
recording the results. 


The amount of aluminium was selected such that 
a full combustion yielded a maximal increase in 
the mechanical energy delivered by the explosive, 
since mechanical energy could be measured and 
served as an indicator of the degree of aluminium 
combustion. In practice, this translated into 
mixing 80% of the stoichiometric amount of 
aluminium calculated from converting all oxygen 
in the explosive to aluminium oxide, which means 
every explosive in this study contained a different 
amount of aluminium mixed into it ranging from 
22% up to 26%. Ideally, the detonation velocity 
should also respond to aluminium reactivity, but 
this cannot be required, since this property is not 
systematic but differs with each explosive. The 
nitrogen content of the explosive has to be above 
approximately 35% in order to generate enough 
working fluid to make the aluminium reactivity 
clearly visible in the mechanical energy delivered 
to metal acceleration and the carbon content should 
be as low as possible. Naturally, the explosive 
also has to contain enough oxygen to burn the 
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D vs. A(Al) for NTO/AI; renormalized BKW parameters and different covolume sets, p = 1.85 g/cm 
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Figure 1. Detonation velocity as a function of aluminum reactivity for the system NTO/AL 
76.5/23.5, p = 1850 kg m®, six different covolume sets for the product gases. The solid line is 


measured detonation velocity. 


aluminium. 


In addition, it is important to use a covolume set 
that contains no mixed theoretical and experimental 
components, if the Becker—Kistiakowsky—Wilson 
(BKW) equation’ of state is used to model the 
product gas thermodynamic functions. Mixed 
covolumes tend to cause even more complication 
in interpreting the results, since several different 
reactivities may match the calculation and 
experiment (Figure 1) with any explosive. By 
mixing theoretical and experimental covolumes 
the dataset may become nonphysical, or at 
least not consistent with the a priori theoretical 
assumptions. In this work, theoretical covolumes 
were calculated by ab initio molecular modelling 
and scaled against solid nitrogen shock Hugoniot 
data.” 


Additional information can be obtained if several 
explosives are used and they are selected such 
that the elemental composition of the explosives 
remains the same while the heat of formation varies 
producing different detonation temperatures. In 
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addition, the explosives have to be processable, 
i.e. insensitive enough for safe handling and for 
pressing metallized charges close to theoretical 
mean density (TMD). 


The above criteria cannot be fully met as far as the 
elemental composition is concerned, but a close 
approximation is reached with the following set 
in an order of increasing detonation temperature: 


3-nitro-1,2,4-triazol-5-one (NTO), nitroimino- 
1-nitrohexahydro-1,3,5-triazine (NNHHT), 
1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX), 
1,3,5-trinitrohexahydro-1,3,5-triazine (RDX) 


and hexanitrohexa-azaisowurtzitane (HNIW, 
CL-20). This set was found through extensive 
screening of the calculated properties of most 
known explosives with Al and by comparing the 
mechanical energy of the explosive at 0% and 
100% Al reactivity at 1 : 20 expansion trying to 
maximize the difference. 
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v-p-isentropes for NNHHT/AI as a function of Al reactivity 0..87 % 
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Figure 2. Pressure vs. specific volume for NNHHT/AI 78/22 as a function of Al reactivity. The 
lowest curve is the 0% reactive isentrope and the uppermost is the 87% reactive isentrope with 
immediate reaction at the C-J plane. Interval between curves: Ad = 20%. 
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Figure 3. Temperature vs. specific volume for NNHHT/AI 78/22 as a function of Al reactivity. The 
lowest curve is the 0% reactive isentrope and the uppermost is the 87% reactive isentrope with 
immediate reaction at the C-J plane. Interval between curves: Ad = 20 %. 
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Theoretical 


When the explosive is chosen correctly, the 
combustion of aluminium affects the expansion 
of the gases strongly and increases the mechanical 
energy available. Al combustion slows down 
the pressure drop during gas expansion and the 
additional energy from Al becomes available at 
lower pressures (Figure 2). At the high pressure 
region, there is only a small difference in the 
pressure as a function of Al reactivity. Most of 
the energy is released while the gases expand to 
approximately 1 : 20 expansion ratio. The band 
of pressure vs. specific volume isentropes (v—p 
isentropes) for most other explosives is much 
narrower than that of the sample explosive 
NNHHT, which is one of the most nitrogen rich, 
insensitive CHNO-explosives available. The 100% 
reactivity isentrope may even lie below that of 
the 0% reactivity at some range of pressures with 
low nitrogen explosives, such as trinitrotoluene 
(TNT). 


The theoretical temperature vs. specific volume 
isentropes (v—T isentropes) are depicted in 
Figure 3 as a comparison of how temperature 
behaves as a function of Al reactivity. As the 
pressure changes relatively little, the energy has to 
stay in thermal form. Therefore, temperature has 
to rise significantly and the effect of aluminium 
combustion is best visible on the v—T isentrope. 


The NNHHT/AI mixture is used as an example 
in the figures, since NNHHT yields the largest 
differences in the wp space due to its high 
nitrogen content and therefore clearer figures can 
be made. With NTO and RDX, the figures look 
similar, but the band of isentropes is somewhat 
narrower. 87% maximum reactivity was used, 
because calculations for any larger reactivity 
with NNHHT would not converge. However, this 
limitation causes no harm to interpretation of the 
measurements, since NNHHT did not ignite the Al 
powder. 


Burn models 


Determining aluminium burn characteristics in an 
explosive is a typical reverse problem, which is 
impossible to solve without a priori information 
about the aluminium burn characteristics. As this 
information is not available initially, it has to be 
estimated and used to model the forward problem 
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instead. By feedback from forward simulations, 
the initial estimate is adjusted until the simulation 
matches the experiments. The adjustment is done 
through one parameter, which adjusts the burn rate 
vs. specific volume function for the aluminium 
powder. 


Four different burn models were used: immediate, 
linear, logarithmic and exponential. The three 
latter models scale the burn degree 4 of the 
aluminium powder against the specific volume, 
weighted with the adjustment parameter 4 called 
the slope of the burn function. The meaning of 
this parameter and the effect of the different burn 
models can be visualized in Figures 4—9 and is 
explained in more detail in the figure captions. v-T 
isentropes are used as examples for figure clarity. 
In this case, NTO/AI was picked as the sample 
system, since this combination yields the largest 
differences in temperature producing the clearest 
figures. The immediate burn corresponds to the 
partial reaction isentropes depicted in the same 
figures. The above elementary functions were 
chosen since any function can be approximated 
by a linear combination of the these three with 
far fewer terms than a polynomial would use, 
potentially opening the possibility of fitting the 
data to experiments even more accurately than 
was possible in this study. 


The paths corresponding to each function with 
each slope were calculated in the energy (£), 
pressure (p) and temperature (7) vs. specific 
volume planes and fourth degree polynomials 
were fitted to these points. The polynomials were 
used in the hydrodynamic modelling to describe 
the expansion of the gases. The software limited 
the function to being a fourth degree polynomial, 
since the hydrodynamic modelling code could 
handle polynomials up to the fourth degree only. 
Otherwise, a linear combination of exponential 
functions or that of different functions would have 
been a much better choice. 


The generated reaction paths yield very different 
types of profiles for 7, p and E, many of which 
are unrealistic and can be filtered out later, when 
the actual simulations are run and the results 
are compared to experiments. This filtering was 
carried out manually, since no automatic method 
was available to do the comparison with the 
recorded experimental data. 
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Linear v-’ -space for NTO/AI 
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Figure 4. Linear burn model paths for NTO/AI in the v-d plane. The slope is defined as the multiplier n in 
the equation 
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where n € N,n € [0,9]. This translates into the distance from the point (0 1) in the horizontal direction 
in the above figure, end point at (3 1). 


Linear burn model particle paths in the T-v-plane for NTO/AI 
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Figure 5. Linear burn model paths in the v-T plane for NTO/AI 76.5/23.5. The x-marks are the BKW 
calculated points and dash—dot lines are partial immediate reaction isentropes. Solid thin lines are fourth 
degree polynomial least squares fits to the BKW-calculated points. Starting point: (—1 7.6), end point: (3 
7.7). 
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Exponential v-’ —-space for NTO/AI 
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Figure 6. Exponential burn model paths for NTO/AI 76.5/23.5 in the v-A plane. The slope is defined as the 
multiplier k in the equation 
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A=k/3,k € Zk € [—10,10]. Note that k = 0 corresponds to the linear path with n = 9. 
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Figure 7. Exponential burn model paths in the v-T plane for NTO/AI 76.5/23.5. The x-marks are the BKW 
calculated points and dash—dot lines are partial immediate reaction isentropes. Solid thin lines are fourth 
degree polynomial least squares fits to the BKW calculated points. Starting point: (—1 7.6), end point: (3 
7.7). Only slopes from —3 to +4 are shown for clarity omitting k = 0, which is equal to n = 9 in Figure 5. 
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Logarithmic v-’ -space for NTO/AI 
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Figure 8. Logarithmic burn model paths in the v-T plane for NTO/AI 76.5/23.5. The slope is defined as 
the coefficient A in the equation 
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Curves with A © N, A € [1,9] are shown in the figure. A = 0 corresponds to the linear path with n = 9. 


Logarithmic burn model particle paths in the v-T-plane for NTO/AI 
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Figure 9. Exponential burn model paths in the v-T plane for NTO/AI 76.5/23.5. The x-marks are the BKW 
calculated points and dash—dot lines are partial immediate reaction isentropes. Solid thin lines are fourth 
degree polynomial least squares fits to the BKW calculated points. Starting point: (—1 7.6), end point: (3 
77, 
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Immediate burn model 


Immediate burn is the simplest, zeroth order 
approximation of Al combustion. 


Assuming part of the Al to react immediately 
at the C-J plane and the rest to be inert or the 
reactions would have frozen immediately after 
the shock front and survived through the entire 
process is not reasonable. If the Al ignites, then it 
most likely burns completely at some stage of the 
explosion, since the temperature will jump from 
Al combustion leading to an easier ignition of the 
rest of the Al. Thus, the reasonable assumptions 
are either no reaction, immediate full reaction or 
continuous burn while the gases expand. 


There is, however, one condition, which might 
realize an immediate burn up to a degree less 
than 100%. If carbon present in the explosive 
gases precipitates on the aluminium and quenches 
combustion, in theory it is possible to reach a 
condition where part of the Al has burned at 
the C-J plane and the reaction zone and then no 
combustion takes places thereafter. To check this 
possibility, immediate reactions were scanned as 
well starting from 0% up to 100% theoretical Al- 
reactivity in 1% steps. None of these calculations 
could match the results for RDX/AI, where the 
Al clearly did ignite, while 0% reactivity of 
course matched results perfectly for the two other 
systems, in which the Al did not ignite. Examples 
of the immediate burn model particle paths in 
the v—T plane are depicted as dot—dashed lines in 
Figures 5, 7 and 9. 


Linear burn model 


In this model, Al burns linearly as a function of 
In v up to a predetermined value of the specific 
volume, where the final burn degree is reached, 
1.€. 


r 


_ Inv | 3 
— l+4n 3+n 


where n © N, n © [0,9]. This model is also 
unrealistic, but represents a _ first order 
approximation. The results given by this model 
did not fit the experimental data for any explosive, 
but were closer than the immediate model for 
RDX/AI. Examples of the linear burn model and 
corresponding v—T isentropes are depicted in 
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Figures 4 and 5. 


Exponential model 


In this model, the particle burn degree was assumed 
to follow an exponential function 
eA ‘Inv —A 


x € 
eA — e-A eeA — e-A 


A=k/3,k © Z,k € [-10,10] as a function of In(v) 
in fact leading to a power law burn rate in the 2 
vs. v space by simplifying the above equation. In 
other words, the slope A in this model corresponds 
to a pressure exponent of the metal powder, since 
the specific volume is directly proportional to 
pressure in a close approximation. It turned out 
that the exponential burn model could match the 
results for RDX/AI with arbitrary accuracy and the 
pressure exponent of the metal combustion under 
extreme pressures was therefore obtained as a side 
product. Examples of the exponential burn model 
and the corresponding v—7 isentropes are depicted 
in Figures 6 and 7. 


Logarithmic model 


In the logarithmic model, the 2 vs. In(v) is given 
by In ((e 4 + In(v))/e “) and is normalized to limits 
[—1 3], e.g. zero at In(v) = —1 and unity at In(v) = 3 
with A € N,A €.'!° This leads to: 


_ in( ae) Ine) 


ne = i 


This model produces rather peculiar particle paths 
in the v-7, v—-p and v-E planes and is therefore 
unrealistic as can be seen in Figures 8 and 9. 
However, even this model did produces better 
results than immediate burn or linear burn for 
RDX/AI showing that the solution was closer to 
reality, but could not fit the data at all points. 


Modelling 


All detonation velocity and isentrope calculations 
were carried out with Fortran BKW7 software 
and using the BKW equation of state (EOS) for 
modelling the product gases. The BK W parameters 
were as follows: a: 0.55093, £: 0.13932, x: 15.631 
and 6: 377.64. These were obtained by fitting 
against experimental data from 140 different 
explosives including aluminized mixtures and 


Journal of Pyrotechnics, Issue 26, Winter 2007 


minimizing the root square mean value of the 
difference to calculated and observed detonation 
velocities and pressures. This way, the error 
marginal could be approximately halved compared 
to the original RDX parameters for the BKW 
EOS. The covolumes and BKW parameter fitting 
are the subject of another paper to be published 
elsewhere. 


The software TDL7 was used for most of the 
hydrodynamic calculations, except the three 
dimensional (3D) check runs, which were carried 
out with Dyna3D22. TDL is a two dimensional 
(2D) code, which turned out not to be a problem. 
The model was run as a true 3D-model once and 
the results were compared to a 2D-run by TDL. 
Differences at the edge location and geometry 
were only subtle at the end of the run (less than 2% 
difference in the position of the outermost point of 
the plates) and it was concluded that the 2D-model 
is more than adequate to model the problem. In 
addition, it ran approximately a thousand times 
faster than the 3D-model, which would have been 
totally impossible to apply to the entire study. 


Lagrangian meshing of the calculation domain 
was applied. There were many reasons to choose 
a Lagrangian meshing scheme: saving computer 
resources and completing the task in a reasonable 
time, no need to mesh the surroundings of the 
charge, better convergence than a Eulerian model 
with the software used and a non-distorted flow 
allowing the Lagrangian scheme to be used. Since 
the simulations were two-dimensional, there 
was also no need to allow the expanding gases 
go beyond the plate boundaries, which limited 
the amount of mesh distortion and allowed a 
Lagrangian scheme to be used. 


The Lagrangian meshing scheme proved the most 
robust and the entire length of the test charges 
could be modeled, while with Eulerian meshing 
the calculations usually crashed by arriving at 
half length of the charge. This was due to strong 
turbulence and to discontinuities of the geometry, 
which the Eulerian model attempted to solve, but 
failed, since the mesh cell size was many orders 
of magnitude above the Kolmogorov scale. It 
is impossible to run any flow dynamics at these 
flow rates using direct numerical simulations, but 
a turbulence model would have been necessary. 
As the codes used do not include one, they were 
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limited to rather short time spans only. Meshing 
was chosen to be rectilinear, as this geometry 
proved most robust and the calculations could be 
run furthest in time. The boundary types and mesh 
density are shown in Figure 10 and Table 1. 


In addition, the 3D-model could not be made 
to converge with aluminized explosives, but 
usually failed after the detonation had propagated 
approximately 50 mm to the charge. The exact 
reason for this behavior was not found, but 
the meshing was suspected to be too coarse for 
the 3D-model to handle non-conventional gas 
expansion schemes for the explosive leading 
to stronger turbulence. The architecture of the 
software (memory limitations) and computational 
resources available prevented increasing the mesh 
density beyond the density that was used in this 
study (Figure 10). 


Calculation blocks for the hydrodynamic problem 


<= Copper plate, 2 x 150 mm 
2 cells in X direction, 
240 cells in Y direction. 


100 


=< Model explosive, 0 cells in X, 239 in Y direction. 


Y/mm 


Initiating explosive, 20 cells in X, 1 in Y direction. 


(not to scale) 


Figure 10. Meshing blocks and number of cells 
for the hydrodynamic problem. All blocks have 

uniform rectilinear meshing with the same cell 

dimensions to match the boundaries. 
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Table 1. Boundary types for the different mesh blocks defined for the hydrodynamic simulation. The 


blocks are depicted in Figure 10. 


Block Left boundary Right boundary Top boundary Bottom boundary 
Initiating explosive Axial Material Material Free surface 
Model explosive Axial Material Continuum Material 

Copper plate Material Free surface Continuum Free 


The output data from TDL were analyzed 
and visualized with Application Visualization 
System (AVS) software on the Centre for 
Scientific Computing supercomputers (CSC, 
Espoo, Finland). AVS was also used to overlay 
the calculations with the experimental X-ray 
radiograph images. 


The shock wave was modeled with multiple shock 
Forest Fire (MS-FF) model. The FF constants were 
taken from the literature for RDX7 and determined 
for the other explosives. The shock model scales 
the explosive burn rate with pressure and allows 
multiple shock waves to pass through the possibly 
partially reacted explosive. Thus, partially reacted 
explosive can be present, shock reflections may 
change the detonation front structure and transverse 
detonation wave components may be present. 


The Forest Fire constants were determined for 
each explosive mixture. The distance to detonation 
as a function of initiating shock pressure was 
measured and the data were fitted to the Forest 
Fire model using a 12th degree polynomial fit. The 
measurements were carried out as card gap tests 
and the location of detonation start was determined 
with short circuit probes and an oscilloscope, 1.e. 
the acceleration distance for the explosive was 
measured for each shock strength. Acrylic sheet 
was used as the shock dampener between the 
donor charge (RDX : wax 95 : 5) and the acceptor 


charge. Shock pressures were taken from tables” 
as a function of acrylic sheet thickness. 


The end point in v-A space in the y-direction was 
not varied, since not much energy conversion is 
taking place at the low pressure end. Therefore, 
it is not critical whether the end point of the 
Al reaction path is at In(v) = 2 or In(v) = 3 or 
somewhere in between. Since the expansion ratio 
in the experiments at the end of the walls could not 
be determined accurately, it was decided to fix the 
specific volume at the end of Al-burn to In(v) = 3. 


The gas equation of state parameters (polynomial 
coefficients determined by the Fortran BKW 
-runs for each reaction profile) were varied and 
the calculations were repeated for each set of 
coefficients. The final reactivity of Al was varied 
from 0 to 100% in steps of 1% and ten different 
slopes were used for all but the exponential burn 
profile, which was modeled in 20 steps. The starting 
point was was not varied in the final runs, i.e. 
aluminium was assumed to ignite at the C-J plane 
starting from zero reactivity, which turned out to 
be true in all cases. Any initial reactivity that was 
initially tried always resulted in an overestimated 
plate position at the high pressure region and 
underestimated the plate position at low pressures. 
All this yielded approximately 4000 simulations 
total necessitating the 2D-approximation. 


Table 2. 
Explosive D/ms', calculated D/ms ', measured aie m-* % of TMD 
NTO 7890 8220 1760 92 
NNHHT 6470 6570 1300 76 
RDX 8630 8750 1800 99 
NTO/AI 76.5/23.5 7340 6780 1850 90 
NNHHT/AI 78.0/22.0 6470 6380 1540 83 
RDX/AI 73.8/26.2 7930 7740 1890 78 
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Experimental 
Materials 


Ecka AS 71 from Eckart Werke, Fiirth, Germany, 
nominally spherical aluminium powder was 
selected as the metal to be studied. This material 
was the finest powder available with a high 
aluminium content (>99%) with no carbonaceous 
impurities on its surface and had particles as 
spherical as possible. The specific surface area 
by BET (nitrogen—helium gas adsorption) was 
0.87 m e which indicates that the powder 
has considerable surface roughness, but still had 
Figure 11. Tool set for pressing the explosives excellent rheological properties and high apparent 
into rectangular bars. density. 


Two of the explosives (NTO and NNHHT) were 
synthesized and novel routes were developed for 
both to make kilo scale synthesis in a laboratory 
possible. Pure RDX was obtained commercially 
from Bofors, Karlskoga, Sweden. The syntheses 


Side view Top view 
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Wooden cap holder’ 
ee Hexogen booster 
Ey Explosive 
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; a Figure 13. Test charge installed and ready to fire. 
Figure 12. A schematic view of the test charge. 
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do not belong into the scope of this paper and will 
be published elsewhere. 


Charges 


The explosives were pressed into rectangular bars 
in a special tool constructed for the purpose shown 
in Figure 11. Copper plates were then attached on 
two opposite sides and an acrylic sheet was placed 
on the two other sides to keep the charge together 
and to increase confinement. An RDX booster (a 
20 x 20 x 20 mm cube) was placed on the top and 
a blasting cap holder was made out of plywood 
and installed on the top of the charge (Figures 12 
and 13). 


Calculations and experimental setup 


The explosive charges were designed as a 
compromise between producing as good a 
measurement resolution as possible and the easiest 
possible hydrodynamic simulations (Figures 10 
and 13). These requirements are complementary 
and a compromise had to be made. 


In theory, a cylindrical charge would be ideal for 
the modelling. Only a 2D-model is required to 
describe the problem accurately and the Jones— 
Wilkins—Lee equation of state (JWL-EOS) can be 
used, since it has been calibrated for cylindrical 
charges. On the other hand, a cylindrical charge 
is very problematic for experiments and practical 
calculations. It sends fragments in all directions, 


optical fibres 


which necessitates shielding of the X-ray cassette 
with heavy armor greatly attenuating the signal 
and reducing contrast, which is bad even without 
extra shielding. A cylinder charge can track the 
gas expansion only up to about | : 10, since the 
material has to stretch and breaks up at about 
this expansion ratio. Material stretching has to 
be taken into account in the simulations making 
the meshing much more difficult due to high 
distortions and Lagrangian meshing can hardly 
be used, which may be problematic due to the 
extensive computer resources needed and due to 
convergence problems. The cylinder wall in an 
X-ray radiograph is always smudged and has very 
low contrast because of geometric effects even 
without a fragment shield. 


The only advantage is that a cylinder wall position 
can be photographed, since the explosive gases 
do not obscure it until it breaks up. However, 
photography is still problematic due to shock 
waves that are formed in the surrounding air 
distorting the image and necessitating modelling 
with ray-tracing to remove the effect of the shock 
wave front in air distorting the recorded geometry. 
Alternatively, the explosive charge has to be 
placed in a vacuum to prevent gas shock waves 
from forming, which presents further engineering 
challenges. 


To avoid the problems above, a rectangular 


oe 


— 
—_ 


™ . 
witness plate 


~ 


_— 
be 
_ _ 
_ film cassette 
= 
— E] 
_. : 
_ lasting cap 1 
_ explosive charge 
: 500 mm 
: _ 1 
— 
~~ 
ae 


~ 


~. shields (plywood) 


~~ 
~~ 
~™ 
~ 
™~ 


~ 


Figure 14. Flash X-ray radiography setup (not to scale). The X-ray tubes and equipment were located 
inside of a bunker behind a 20 mm polycarbonate shield. Acceleration voltage: 180 kV, capacitor energy: 


400 J per tube, pulse width: less than 2 ns. 
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Figure 15. Flash X-ray radiography setup. Tubes 
are situated on the left inside the bunker, the 
charge stands in the middle and the man left 
hand points the film cassette. 


geometry was chosen (Figures 12 and 13). This 
geometry necessitates a 3D-model in theory, but 
it was concluded that 2D-model approximated the 
case to a precision of 2% at worst and in the high 
pressure region, no difference was seen. This small 
deviation is absorbed in the experimental errors, 
such as not having the charge exactly in parallel 
with the X-ray cassette and no 3D-modeling is 
needed. 


Inage filenane : findet3.$1 © Shot Number : 3 Camera: @ Tine | 00:00:00.0000000 Date : 09/06/1994 


Figure 16. The detonation of NTO/AI 76.5/23.5 
photographed with 4 frames overlaid. No sign 
of Al ignition is seen. Note how well geometry 
is preserved within the gas cloud from the cubic 
RDX booster. 
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The JWL equation of state could not be used for 
rectangular geometry, however, this was not a 
problem, since the JWL-EOS does not include 
explicit chemistry. As tracking chemical reactions 
was the goal of this study, JWL-EOS was out of 
question in any case. In a rectangular charge, the 
walls hardly have to stretch at all and distortions 
are so small that a Lagrangian mesh can easily be 
used saving huge amounts of computing time and 
increasing the time limit, where the hydrodynamic 
model fails. 


The measurement precision of a rectangular charge 
with X-ray radiography is much higher, since the 
contrast is sharp due to much less curvature of 
the walls compared to a cylinder. In theory, the 
precision is also higher with photographic methods 
for the same reason. Regardless of these measures, 
advanced digital edge detection algorithms and 
image processing were necessary to find the wall 
positions accurately from the radiographs. 


A rectangular charge sends metal fragments only in 
two directions and not towards the X-ray cassette, 
since that side of the charge is layered with 
plexiglass (Figures 12—15). Plexiglass fragments 
can be stopped with X-ray transparent materials, 
such as plywood and no reduction of contrast is 
caused. 


The plate positions were recorded using flash 
X-ray (Figures 14 and 15) equipment by Hewlett- 
Packard using 180 kV acceleration voltage and 
400 J capacitor energy for each pulse. X-Ray 
radiography was essential, since the entire 
wall profile and the location of the detonation 
front have to be known at a known moment of 
time. Hence, short circuit gauges etc. and point 
measurement methods are not applicable. Digital 
fast photography with a Hadland Photonics SVR 
CCD-camera was also attempted, but it failed to 
track the wall positions due to too low contrast 
(plates not visible) or obscuring by explosive gases 
or both (Figures 15—17). Detonation velocity was 
recorded with an optoelectronic counter (Explomet 
by Kontinitro) and optical fibers were attached at 
two locations on the explosive surface as shown 
in the schematics (Figure 12). Detonation velocity 
was also obtained from the photographs, but with 
less accuracy. 


A high power flash bulb (Hadland Photonics 
experimental prototype) and a white screen were 
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> 
4 Camera: @ Tine : 00:00:00.0000008 Date : 69/06/1994 


Figure 17. The detonation of NNHHT/AI 78/22 
photographed with 4 frames overlaid. No signs 
of Al ignition are seen here either. The gases are 
more luminous due to higher carbon content and 
higher temperature than with NTO/AI. 


installed behind the charge to provide a luminous 
background against which the plates were expected 
to be seen as shadows. Unfortunately, this did not 
work either, no plates were seen in the photographs 
(Figures 16—18). This was particularly problematic 


= 


Inage filename | findetS.SAL Shot Nunber : 5 Canera: @ Tine : 00:00:00.0000000 Date : 05/06/1994 


Figure 18. The detonation of RDX/AI 73.8/26.2 
photographed with 4 frames overlaid. The gas 
cloud is highly luminous and opaque showing the 
aluminum ignites and burns. 
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with RDX/AI, which ignited the Al and generated a 
highly luminous and opaque gas cloud (Figure 18). 
Even though the shape of the cloud does reflect 
the position of the walls, the position cannot be 
determined with high enough precision. 


Flash X-ray setup was done similarly (Figures 14 
and 15). Two X-ray tubes triggered in a sequence 
were used and the charge and the latter tube were 
aligned to record an image at right angles to the 
film plane with respect to the longitudial and 
perpendicular axes of the charge, centered at the 
charge middle point. The first X-ray tube was then 
off-axis in the longitudial direction only below the 
first tube and the angular error thus created was 
known. The angular correction was made at the 
visualization step, i.e. the calculated images were 
distorted by the same amount as the experimental 
images, because digitally generated visualization 
data were easier to manipulate than the experimental 
radiographs. The tubes were triggered in the order 
that produced the more distorted image first and 
then the less distorted, center aligned tube was 
fired. This was done because the latter image is 
more important revealing the situation after larger 
gas expansion. First images served as a backup, 
should the second X-ray for instance flash miss 
the detonation front and were actually not needed 
at all, since all rounds were successful. 


The distance between the X-ray tubes and the 
charge was 1500 mm and that of the film plane 
to charge was 500 mm. The film cassette was 
embedded inside a shield of 5 to 8 stacked, 10 mm 
thick sheets of plywood. The amount of plywood 
was kept to a minimum and estimated based on the 
charge energy, more for RDX and RDX/AI, less 
for NTO/Al. RDX/AI still punctured the shields, 
which can be seen as an artifact on the radiograph 
in Figure 24. 


The X-ray radiographs were taken on regular 
X-ray film and the images were digitized at the 
Center for Biotechnology at the University of 
Turku with an image analysis system designed for 
imaging protein sequencing gels with low contrast 
and luminosity. The equipment was modified by 
exchanging the UV-lighted table for a visible light 
table. A digital camera connected to a computer 
was used to image the radiographs and 16 to 20 
overlays were used to increase contrast and filter 
out noise. A Prewitt filter was used as the edge 
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detection and run on the images after digitization 
for further improvement of the precision. The 
digital camera resolution was 640 x 480 pixels 
and limited the maximum resolution of the images 


presented in this study. 


Results and discussion 


Two charges each of pure NTO, NNHHT, RDX 
and mixtures NTO/Al, NNHHT/AI and NTO/AI 
were detonated and recorded with the fast camera 
and with flash X-ray. The images and the results 
from the simulations are shown in Figures 16—24. 
Detonation velocity data are presented in Table 2. 


Aluminium with NTO 


For NTO/AI the calculation yields a 1850 K 
C-J temperature, if all Al is inert. Such a low 
temperature will not promote an ideal detonation, 
so the observed detonation velocity is lower than 
the calculated one as calculations do not account 
for all the possible losses. The Al has to be inert, 
since the measured detonation velocity is still 
considerably higher than the calculated one, if 
100% Al is assumed reactive in the calculation. 
Note that in the case of NTO/AI the detonation 


Figure 19. The detonation of pure NTO flash 
X-rayed. From the left: calculated flow pattern, 
measured radiograph (two frames overlaid) 

and the two previous images overlaid. The fit is 
perfect showing the 2D-calculation is sufficient to 
describe the case. 
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Figure 20. The detonation of pure NNHHT flash 
X-rayed. From the left: calculated flow pattern, 
measured radiograph (two frames overlaid) and 
the two previous images overlaid. Again, the fit 
is perfect. The crack in the left plate is due to a 
small void in the explosive charge on that side. 
NNHHT was extremely difficult to formulate to a 
crack-free, homogeneous charge and the above 
radiograph is the best available. 


Figure 21. The detonation of pure RDX flash 
X-rayed. From the left: calculated flow pattern, 
measured radiograph (two frames overlaid) 
and the two previous images overlaid. Even in 
this case the fit is perfect. The extra artifacts 
originate from a loop in the ignition cable, 

that had fallen between the charge and the film 
cassette. 
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Figure 22. The detonation of NTO/AI 76.5/23.5 
flash X-rayed. From the left: calculated flow 
pattern, measured radiograph (two frames 
overlaid) and the two previous images overlaid. 
The fit is perfect when no Al is allowed to react at 
any stage. 


Figure 23. The detonation of NNHHT/AI 78/22 
flash X-rayed. From the left: calculated flow 
pattern, measured radiograph (two frames 
overlaid) and the two previous images overlaid. 
Again, the fit is perfect when no aluminum is 
allowed to react. 


velocity vs. reactivity has a negative slope and 
a lower velocity indicates a higher degree of Al 
reaction. As it is unlikely that the low detonation 
temperature will ignite Al the observed, lower 
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Figure 24. The detonation of RDX/Al 73.5/26.5 
flash X-rayed. From the left: calculated flow 
pattern, measured radiograph (two frames 
overlaid) and the two previous images overlaid. 
In this case, a perfect fit can be made when 17% 
of the aluminum is allowed to react according 
to the exponential burn model with a slope +1 


(Figures 6 and 7) corresponding to an exponent 
+1/3. 


than calculated 0% reactive velocity is due to 
unaccounted losses and not due to partial reactions 
of aluminium. 


Loose Al powder particles were found in the 
surroundings after the NTO/AI charge had been 
fired. Thus, there is strong experimental evidence 
no Al reacted with NTO. Moreover, by observing 
the photograph of the detonation in Figure 16, it 
is clear no Al ignited. Otherwise, the detonation 
products would not have been transparent. 


The hydrodynamic model of NTO/AI matches 
the measurements perfectly when all of the 
Al is assumed inert throughout the process. If 
100% reactivity is assumed, the model greatly 
overestimates the copper plate energy. If any 
of the burn models are used, the energy is also 
overestimated or the copper plate profile is not 
reproduced (Figure 25). 


Journal of Pyrotechnics, Issue 26, Winter 2007 


Figure 25. The detonation of NTO/AI 76.5/23.5 
with 32% Al reactivity using the logarithmic 
burn with a slope +3. Note how much the plate 
position is overestimated. 


Aluminium with NNHHT 


With NNHHT/AI, the theoretical detonation 
temperature is 2300 K at the density of the 
experimental charge with non-reactive Al. As 
a result, the measured NNHHT/AI detonation 
velocity is closer to the calculated one, although 
still somewhat lower. The detonation velocity 
indicates some of the Al reacts, but a conclusion 
cannot be made based on this information only. 
It is more likely that losses reduce the detonation 
velocity to lower than calculated values as with 
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Figure 26. The detonation of NNHHT/AI 78/22 
with 48% Al reactivity using linear burn with 

a slope 0. Again, the plate position is greatly 
overestimated. 


NTO/AI, but the difference is less due to higher 
temperature and faster reactions promoting less 
unideal detonation characteristics. As the slope 
of the A-D dependency is negative in this case as 
well, a match is found at some point between A = 0 
andj= 1. 


Unburnt aluminium was found in the surroundings 
after the NNHHT/AI explosion as well. The nearby 
X-ray film cassette was coated with a film of Al, as 
if silver paint had been sprayed onto the surface. 
The Al probably melted in the detonation, but still 
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Figure 27. The detonation of RDX/AI 73.5/26.5 
with 96% Al reactivity using exponential burn 
with a slope of —10. If this was possible to 
realize, the plate push ability of the explosive 
would be greatly enhanced. 


did not ignite. The higher luminosity of the higher 
energy and higher carbon content booster charge 
indicates the glow in the gas cloud originates from 
glowing carbon particles (Figure 17). 


The hydrodynamic simulation yields 0% 
reactivity in this case as well, i.e. inert Al. Any 
other reactivity or any other burn model yields 
unmatched experimental and calculated results 
(Figure 26). 
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Aluminium with RDX 


RDX was the only explosive that was able to 
ignite the aluminium. The calculated detonation 
temperature of RDX/AIl was 2500 K with inert 
aluminium and heat losses taken into account, 
which is high enough to ignite the aluminium. 
Therefore RDX and RDX/AI, are nearly ideal 
explosives with fast reactions and a thin reaction 
zone and their detonation properties can be readily 
modelled. 


The detonation velocity of RDX/Al is somewhat 
lower than calculated, but the difference is smaller 
than that of either NTO or NNHHT and falls within 
the experimental error. Again, losses due to inert 
Al are the likely cause, although the detonation 
velocity also fits to about 20% reactivity of Al. 
Surprisingly, that happens to be almost exactly the 
reactivity that is detected from the gas expansion 
and hydrodynamic model, so in this case also 
the detonation velocity matches the observed 
reactivity, although the aluminium reacts far away 
from the C-J plane and explosive reaction zone. 
Despite a good match, this information should not 
be trusted. The reactivity based on the detonation 
velocity is highly dependent on which covolume 
set was used in the calculations and any reactivity 
can still be found by using different covolume sets 
for the product gases in the calculations. 


With RDX/AI, the Al ignites and even contributes 
to the mechanical energy to some degree. The gas 
cloud is opaque and highly luminous proving the 
aluminium is burning (Figure 18). In the case of 
RDX, there is more carbon present in the gases 
than with other studied explosives even without 
the Al reacting. However, the pure RDX gases 
are not totally opaque, as can be seen from the 
other photographs (Figures 16-18). In all of them, 
the RDX booster causes more luminosity than 
either NTO or NNHHT, but the gases are still 
partially transparent. Hence, the opacity cannot be 
explained by carbon only. 


Hydrodynamic modelling of the RDX/AI reveals 
much more information. If Al is let burn while 
the gases expand, the resulting hydrodynamic 
calculation can match the plate position perfectly 
(Figure 24). More reactivity yields non-matching 
results (Figure 27). The best fit is obtained when 
Al is reacting according to the exponential reaction 
scheme with a slope +1/3 (Figure 24), starting 
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Exponential v-’ —space for hexogen/Al, slopes —10 to 10 


In (vi(cm*/g)) 


Figure 28. Exponential burn model v-A space for RDX/AI 73.5/26.5, end point at (3A=0.17). The slope 
+1 corresponding to an exponent +1/3 is depicted as a thick line and represents the reaction path that 


reproduces the RDX/Al charge geometry. 


from zero reactivity at the C-J plane and ending at 
17% reactivity at In v = 3. Any other combination 
of reaction paths yields non-fitting results with 
+1% tolerance in | and +1/3 tolerance in reaction 
path slope. The v-’ space and the corresponding 
v—-T space are depicted in Figures 28 and 29. 


The Al burning during gas expansion at lower 
pressures contributes to the mechanical energy, but 
not to the shock energy of the explosive. Even then, 
only 17% of the Al is taking part in the reactions 
during the time scale of interest, i.e. during a short 
enough time to affect the mechanical energy of 
the expanding gases. The rest of the Al burns later 
creating an airblast after the detonation, but does 
not show within the scope of the hydrodynamic 
model. No unburnt aluminium was found in the 
surroundings after the detonation of the RDX/Al 
charges. 
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Ignition and combustion of 
aluminium 


Detonation temperature was found to be the 
limiting factor for aluminium ignition. There is 
even previous indication that this might be the 
case.*? If the temperature is not high enough 
(2500 K), Al does not ignite at all. If the ignition 
temperature is reached, as in the case of RDX/AI, 
the Al still does not react at the C-J plane but burns 
afterwards with an orders of magnitude slower rate 
than the explosive itself as also observed in other 
studies.'”*! The burn rate of aluminium has been 
found to be diffusion limited through the oxide 
shell formed at the surface.' At high pressures, the 
oxide shell is always formed, since the pressure 
prevents the metal or the oxide from evaporating 
and a similar, diffusion limited rate can be assumed 
as at lower pressures. 


At low pressures, aluminium has been observed 
to have a pressure exponent of unity,’ although 
the study was carried out on activated aluminium 
powders, where the diffusion rate was accelerated 
and a higher exponent was expected. 
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Exponential v-T —space for hexogen/Al, slopes —10 to 10 


In (v/(cm?/g)) 


Figure 29. Exponential burn model paths for RDX/AI 73.5/26.5 in the v-T space with maximum 4 = 0.17. 
The path corresponding to an exponent + 1/3 is emphasized. The x-marks are BKW calculated points, 
dash—dot lines are partial immediate reaction isentropes and the solid thin lines are fourth degree least 


squares polynomial fits to the calculated points 


In the explosive mixture a low positive value of 
+0.33 was found. This is likely to be due to the 
fact that at high pressures aluminium oxide and 
aluminium metal no longer evaporate generating 
at least a molten aluminium oxide shell on the 
particle. While such a shell does not prevent 
combustion, it does slow it down, since the 
oxidizer has to diffuse through a liquid barrier 
before reaching the aluminium metal below it. 
This effect should reduce the pressure exponent, 
since increasing the oxidizer concentration (e.g. 
pressure) outside of the particle surface has less 
effect on the particle oxidizing rate due to a liquid 
barrier between the fuel (Al) and the oxidizer 
(product gases from the explosive). At low 
pressures this barrier is removed, since aluminium 
oxides and the metal itself can evaporate and get 
mixed with the oxidizing gases and the burn rate 
becomes directly proportional to pressure, e.g. the 
pressure exponent approaches unity. 
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Conclusions 


Aluminium ignition and combustion 
characteristics in an explosive mixture were 
successfully determined and measured. The 
procedure is difficult and tedious and not suited 
for routine work, but it was nevertheless shown 
to be possible to carry out. The procedure can 
be improved by introducing an automatic global 
optimization algorithm fitting the calculations and 
the experiments and providing automatic feedback 
to the hydrodynamic model. 


Aluminium was shown not to contribute to the 
shock energy of the explosive; on the contrary, it 
reduced the shock energy from the theoretical zero 
reactivity value more, the lower the detonation 
temperature of the explosive was and caused 
unaccounted losses and non-ideality to the 
detonation. 


Aluminium ignition temperature even under 
detonation was shown to be the decomposition 
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temperature of aluminium oxide, e.g. 
approximately 2500 K. If the explosive does 
not reach this temperature with inert aluminium 
including heat losses to the Al, the aluminium 
powder will not ignite at all. This finding was 
somewhat surprising considering the fact that 
the aluminium powder passes a detonation shock 
front with extreme pressures, pressure gradients 
and shear forces, which should break up the oxide 
shell, heat up the particles considerably due to 
adiabatic compression and thus greatly enhance 
the probability of ignition. 


Aluminium was found to burn in a propellant-like 
manner during the gas expansion and its burning 
rate can be best described with the traditional burn 
law of the type r = Ap". The pressure exponent 
was found to be +0.33, i.e. slightly positive. No 
significant reactions of aluminium were detected 
at the C-J plane or the pure explosive reaction zone 
thereafter. If, however, the aluminium combustion 
is considered a reaction zone, the zone extends 
from the detonation front down to ambient pressure 
and can span very large distances. 


A maximum reactivity of 17% in RDX/AI was 
observed while aluminium did not ignite at all 
in NTO/AI and NNHHT/AI. In order to reach 
higher degrees of reactivity within the pressure 
range, where most of the mechanical energy of the 
explosive is released, considerably finer powder 
needs to be used than was applied in this study. 
The efficiency of Al combustion even in the ideal 
case is still low: only about 30% of the added 
energy can be converted to mechanical energy 
while the base explosive has a considerably 
higher thermal efficiency. More mechanical 
energy can be produced if working fluid from the 
surroundings can be introduced, such as water 
mist in an underwater explosion, but such energy 
is only available at low pressures through a blast 
wave and does not contribute to the primary energy 
release mechanisms of the explosive. Aluminium 
energy can be best utilized with nitrogen rich 
explosives, which provide large amounts of 
relatively inert working fluid to convert thermal to 
mechanical energy. Nitrogen rich explosives are 
also necessary to determine the reaction profile of 
Al under detonation conditions unambiguously. 
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Ignition And Burning Rate Characteristics Of Pyrolant 
Black Powder 
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Abstract: Composite propellants contain binder as fuel which connects oxidizer particles and metals. The 
burning rate of a propellant is affected by the concentration and type of binder. Pyrolant is mixed with 
oxidizer particles and metal particles. Pyrolant does not contain binder, so burns smoothly. The burning 
rate of black powder (BP) which is a kind of pyrolant is higher than that of AP composite propellant. The 
shorter the reaction time in the gas phase near the burning surface of black powder becomes, the higher the 
burning rate becomes. The burning rate of BP is inversely proportional to the ignition delay time. 


Keywords: Black powder, pyrolant, ignition delay time, burning rate 


Introduction 


Pyrolant consists of mixed oxidizer particles and 
metal particles. Oxidizers such as ammonium 
dinitramide (ADN), ammonium perchlorate (AP), 
ammonium nitrate (AN), hydrazinium nitroformate 
(HNP), potassium nitrate (KN), cyclotrimethylene 
trinitramine (RDX), etc. are used. Metals such 
as aluminum (Al), boron (B), magnesium (Mg), 
titanium (T1), zirconium (Zr), etc. are used. Several 
kinds of pyrolants are used for pyrotechnics, as gas 
generator, ignition material, heat supply, etc. 


These particles are electrified and easily ignited 
by static electricity. The sensitivity of a pyrolant 
is generally high. The ignition delay time of a 
pyrolant is shorter than that of an AP composite 
propellant. The burning rate of a pyrolant is 
easily controlled by changing the concentration of 
oxidizer.'> 


There are two types of pyrolants. One type 
burns homogeneously, oxidizer and fuel melt 
at the burning surface and it is easy to evaluate 
the burning mechanism. The burning rate of 
black powder (BP) is not changed by the sizes of 
potassium nitrate (KNO;), charcoal (C), and sulfur 
(S). Black powder belongs to this homogeneous 
type. The other type burns heterogeneously, not 
melting at the burning surface and the burning 
rate changes when the size of the compositions 
changes, as Mg/TF, B/KNO and it is difficult to 
evaluate the burning mechanism. 
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We used black powder in this study. Standard 
compositions were KNO;/C/S = 75/15/10 wt%. 
Charcoal was made from pine. By changing the 
concentration of these materials, the burning rate 
was controlled. We studied the combustion wave 
structure near the burning surface with small 
thermocouples, and the relationship between 
burning rate and ignition delay time. 


Theoretical investigation 
(1) Heat balance at burning surface 


The burning rate, 7, is defined by the following 
equation:*” 


rp tC, (TZ, -h)-Q.3 


LT. 
SA a 
(1) 
fe) =U, (T +Ty) 
rp 
=—(t,,+T.) 
g 
ro.RT. 
= Pe * (7, +7) 
‘ (2). 


In equation (2), we use the mass balance equation, 
and equation of state, 
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rp, = ePg 
P 
Pe Rep 
g 
10, 104K 
g p; P 


(3) 


where p, = density of pyrolant, C, = specific heat 
capacity, 7, = surface temperature, 7) = initial 
temperature, O, = heat of reaction at the surface, 
A, = heat conductivity in the gas phase, ®,, = 
temperature gradient near the burning surface in 
the gas phase, 7;= flame temperature, 6 = reaction 
length, u, = gas velocity, t,, = physical process 
time, as mixing and diffusion, 74, = chemical 
reaction time, p, = gas density, R = gas constant, 
and P,=pressure. The adiabatic flame temperature 
T; is calculated with NASA-CEA.® In this study 
P,= 0.1 MPa, p,, C,, and R are constant. If 7, and 
Q, are constant, the burning rate will be obtained 
following equation (4): 


jf’ 
ae ae © 
T; (Typ +T 4) Top ey 


(4) 
where a@ is a constant. If 7, < 7; and 1, = 0, 
equation (4) will become 


oo 
cy (5) 


The burning rate is determined by the chemical 
reaction time. This chemical reaction time is 
obtained from equation (6): 


E 
Ty =Zy, exp(—n) 
sb (6) 


where Z,, is a constant and Ey, is the activation 
energy of the pyrolant. 7,, is obtained following 
equation (7): 


T, +T. 
= a. 

(7) 
where 7; is the flame temperature and 7, is the 
temperature of the burning surface. 

(2) Ignition characteristics 


The surface temperature of the pyrolant increases 
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with passing time when pyrolant is put into the 
furnace, and reaches decomposition temperature, 
and gas is generated at the surface of the pyrolant. 
The generated gas reacts near the surface and heat 
is conducted from the reaction gas to the surface. 
A luminous flame appears and the pyrolant 
ignites. Ignition delay time tj, can be separated 
into two terms, the first of which is physical delay 
time t,, which is the time to reach the pyrolant 
decomposition temperature, and the second is the 
chemical delay time z,, which is the time from 
decomposition to appearance of a luminous flame. 
The ignition delay time is obtained following 
equation (8): 
Tig =T, +7, (8) 
Tig 1s measured experimentally, and the physical 
delay time t, is defined as the minimum ignition 
delay time. The chemical delay time is obtained 
following equation (9): 
a aa (9) 
The chemical ignition delay time t, is the same as 
the chemical reaction time t,, of equation (6). 


E 
TR 

‘ (10) 
The decomposition gas reacts in the gas phase 
near the burning surface when it burns or ignites, 
therefore Z, = Z,,, E,= Eg, and T,= Ty,. If physical 
ignition delay time t, < 1, we will obtain the 
following equation: 


7, = Z, exp( 


1 
2log(—) « log(z,,) = log(z,) = log(z,,) 
7 (11) 


Burning rate is inversely proportional to ignition 
delay time. 


Experimental 
(1) Thermal decomposition 


The thermal decomposition characteristics were 
obtained with thermo-gravimetric (TG) and 
differential thermal analysis (DTA). The heat 
of reaction at the surface (Q,) was measured 
using differential scanning calorimetry (DSC). 
These experiments were conducted in N, gas 
atmosphere at 0.1 MPa and the N, gas flow rate 
was 50 ml min '. The heating rate was 10 K min ! 
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and the sample mass was 3 mg for TG, DTA and 
DSC experiments. The mean diameter of KNO, 
particles was 30 um, that of S was 30 pm, and that 
of C was 10 um. 


(2) Burning rate 


To investigate the effect of sulfur on burning 
rates, the burning rate was measured with varying 
concentrations of sulfur in an atmosphere at 
0.1 MPa, with a KNO,/C ratio of 75/15. The 
concentration of oxidizer was changed and the 
effect of oxidizer on burning rate was measured, 
with C/S= 15/10. The pellet was 10 mm in diameter 
and 10 mm in length and cigarette burning type 
was used. The burning rate was measured with a 
video camera. The measurement was repeated 3 
times for each set of conditions. 


(3) Ignition delay time measurement 


The ignition delay time was measured changing 
the concentration of sulfur and temperature in the 
furnace at atmosphere. The furnace inner diameter 
was 36 mm, and depth was 36 mm. The furnace 
bottom temperature was changed from 473 K to 
1000 K. The sample weight was about 10 mg and 
samples were dropped from the upper side of the 
furnace. After dropping the sample, a blue flame 
appeared first from the combustion of sulfur, and 
a luminous flame followed. The time between the 
drop and appearance of the luminous flame was 
measured with a high speed video camera. The 
measurement was repeated 30 times for each set 
of conditions, and the time for 50% probability of 
ignition was defined as the ignition delay time. 


(4) Temperature measurement near burning 
surface 


The burning rate of BP was very high, and it was 
difficult to measure the temperature gradient near 
the burning surface. To measure the temperature 
gradient, the pressure was set below atmospheric 
pessure, however at pressures lower than 20 kPa 
BP did not burn. We measured the temperature 
gradient to evaluate the effect of sulfur 
concentration at pressures between 30 kPa and 
70 kPa. In this pressure range, burning rate was 
still high. The pellet was 10 mm in diameter and 
10 mm in length and cigarette burning type was 
used. The line diameter of the thermocouple was 
12.5 um and it was Pt/Pt-10%Rh. The samples 
were without sulfur (KNO;/C = 75/15), and with 
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sulfur (KNO,/C/S = 75/15/10 wt%). 


Results and discussion 
(1) Thermal decomposition characteristics 


Figure 1.1 shows the results of DTA and TG 
without S (KNO,/C = 75/15). Without sulfur a 
large decrease started at 670 K and ended at 760 K 
and indicated the reaction of black powder. The 
endothermic peak of KNO; appeared at 607 K and 
was the melting point. The exothermic reaction 
peak appeared 728 K. 


Figure 1.2 shows the results of DTA and TG at €, 
= 0.25 (KNO,/C/S = 62.5/12.5/25 wt%). Sulfur 
melted at 433 K and evaporated at 533 K. The 
endothermic peak of KNO; appeared at 582 K and 
was the melting point. The exothermic reaction 
peak appeared at 705 K at €, = 0.25. KNO, and S 
melted at the burning surface. The decomposition 
temperature decreased with increasing ¢,. 


With sulfur two large decreases appeared: the first 
decrease started at 433 K and ended at 533 K. The 
decrease caused by the sulfur evaporation was 
20.8% at € = 0.25 and not all the sulfur content 
evaporated. Sulfur was considered to exist in the 
compound generated by the reaction between sulfur 
and KNO; or charcoal. The heat of reaction at the 
surface was about 2100 kJ kg ' with DSC and also 


p \ 
<——>, 


Decomposition of BP 


300 550 800 
Temperature, K 


Figure 1.1. Results of DTA and TG without S. 
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Evaporate of sulfur 


Decomposition of BP 


300 550 800 
Temperature, K 


Figure 1.2. Results of DTA and TG with €, =0.25 
(25 wt%). 


did not change with the sulfur concentration. 
(2) Burning rate 


The relationship between burning rate and 
concentration of sulfur is shown in Figure 2. The 
calculated adiabatic flame temperature and C* are 


also shown. C* is the characteristic velocity and 
1 
is proportional to T',2. 
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Figure 2. Burning rate characteristics. 
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Figure 3. [gnition characteristics. 

The characteristic velocity was about 787 m s | 
which was lower than that of AP composite 
propellant. The adiabatic flame temperature was 
over 1500 K at €, < 0.20 and the flame temperature 
was higher than the decomposition temperature. 
Burning rate increased with increasing ¢, and 
reached a maximum burning rate of 151 mm s ' 
at €,= 0.25. Burning rates were not proportional 
to the adiabatic flame temperature. Burning rates 
increased from 2.2 mm s! to 10.5 mm s ', when 
the concentration of oxidizer was changed from 
85 wt% to 60 wt% at a ratio of C/S = 15/10. In 
this case the concentration of sulfur was increased 
from 6 wt% to 16 wt%. The sulfur melted and 
evaporated and reacted in the gas phase, therefore 
these results showed that the heat conducted 
from the gas phase to the burning surface was an 
important factor for burning rate. 


(3) Ignition characteristics 


The relationship between furnace temperature and 
ignition delay time is shown in Figure 3 at ¢, = 0, 
0.1, 0.25. The concentration of S was 0 (é = 0), 
and below 598 K the sample did not ignite. Ignition 
delay time decreased with increasing furnace 
temperature and also decreased with increase of 
concentration of S. These results showed that 
sulfur affected the ignition delay time. 


The ratio of chemical ignition delay time, 7,+, is 
shown in Figure 4 as parameter of €,. The ratio 
of chemical delay time was smaller than 0.5 at 
furnace temperature >573 K. The ignition delay 
time of BP was occupied by physical processes. 
The ratio of chemical ignition delay time was 
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Figure 4. Ratio of chemical delay time. 


constant at changing sulfur concentration, when 
the temperature was between 573 K and 773 K. 
Equation (11) was modified as follows: 


T, = AT, 
logz, =loga +logr,, 


a = constant 


1 
logr, « logz,, « 2log— 

. (12) 
The ratios of zt, fluctuated when the furnace 
temperatures were 523 K and 823 K, and so 
equation (12) could not be used. 


(4) Heat balance at burning surface 


The temperature profiles near the burning surface 
were measured with a small thermocouple and are 
shown in Figure 5. These temperature histories 
were without S and with € = 0.1 (10 wt%) at 
pressure = 30 kPa. Solid phase is shown at distance 
<0, and zero means at the burning surface. The 
temperature increased smoothly in the solid phase, 
and a large variation occurred in the gas phase. 
The temperature gradients near the burning surface 
were obtained from these temperature histories. 
The surface temperature (7,) and temperature 
gradient (®,,.) at 0.1 MPa were estimated from 


Table 1. Surface temperature (T,), heat of 
reaction (Q,), and temperature gradient near 
burning surface (®,,) 


Parameter é,=0 é€,=0.1 
TJ/K 783 692 
OJkJ K"! 2160 2180 
®.x107/Kmm! 5.0 14.8 
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Figure 5. Temperature histories near burning 
surface. 


these low pressure values when pressure changed 
from 30 kPa to 70 kPa. The heat of reaction (Q,) 
was obtained with DSC. Surface temperature 7,, 
Q,, and ®,, are shown in Table | at a pressure of 
0.1 MPa. The surface temperature was 783 K 
without S, and was 692 K with S and decreased 
roughly 100 K and the variation was smaller than 
that of ®,,. The heat of reaction (Q,) at the surface 
was about 2100 kJ kg ' and also did not change 
with sulfur concentration. The temperature 
gradient near the burning surface in the gas phase 
increased with increasing concentration of sulfur. 
The burning rate was controlled by the reaction 
rate in the gas phase. 


(5) Relationship between burning rate and 
ignition delay time 


The relationship between burning rate and ignition 
delay time is shown in Figure 6 at 0.1 MPa. The 
concentration of sulfur changed and the furnace 
temperature changed from 573 to 823 K. The 
burning rates were obtained by changing the 
concentration of sulfur. Ignition delay time and 
burning rate showed a strong relationship. The 
burning rate was inversely proportional to the 
ignition delay time. 


The surface temperature 7, and QO, of BP did not 
change when the concentration of S changed from 
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Figure 6. Relationship between burning rate and 
ignition delay time. 


0 tol0 wt%. It was possible to use equation (4). 
KNO, and S melted at the burning surface, so 
the physical process time t,, = 0. The ratio of 
chemical ignition delay time was constant, when 
¢, was changed from 0 to 20 wt%, so 1, ~ Arig was 
obtained, and it was possible to use equation (12). 


Conclusions 


The following statements of BP were obtained in 
this study. 


(1) Burning rate increases with increasing sulfur 
concentration at €, < 0.25. 


(2) The ignition delay time decreases with 
increasing furnace temperature and also with 
increasing S content. 


(3) The burning rate is controlled by conductive 
heat from the gas phase to the burning 
surface. 


(4) The burning rate is inversely proportional to 
the ignition delay time. 
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Abstract: An experimental study of laser ignition of Zr/KCIO,, Zr/PbCrO,, Fe/KCIO, Fe/KNO;, B/KCIO, 
and B/KNO; mixtures is presented in this paper. The ignition system is composed of two lenses which focus 
the laser beam on the tablet. A sapphire porthole, juxtaposed with the pyrotechnic composition, protects the 
optical unit but causes losses of heat which result in an increase of the energy threshold. 


This experimental set-up was used to determine the ignition sensitivity of five mixtures. A particular 
behaviour of B/KNO; was observed. Experimental results pointed out an increase of the energy of ignition 
when the laser power grows. The other mixtures exhibit a conventional behaviour of the trend of ignition 
threshold E59 according to the power P of the laser beam: a continuous decrease of the value of Eso is 
observed when P is rising. 


Some parameters linked to the experimental device also have a great influence. Among them the power 
density seems to be predominant but the thermal conductivity of the sapphire windows also plays an 
important role. It is possible to find an optimum value of the laser diameter, for a given power density. 
This experimental and parametric study have shown that, on such devices, it is possible to find optimal 
conditions of the ignition of various pyrotechnic mixtures. 


A modelling, based on a progressive absorption of the laser beam inside a reactive pyrotechnic composition, 
was developed. It takes account of the energy exchanges between the ignition system and the pellet. Two 
interface parameters make it possible to optimize calculations. One relates to the heat transfer on the 
level of the porthole/tablet interface, the other relates to the laser absorption of energetic materials. By 
combining these two terms, it has been possible to corroborate the experimental thresholds on five different 
pyrotechnic mixtures. 


To evaluate the interface parameters, whose implication in initiation is critical, several methods are 
exposed. One of them proposes a numerical calculation based on a random and automatic installation of 
the grains. The result gives Gaussian variables. This method enables simulation of the statistical tests of 
Bruceton and Langlie which are used in pyrotechnics. 


The numerical results show that the sensitivity of the mixtures depends primarily on the propagation of 
heat towards the porthole and the interior of the tablet, as well as on the in-depth absorption of the laser 
in the pellet. 


Keywords: Laser diode, ignition, pyrotechnic mixtures, numerical modelling 


advantages. The laser diodes are highly reliable. 
Because of their small size and light weight, they 
can be easily used in embarked ignition systems. 
The ignition of different pyrotechnic mixtures, 
using this method, was carried out under various 
experimental conditions in the LEES laboratory.* © 


Introduction 


Laser ignition of energetic materials is of great 
interest for practical or fundamental studies and 
applications. Since the emergence of new laser 
devices several studies have appeared.'” 


Among the various methods of lighting, ignition 
by means of a laser diode presents many 
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Among them, boron mixtures and Fe/KCI1O, are 
presented in this paper. Ignition of boron mixtures 
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was investigated as soon they could be used in the 
field of propulsion and studied previously in this 
context.’ * 


Firstly the experimental device and test procedure 
are described. The influence of the laser beam 
diameter is demonstrated and a comparison 
between KCIO, and KNO, oxidizer is presented. 


In order to explain the influence of all parameters, 
a numerical modeling of the solid state ignition 
process will be described. This model needs 
the numerical values of equivalent thickness of 
air at the porthole/pellet interface and the depth 
of absorption of the laser beam at the surface 
of the solid. By means of an optimization, the 
determination of these parameters is possible 
and allows the numerical simulation of various 
pyrotechnic mixtures during ignition. 


The validation of this model is then achieved by 
means of a comparison with experimental data 
coming from studies carried out in our laboratory. 


Experimental 


This type of initiator uses a system of two lenses 
to focus the divergent beam on the surface of the 
pyrotechnic composition (Figure 1). A detailed 
description of the experimental set up is given in 
ref. 3, 4 and 9. 


The laser used in this work was a diode provided 
by Opto Power company and referenced OPC- 
A007-1047-FC. The maximum output power 
of this laser diode is 7 W and its wavelength is 


Sapphire Window 


A=815 nm 


Optical fibre ® = 1 mm NA=0.14 _L, 


PMMA cylinder 


820 nm. It delivers the luminous radius at the end 
of a sheathed optical fiber. A sapphire window, 
of 1 mm thickness, separates the optical system 
from the pellet to avoid damage to the system 
during combustion. The end of the optical fiber 
is positioned on the focal point of the first lens, 
whereas the pellet is located on the focal point of 
the second one in front of the sapphire window. 
The choice of the diameter of the laser spot is 
limited to four values (10, 16, 25 and 25 mm), 
according to the focal distances of the two lenses. 


Experimental results 


A complete experimental study was carried out 
and reported in ref. 4. In this work, the influence of 
the following parameters on the energy threshold 
of ignition was presented: 


* the mass percentage of boron 
¢ the diameter of the spot 

e the power density 

e the oxidizing agent 

e the porosity ¢ 

e the percentage of opacifier. 


Among these parameters, the choice of the diameter 
of the laser beam and the nature of oxidizing agent 
will be presented in the following section. 


Results of ignition for B/KNO, 


The result of ignition of B/KNO, pellet with a 
porosity of 21.5% is presented in Figure 2. The 
influence of mass percentage of boron is very 


Infrared Sensor and filter 


Pyrotechnic 
==" composition 


Brass closure 


Steel confinement 


Figure 1. Schematic diagram of the experimental set-up. 
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Figure 2. Measured energy thresholds for 
various amounts of boron (porosity € is 21.5%). 


iN 
oO 


important and the Ey energy values decrease from 
400 Jcm* to 5Jcm~ when the percentage of 
boron increases from 15 to 40%. For the following 
studies the choice of 20% of boron for all pellets 
is assumed. 


Another important experimental parameter is the 
laser beam diameter. In our case four diameters 
have been used and the result of the experiments 
is reported in Figure 3. These results indicated the 
existence ofa minimum energy fora constant power 
density equal to 2.2 kW cm”. Another result, not 
presented here, concerns the standard deviation. 
This parameter continuously decreases when the 
laser beam diameter grows. The compromise of 
180 um diameter offers a better performance for 
the use of energy at the pellet surface and this value 


Energy (mJ) 
5 


100 200 300 400 500 
Diameter (ym) 

Figure 3. Energy thresholds for various 

diameters of the laser beam (P/S = 

2.45 kW cm). 
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Figure 4. /nfluence of power density on the ignition 
thresholds for B/KNO; and B/KCIO, mixtures (30 wt% 
of boron, porosity € is 20.7%). 


is kept for the following experimental studies. 


Finally with the selected values for percentage of 
boron and diameter, the influence of the power 
density P/S'is presented in Figure 4. A strange result 
is obtained: the energy thresholds of ignition firstly 
decrease for power densities lower than 2 kW cm 
and then increase for higher values. This result has 
been described in detail in ref. 10. Nevertheless it 
can be verified, in Figure 4, that a classical result 
for the B/KCIO, mixture is obtained. Moreover, 
the ignition thresholds of boron, if mixed with 
KCIO, oxidizer, are only one-third those of the B/ 
KNO, mixture in the same conditions. In the same 
way, we always obtained a regular decrease in the 
energy thresholds of ignition with power density 
on various pyrotechnic mixtures: Zr/KCIO,, Zr/ 
PbCrO,, Fe/KC1O,, Fe/KNO3. These pyrotechnic 
compounds were studied by Radenac* in our 
laboratory. 


In order to link the influence of all parameters, 
a numerical model has been completed. The 
following section presents the main _ results 
obtained with such a computing tool. 


Numerical simulation 


The principle of the numerical modeling is 
detailed in ref. 11. It relies upon a thermal balance 
of the pellet with a heat source coming from 
the chemical reaction (Arrhenius’ Law with nth 
order reaction). The heat exchange is considered 
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in the sapphire porthole. The description of 
the heat exchange at the interface between the 
pyrotechnic pellet and sapphire medium is taken 
into account by the equivalent thickness of air e 
and the depth 6 of absorption of the laser beam. 
In order to use this numerical model, input data 
are requested; most of them are obtained by means 
of thermal analysis measurements but two inputs 
are unknown and play an important role in the 
condition of ignition : The interface parameters e 
(i.e. the equivalent thickness of air at the frontier 
between the porthole/pellet and the depth 6 ~ yw"! 
of absorption of the laser beam. In the case of 
Fe/KCIO, mixtures, which have been previously 
studied in our laboratory [Radenac,’ ICT), these 
parameters were obtained in three different ways: 


1. analytical approach, 
2. optimization 


3. numerical modeling of the pile up of grains of 
reducer and oxidizer in the pellet with respect 
to the volume of air (porosity ¢) 


The first approach is presented in ref. 9 and gives 
the following results :e=27.1 um andoé = 18.7 wm 
(u = 5.3 10*m'') which are respectively the 


62 wt % Fe, 38 wt% KCIO,, porosity 26 % 
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equivalent thickness of air at the interface and the 
depth of penetration of the laser beam inside the 
volume of the pellet. 


By use of an optimization computing tool, by 
comparison with experimental data on Fe/KCIO, 
mixture, the following values have been found : 


e=8.5 um and uw =5.5 x 10* wm |. 


Using the numerical modeling of the pile-up of 
grains of reducer and oxidizer in the pellet, the 
following values have been found : e=9.5 um 


and uw =4.9 x 10* um. 


With the second set of values of interface 
parameters, the computation of the conditions 
of ignition of several pyrotechnic mixtures, 
previously studied in our laboratory was carried 
out. 


The main results are reported in Figure 5. It can 
be seen that the best agreement is obtained for Fe/ 
KCI1O,, and in this case the mean quadratic error 
between experiments and modeling is only 3%. 


For B/KCIO,, the fit is medium and the less 
acceptable mean quadratic error between 
experiments and modeling is 9%. 
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Figure 5. Comparison of numerical simulation and experiments on various pyrotechnic mixtures. 
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Table 1 Comparison between Bruceton experiment and two simulation tests (Fe/KCIO, mixture, P/S = 


1.53 kW cm ?). 


Time of energy threshold (ms) 


Experiment (Bruceton) 


= 5.25 
is 5.2 


A larger deviation is obtained in the case of the Zr/ 
PbCrO, mixture with a mean quadratic error of 
21%. The Zr/KCIO, mixture is represented with a 
medium error (15%). 


All these experiments were carried out with one 
set of interface parameters. This pair of values 
was optimized for the case of Fe/KCIO,. This is 
the reason why the best agreement between the 
numerical simulation and the experimental results 
is obtained for this mixture. The agreement between 
experiments and computation could be improved 
if, for each mixture, the interface parameters were 
determined by this method. 


But it is also possible to estimate the interface 
parameters by means of numerical simulation 
of the structure of the pellet. This simulation 
relies on a random simulation of stacking of 
grains of reducer and oxidant representative of a 


Calculation (Bruceton) Calculation (Langlie) 


6.2 6.54 
5.35 Di3d 
5:3 4.16 


pyrotechnic composition. The method is explained 
in ref. 12. This tool enables us to obtain the 
value of the interface parameter for any mixture 
without experiment (and experiment of ignition). 
In Figure 6, oxidant grains are in grey and with 
sizes in accordance with the real granulometry. 
Reducer grains are drawn in black; in this case the 
particle size is constant. The remainder volume (in 
white) corresponds to the volume of air, that is to 
say the porosity. All proportions of all components 
are respected. Thanks to this method, it is thus 
possible to obtain a statistical distribution of the 
interface parameters. 


Then, the statistical distributions of the interface 
parameters are used as input data in a numerical 
simulation program of the ignition in the 
condensed reactive system. The random aspect of 
the real conditions of an ignition of pyrotechnic 
composition is then introduced. The comparison 
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Figure 6. Numerical simulation of the structure of the pellet. 
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with the experiments authorizes choices in the 
manner of carrying out tests, in the field of 
pyrotechnics, in order to determine the probability 
of ignition under given conditions of energy 
deposition. 


Moreover, even for fixed parameters, the results of 
two consecutive calculations are never the same. 
This points up the influence of the heterogeneity 
of the mixture on the results. Thus it is necessary 
to carry out averages, 1000 tests on the standard 
mixture: 6 = 20.4 + 3.4 um and e =9.5 +2 um. 


It is then possible to carry out a virtual Bruceton (or 
Langlie) test because each numerical simulation 
is performed with a new pellet with it own values 
of interface parameters. 


The result of such a comparison is presented in 
Table 1. Different times which correspond to three 
energy thresholds (E£y;, Es) and E);) are obtained 
using experiment or the numerical method. 
Whatever the test used, the medium time fs, is in 
a very narrow range and close to the experimental 
value. The higher energy level is slightly higher 
in the case of the virtual Langlie test. The more 
significant discrepancy is found for the Ey); energy 
thresholds which is lower when computed by 
means of a virtual Langlie test. Our interpretation 
is that there is an over-estimation of Ey; energy 
thresholds when a Bruceton method is used. 


Conclusion 


A specific laser diode ignition system was 
presented in this paper. Some results in the case of 
boron mixtures were delivered. For this reducing 
agent, an abnormal behavior was observed when it 
was mixed up with KNO,; oxidizer. The evolution 
of the ignition thresholds as a function of the 
power density exhibits a minimum for a value of 
P/S of 2kW cm. This particular characteristic 
was never observed for the five other pyrotechnic 
mixtures previously studied in our laboratory. The 
detailed explanation was reported in ref. 10 and is 
based upon the presence of BN particles produced 
during the combustion which form an opaque 
cloud toward the laser beam. 


A specific numerical modeling of the ignition 
process has been completed. Among the various 
parameters used in this code, particular attention 
was paid to the interface parameters between the 
surface of the pellet and the optical window which 
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is in direct contact with the energetic material. 
The values of these parameters were found in 
three different ways: an analytical approach, 
an optimization process and, finally, a random 
simulation of the structure of the pellet. In this last 
case, the pile structure of the pellet was used in order 
to determine the value of the interface parameter 
for each shoot. It was then possible to obtain a 
given threshold of ignition which corresponded 
to a precise heat balance at the surface pellet. 
As it is possible to compute deviation of ignition 
thresholds around medium values F's, it is also 
possible to compute statistical law of probability 
of inflammation versus energy. The numerical 
modeling of ignition is a determined code. But ifa 
combination of this code with random simulation 
of the pellet structures is performed, pyrotechnic 
tests like Bruceton or Langlie are potentially 
achievable. Virtual pyrotechnic tests were 
executed and compared with our experiments. A 
very good agreement was observed for E's) and 
Eos thresholds. A less acceptable result for the Eo; 
energy of ignition is imputed to the Bruceton test 
itself. 


Inconclusion, itis possible, withthe numerical tools 
presented in this paper, to simulate conditions of 
ignition with a laser diode on various pyrotechnic 
mixtures. The validations were carried out on five 
different pyrotechnic compounds leading to a good 
agreement between experiment and calculations. 
The difference between these values was in the 
range 3 to 21%. This new method is potentially 
applicable to other ignition tests in the fields of 
pyrotechnic. 
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Synthesis, Structural and Reactive Characterization of 
Miscellaneous Nanothermites 


Denis Spitzer and Marc Comet 
Institut franco-allemand de Saint Louis (ISL) ; Laboratoire ISL / CNRS NS3E (FRE 3026) ; 68301 Saint Louis, 


Abstract: /n this paper, we report the elaboration and the characterization of several kinds of nanothermites 
made at the French German Research Institute of Saint-Louis (ISL). Three types of materials are presented 
herein: tungsten trioxide based nanothermites obtained by physical mixing of this metallic oxide with 
aluminum nanoparticles, tungsten trioxide based nanothermites elaborated by coating WO; nanoparticles 
with aluminum using a chemical process; and molybdenum oxide based nanothermites composed of 
Al,Mo,O, nanostructured phases and aluminum nanoparticles. In the light of these examples, we have 
identified general trends concerning the structural and reactive behavior of the new materials. 


Keywords: nanothermites, tungsten oxide, molybdenum oxide, aluminum coating 


Introduction 


Thermites are non-nuclear materials with a very 
high bulk energetic potential (10-25 kJ cm’). 
Thermite mixtures are composed of a metal oxide 
and an oxophile metal used as a strong reducer.' 
Contrary to explosive energetic compositions, 
thermites have a high density, and slow kinetics 
of combustion. Moreover, thermites have lower 
sensitivity than explosive energetic compositions 
regarding thermal and mechanical stress. Finally, 
they do not detonate, and generate small quantities 
of gaseous products. Decomposition of thermites 
is a redox reaction between two solid phases. The 
combustion kinetics are governed by the mass 
transfer between the oxidizer and the fuel phase. 
For this reason, nano-sized energetic compositions 
are of great interest because of the unconventional 
properties resulting from their high ratio of surface 
to volume atoms. A simplified model was used to 
determine the particle sizes for which the reactivity 
is expected to be changed.” From this model, we 
found that “nano-effects” appear when particles 
are smaller than some hundreds of nanometres. 
They are dramatically increased for particles 
having a diameter below twenty nanometres. 


To elaborate nanothermites having energetic 
performances which significantly differ from 
classical thermites, it is necessary to use metallic 
oxide and metal particles having sizes in the range 
defined above. 


As aluminum is the most efficient metal to 
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formulate thermites, aluminum nanoparticles 
are commonly used to elaborate nanothermites. 
These particles exhibit a “cherry” structure with 
an aluminum core surrounded by an alumina shell. 
The use of aluminum nanoparticles smaller than 
fifty nanometres to elaborate nanothermites is not 
reliable because their alumina content is too high. 
An elegant way to solve this problem consists 
in coating metallic oxide nanoparticles with 
aluminum deposited by a chemical process. 


Metallic oxide nanoparticles are easier to elaborate 
and to stabilize at the nano-scale. For instance, an 
original sol-gel process can be used to structure 
molybdenum oxide particles with different alumina 
contents. The alumina content is used to control 
the size and the nature of Al,Mo,O, particles and 
the reactivity of nanothermites containing these 
phases. 


Time resolved cinematography (TRC) is a 
method which is routinely used at ISL to study 
the combustion of energetic materials. From an 
experimental point of view, the energetic material 
is ignited by a laser beam (CO2; 10 W; 10.8um) 
which is focused by the mean of a lens disposed 
on an optical bench. The ignition delay time is the 
duration between the laser impact on the energetic 
material and the beginning of its combustion. 
The ignition delay time can be correlated with 
the density of surface energy necessary to initiate 
the material. The combustion rate is measured by 
ultra fast cinematography using a Photron camera 
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which can capture till 125,000 frames per second. 


Tungsten trioxide based nanothermites 
obtained by physical mixing* 


Tungsten trioxide (WO;) based nanothermites 
were elaborated by physical mixing of the oxide 
with aluminium nanoparticles. WO; and Al 
nanoparticles are commercially available and were 
purchased from Aldrich and Nanotechnologies 
respectively. From an experimental standpoint, the 
mixing was carried out by dispersing the particles 
in a liquid phase with ultrasonic stirring. In this 
case, diethyl ether was chosen as the dispersive 
medium because this solvent can be easily 
removed by evaporation. In addition, diethyl 
ether allows elimination of most of the adsorbed 
water which is trapped within the alumina shell 
covering the aluminum nanoparticles. Therefore, 
this method leads to dried nanothermite powders, 
which must be kept in a rigorously anhydrous 
atmosphere. As a general rule, adsorption water 
does not have a significant influence on the ageing 
of the aluminum nanoparticles, except if the 
material undergoes thermal shocks resulting in the 
condensation of atmospheric moisture. However, 
it was noticed that adsorbed water significantly 
desensitizes WO,;-based nanothermites towards 
friction stresses. 


Most of the publications dealing with the reactivity 
of nanothermites focus on the influence of the size 
of the aluminum nanoparticles. Although this 
approach seems attractive, it is limited by the fact 
that aluminum nanoparticles contain more and 
more alumina as they become smaller.‘ In thermite 
mixtures, alumina is considered as an inert 
compound which limits the diffusion of oxygen 
atoms from the metallic oxide to the aluminum. 
Hence, decreasing the aluminum particle size to 
improve the reactivity is not realistic for particles 
smaller than fifty nanometres. Therefore, we 
decided to use the same aluminum nanoparticles, 
having a 45.7 nm aluminum core surrounded by a 
2.9 nm alumina shell, in all experiments. 


Surprisingly the incidence of the size of the 
metallic oxide particles on the reactivity of these 
materials has been disregarded. However, it can be 
stressed that oxide tuning is the key to controlling 
the performances of nanothermites in the future. In 
order to demonstrate that the size of oxide particles 
plays a major role in nanothermite reactivity, three 
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0 x 1,000 


Figure 1. Thermites obtained by physical mixing 
of tungsten trioxide particles having a mean 
diameter of 30 nm (A.) and 20000 nm (B.) with 
aluminium nanoparticles (diameter ~ 50 nm). 


WO, grades were used. These materials, which 
were provided by Aldrich, had mean particle sizes 
of 30, 50 and 20000 nm respectively. 


The microstructure of WO,/Al nanothermites 
was observed by Scanning Electron Microscopy 
(SEM). In the mixtures containing only 
nanoparticles, it is difficult to discriminate 
particles according to their chemical nature 
(Figure 1A). Conversely, it appears clearly that 
aluminum nanoparticles homogeneously cover the 
surface of WO, micrometric particles (Figure 1B). 
These micrographs are a good illustration of the 
fact that the nanostructuring of each component 
is necessary to reach a satisfying contact between 
the reactants. 


WO,/Al nanothermites which contain only 
nanoparticles have an impressive reactivity. The 
fireball generated by the deflagration is so hot that a 
bang due to overpressure is heard. The combustion 
rate can reach 7.3 ms '. This value is extremely 
high compared to classical energetic materials. For 
example, a mixture of potassium chlorate and sugar 
tested in similar experimental conditions exhibits 
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Figure 2. Scanning electron micrographs of pure WO; nanoparticles (A.) and WO,/Al nanocomposite 
(B.). Transmission electron micrograph illustrating the coating of WO; nanoparticles by chemically 


deposited aluminum (C.). 


a combustion speed two orders of magnitude 
lower. The ignition delay time of some WO,/ 
Al nanothermites is less than two milliseconds 
with a CO, laser source (A = 10.6 tm) delivering 
a 10 W power on a surface of approximately 
1 mm’. However, dried WO,/Al nanothermites 
are extremely sensitive to friction stress (<4.9 N). 
This sensitivity is significantly reduced when the 
aluminum nanoparticles contain small amounts 
of adsorbed water. WO,/Al nanothermites are 
relatively insensitive to impact stress (42.2 J). TRC 
experiments showed that WO,/AI nanothermites 
aged at 200 °C in air for several hours keep their 
reactive properties. Nevertheless, the repeatability 
of energetic properties is deteriorated by thermal 
ageing. 

Tungsten trioxide based nanothermites 
obtained by chemical coating 


Tungsten trioxide based nanothermites can 
be remarkably desensitized by embedding the 
metallic oxide nanoparticles with aluminum. This 


process was patented by ISL. 


Scanning electron microscopy (SEM) reveals that 
pure WO; particles (Figure 2A) are significantly 
smaller than WO,/AI particles (Figure 2B). The 
coating was proved by transmission electron 
microscopy (TEM) which clearly shows. that 
aluminum is located around tungsten trioxide 
particles and covers them (Figure 2C). 


Nitrogen adsorption corroborates this result. 
Indeed the experimental specific area of the WO;/ 
Al composite is higher than the specific area of 
WO, particles. A mathematical approach allows 
demonstration that the density of aluminum 
prevails against the increase in diameter. For 
instance, in the case of a WO; spherical particle 
(diameter ~ 30 nm) which is covered by a 
homogeneous aluminum shell, the specific surface 
area will be smaller than that of WO; provided that 
the aluminum thickness does not exceed twenty 
nanometres. 


Figure 3. SEM observation of the microstructure of an agar/APM composite gel (A.); TEM picture 
representing the string-like structure of elementary strands (B.). 
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Figure 4. Qualitative (SEM) and quantitative 
(BET) characterization of the microstructure 
of Al.Mo,O, phases according to their alumina 
content: 0 wt% (a), 11.5 wt% (b), 17.0 wt% (c) 
and 23.3 wt% (d). 


During their thermal decomposition, these 
materials release the same amount of energy 
as the corresponding nanothermites elaborated 
by physical mixing. However, the sensitivity to 
mechanical stresses is significantly decreased. 
This can be explained by the fact that composite 
particles behave towards each other like metallic 
particles. 


Al,Mo,O, based nanothermites 5 


Mixed oxide phases (Al,Mo,O.) were structured 
at the nano-scale in order to get nanothermites 
with tuned reactive properties. For this purpose 
an ammonium paramolybdate (APM) aqueous 
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solution was gelified by using agar, a polysaccharide 
of botanic origin. The resulting hydrated gels 
were washed with acetone in order to remove 
water. Water extraction induced the solidification 
of the inorganic phase (APM) within the organic 
matrix (agar). These materials are highly porous 
crumbly monoliths called “composite gels”. They 
have a three-dimensional porous structure formed 
of interlinked strands of average section equal to 
45 nm (Figure 3A). APM particles are located 
within the agar strands, in which they define a 
string-like structure (Figure 3B). A mathematical 
approach based on specific area measurements was 
used to optimize the composite gel formulations. It 
results from these calculations that the agar content 
of an ideal composite gel must range between 40 
and 60 wt%. For these formulations, APM does 
not crystallize out of the agar matrix. Moreover, 
APM particles have small sizes (20-34 nm) and 
they are sufficiently spaced in the agar strand 
(2.0-3.5 nm). 


Composite gels underwent an oxidative treatment 
which allowed burning off the agar and conversion 
of APM into molybdenum oxide. A direct thermal 
oxidation of composite gels results in a structural 
collapse, leading to molybdenum trioxide (MoO;) 
micrometric particles. A chemical treatment 
was developed to prevent this phenomenon. It 
consists in treating composite gel powders with an 
anhydrous aluminum trichloride (AICI) solution 
in diethyl ether (Et,O). The water adsorbed by the 
composite gels readily reacts with AICI, to give 
an oxychloride layer at the surface of gel strands. 
During the thermal treatment, this layer prevents 
the collapse of the structure by acting as an in situ 
nano-sized crucible. Most of the chlorine atoms 
were removed and several chemical species were 
formed. The formula Al,Mo,O, refers to these 
species, which mainly consist of amorphous 
alumina associated with crystallized molybdenum 
dioxide (MoO,) and aluminum molybdate 
(Al,(MoO,);). The under-oxygenated MoO, phase 
probably forms in the reducing environment of the 
core of the agar strand. The aluminum molybdate 
most likely appears at the surface of the agar strand, 
where the APM particles are in contact with AICI, 
and oxygen from air. The final microstructure 
of Al.Mo,O, phases directly depends on the 
concentration of the AICI,Et,O solutions used 
to impregnate the composite gel (Figure 4). The 
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Figure 5. Evolution of the ignition delay time 

(6) and of the combustion rate (V,) of Al,Mo,O, 
nanothermites depending on the alumina content 
of the Al,Mo,O, phase. 


transition from micro- to nano-scale occurs for 
alumina weight contents higher than 8.5%. An 
alumina content of more than 17.0% results in the 
densification of the material. 


Nanothermites were elaborated by physical mixing 
of the Al,Mo,O, nano-sized phases with aluminum 
nanoparticles. From an experimental standpoint, 
the dispersion was performed by ultrasonic stirring 
in hexane. The powdery materials were pelletized 
and characterized by time resolved cinematography 
(Figure 5). Their ignition delay times (2.8—17.8 ms) 
and their combustion rate (1.7-103 cm s') can 
be adjusted through the formulation of Al,Mo,O, 
phases. It can be noticed that the reactivity of these 
nanothermites is correlated to the structure of the 
Al,.Mo,O, phases. Indeed, the higher combustion 
rate corresponds to the micro to nano transition 
while the smaller ignition delay time is observed 
just before the Al,.Mo,O, densification. 


Although they have better energetic performances, 
Al,.Mo,O, nanothermites are more insensitive 
to thermal and friction stresses than a thermite 
elaborated with pure micron-sized molybdenum 
trioxide. This insensitivity is probably due to the 
presence of alumina, which limits the contact 
between the reactive species. 
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Conclusions 


The energetic performances of nanothermites 
strongly depend on the elementary microstructure 
of the metallic oxide (WO ;, Al,Mo,O-) and the 
reducing metal (Al) they are made of. 


Concerning aluminum, the main obstacle is the 
“cherry” structure of the particles which limits 
the amount of useful matter and prevents contact 
with the metallic oxide. To solve this problem, 
the unique solution is to coat metallic oxide 
nanoparticles by chemically deposited aluminum. 


Concerning metallic oxides, a decrease of the 
WO; particle size dramatically shortens the 
ignition delay time and significantly accelerates 
the combustion rate of WO;/Al nanothermites. 
Alumina can be used to tune the structure of 
molybdenum oxides in the form of Al,Mo,O, 
phases. These nanostructured materials allow 
elaboration of highly insensitive Al,Mo,O,/Al 
nanothermites and control of their reactivity. 
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Abstract: HNF hydrazinium nitroformate, H;NNH,*+C(NO,)* , is a water soluble salt. It has a positive 
oxygen balance of 13%, a much more positive enthalpy of formation than ammonium perchlorate (AP) and a 
much higher heat of explosion than AP, 5579 J g‘ against 1972 J g '. Additionally it has no chlorine and all 
the problems with hydrogen chloride formation can be avoided when used as oxidizer in rocket propellant 
formulations. All these advantages together could make HNF an oxidizer with better performance than AP. 
One inherent disadvantage may be the lesser thermal stability of HNF. Therefore an extensive investigation 
was performed on the thermo-chemical stability of HNF: Three sample lots of HNF have been investigated 
at ICT. They were provided by APP BV, The Netherlands. The thermal stability was determined by the 
following methods: 


* autoignition temperature with 0.2 g at 5 °C min! heat rate in a Wood’s metal bath 
¢ vacuum stability test (VST) 

° heat generation rate as function of time and temperature 

¢ mass loss as function of time and temperature 

* adiabatic self heat rate. 


Lots 2 and 3 have been characterised by heat generation rate at 60 °C, 65 °C, 70 °C and 75 °C and in short 
by mass loss. Lot 1 was extensively used for mass loss determinations in the temperature range 50 to 80 °C. 
HNF shows high heat generation rates. All curves from both methods indicate self accelerating behaviour. 
They have been described with autocatalytic reaction kinetic models. The Arrhenius parameters have been 
determined for lot 1 from mass loss data and for lots 2 and 3 from heat generation data. The activation 
energies for the intrinsic decomposition reaction are 166, 139 and 132 kJ mol! and for the autocatalytic 
reaction 159, 128 and 117 kJ mol’ in the order lots 1, 2, 3. The kinetic data are compared and discussed. 
Data for lifetime at different temperatures are given in terms of the times to reach preset values of mass loss 
and energy loss. 


Keywords: autocatalytic, kinetic, HNF 


Introduction 


The chemical formula of HNF is shown in 
Figure 1. It is a salt between the nitroform acid 
and hydrazine. One may name the substance as 
hydrazinium trinitromethanate and by ignoring 
its salt character hydrazine trinitromethane. It 
is easily soluble in water. Especially at elevated 
temperatures the salt may be in an equilibrium 
state with its pre-dissociated form, shown in 
Figure 2, as is assumed also for ammonium nitrate 
(AN) and ammonium dinitramide (ADN). In 
Table 1 some data are compiled to characterize 
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and compare HNF with other oxidizers. 


HNF has been investigated by thermal analysis 
several times already.” * These papers concentrate 


HaN—NH3 =< i NOD 
\ O> 
Figure 1. Chemical formula of HNF, hydrazinium 


nitroformate. 
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Table 1. Performance data of the oxidizers ADN, AP, AN and HNF-. The values of the heats of explosion 
and gas volumes have been calculated by the ICT-Thermodynamic-Code. The other data are from the 


ICT-Thermochemical Data Base! 


Substance Molar mass/ O, balance Enthalpy of formation On Ie" Density/ Gas vol®/  Ty°/°C 
gmol ! (%) e/ gem~ cm? g! 
kJImol'! Jg! 
ADN 124.056 +25.8 —149.8 —1207.5 3337 1.81 592 92.9 
AP 117.489 +34.0 —295.8 —2517.7 1972 1.95 533 130; D 
AN 80.043 +20.0 —365.6 —4567.5 2479 1.73 459 169.9 
HNF 183.081 +13.1 —76.9 —420.0 5579 1.91 568 124;D 


* Ocy heat of explosion, value with reaction water as liquid (25 °C). ° Gas volume at 25 °C, 1 bar without water, for 
thermodynamically controlled combustion. © D in the column of 7}, (melting temp.) means decomposition at the given 
temperature, which may vary with the sensitivity of the observation. 


;? 
NO — ~ NO, 
H»N— NH3 ( 2 ——» | H2N iG 
ne NO 


Figure 2. Equilibrium between salt form and pre-dissociated form. 


mainly on decomposition products, decomposition 
mechanisms and safety aspects. In some data® on 
activation energies for the decomposition of HNF, 
measurements (TGA, DSC) are given obtained by 
TA, and typical TA data evaluation was performed. 
These TA models do not describe the data well. The 
activation energies given? are in the range of 217 
and 235 kJ mol obtained by thermal analysis (TA). 
In the present paper the data measured between 50 
and 80 °C have been described with autocatalytic 
models, which give very good reproductions of 
the isothermal measurement curves. 


Substances, measurement methods 
and conditions 


Three lots of HNF were shipped to ICT from the 
company APP (Aerospace Propulsion Products 
BV), RT Klundert, The Netherlands. The colour 
of all three samples was intensely yellow-orange. 
All lots were of type S according to the APP 
classification. The particle shape was different 
between HNF lots 1/2 and HNF lot 3. HNF lots 1/2 
were more needle-like and coarser than HNF lot 3, 
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which appeared nearly equal in all three geometric 
dimensions. The mean particle size of lot 1 was 
about 240 um with a length to diameter ratio of 
7.4, and for lot 3 these were found to be 76 um 
and 2.8. Table 2 lists some data. The HNF was 
used as delivered. The purity is, according to APP 
information, between 98.8 mass% and 99.6 mass% 
based on acid content (not further specified) and 
between 97.2 mass% and 99.6 mass% based on 
hydrazine content. To determine the thermal 
stability, several test methods and measuring 
methods were applied: 


* autoignition temperature at 5 °C min’ heat 
rate, 200 mg sample size; 


* vacuum stability test (VST), 2 g sample size; 


¢ heat generation rate (HGR, dQ/dt) as a 
function of time and temperature, determined 
with a TAM™ (Thermal Activity Monitor), 
produced by Thermometric AB, Sweden; 
3 ml ampoules, sample amount between 1.2 
and 1.4 g; 


HNF lot 1 75 °C 
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Table 2. Characteristic data determined with HNF lot 1. Experimental oxygen content as difference to 
100 mass% from measured CNH data. S means limit energy for impact (impact sensitivity); R means limit 


force for friction (friction sensitivity) 
Elemental composition [mass%]| 


Cc H N O 


Measured at 
ICT 7.01 2.78 37.64 


Theoretical 
values 


APP-data®’ — _ — — 


HINF let2. 60°C, 65 °C, 70 °C, 75 °C 
HNF lot 3. 60 °C, 65 °C, 70 °C, 75 °C 


e mass loss (ML) as function of time and 
temperature, determined in PID controlled 
and isothermally operated aluminium block 
ovens, sample amount | g and 2 g per glass; 


HNF lot 1 50°C, 60°C, 65 °C, 70 °C, 75 °C, 
80 °C 

HNF lot 2 65 °C, 70°C,'75 °C 

HNF lot 3 65 °C, 70 °C, 75°C 


¢ adiabatic self heat rate, determined with an 
ARC™ (Accelerating Rate Calorimeter), 
200 mg sample amount; 


With these methods two tests on the samples are 
achieved: 


* gas generation by VST and mass loss 


¢ neteffect of the sum of heats ofreactions during 
decomposition by autoignition temperature, 
heat generation rate, adiabatic self heating. 


These different investigations into the 
decomposition behaviour of a sample confirms the 
results, if they are consistent.® 


(52.57) 


6.560 2.753 38.254 52.433 


S/Nm RIN Density/ Heat of 
gcm° comb/J g! 
2 16 1.846 5796 
— — 1.87-1.93 (5579) 
2-5 18-36 1.86 5824 


The autoignition temperature was determined with 
a standardized apparature at 5 °C min | heat rate 
in a Wood’s metal bath. A sample mass of 0.2 g 
was used in special vials. The vacuum stability 
test, also a standardized test, was performed with 
the special conditions adapted to HNF by TNO- 
PML: test temperature 60 °C, test time 48 h. 
The sample amount was 2 g. The gas generation 
was determined in the vacuum apparatus with 
an integrated mercury column manometer. The 
last two columns of Table 3 show the values 
transformed from 60 °C, 48 h to the standard 
conditions of 90 °C, 40 hand 100 °C, 40h, whereby 
a factor of 4 for a 10 °C temperature change has 
to be used because of the values of activation 
energy found. The limiting value of gas generation 
volume under standard conditions is 1.2 ml g'' for 
both test conditions for the assessment 'stable'. 
Even at 90 °C this limit is far exceeded. For 80 °C 
no official test condition exists. Scaling the 90 °C 
conditions to 80 °C results in a limiting value of 
0.3 mlg¢!. 


Heat generation rate and heat generation 
Measurement data 


The heat generation rate (HGR, dQ/dt) was 
measured with a TAM™ microcalorimeter in the 


Results standard glass ampoules. Because HNF has a high 
Stability tests positive oxygen balance one may assume that 
Table 3. Stability test data. 
Autoignition temp./. VST [ml 7] Extrapolated Extrapolated Extrapolated 
HNF lot 7 A a4 3 g [ml g '] 80°C, [ml g '] 90 °C, [ml g¢ '] 100°C, 
Cat5°C min 60 °C, 48h, 2 g 40h 40h 40h 
HNF lot! 131 0.082 4.4 17.5 
HNF lot2 129 0.084 4.5 17.9 
HNF lot3 129 0.072 3.8 15.4 
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1200 


dQ/dt [yW/g] 


me 75°C, air, lot 1 


800 75°C, air 


ce | 70°C, air 


400 | 65°C, argon 
65°C, air 


200 


0 5 10 15 


HNF lot 2 


time [d] 


25 30 35 40 


Figure 3. HNF lot 2. Heat generation rate (HGR) as function of time at temperatures between 60 °C and 


75 °C. HNF lot 1 is included for comparison. 


measurements with air inside the ampoule should 
not greatly influence the HGR. Some measurements 
were also made in argon atmosphere. This was 
achieved by closing the ampoules in a high 
performance glove box from M. Braun Inertgas- 
Systeme GmbH, D-85748 Garching, Germany, 
type MB-200—MOD, operated with argon. The 
comparison of the curves indicates that the 


200 


Q [u/ 

180 | CW 
75°C, air, lot 1 

160 


140 75°C, air 70°C, air 


120 


100 70°C, argon 


assumption is valid, see Figure 3 for lot 2. HNF 
shows complicated decomposition behaviour 
after some initial period. This is found in the 
HGR curves in a significant way. After some time 
period a sudden acceleration in heat generation 
activity can be seen, Figure 3, lot 2 and Figure 5, 
lot 3. After a further period mitigation occurs and 
the curve again has a more smooth shape, but 


HNMF lot 2 


65°C, air 


65°C, argon 


0 5 10 15 


60°C, air 


time [d] 


25 30 35 40 


Figure 4. HNF lot 2. Heat generation (HG) as function of time between 60 °C and 75 °C. HNF lot 1 is 


included for comparison. 
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H 


dQ/dt [yW/g] 


70°C, argon 


iv 


70°C, air 


| 
| 75°C, air 


0 10 


HNF lot 3 


65°C, air 


15 


65°C, argon 60°C, air 


y 


time [d] 


20 25 30 


Figure 5. HNF lot 3. Heat generation rate (HGR) as function of time between 60 °C and 75 °C. The 
change of the atmosphere gives no consistent picture with lot 3. 


the HGR values then increase more rapidly. The 
interpretation is that the intermediate activity has 
produced further autocatalytically effective species 
or at least raised their concentration peakwise. 
This effect can be seen especially at temperatures 
causing a medium decomposition rate, mainly 
at 65 °C and 70 °C. But close inspection of the 
curves at 60 °C reveals this effect also. At 75 °C 
the decomposition rate is already so high that the 


180 


Q [J/g] 70°C, argon 
Vu 


160 


140 75°C, air Covel 


120 


100 


10 


described activity increase is flattened and hard to 
see. Remarkable are the very high values of the 
measured heat generation rates compared with 
values of a double base propellant, which range 
between 15 and 60 uW g | at 75 °C. Figures 4 and 
6 show the heat generation of the two HNF lots. 
Figure 7 compares lots 2 and 3 via heat generation 
(HG, Q). HNF lot 3 is, during the initial conversion 
period, always somewhat more stable than HNF 


HM lot 3 


65°C, argon 


15 


60°C, air 


time [d] 


20 25 30 


Figure 6. HNF lot 3. Heat generation (HG) as function of time between 60 °C and 75 °C. 
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t,T dA(t,T dM , (¢,7 L 
290] (S| tan) (P| a) 
Oa(t.) = (“AHpa) A(O) _ 
a _ (2) --(Me?) 
dt , dt J, dt r 
(3) 


O,.(6,T) =O, (t,T)/ Ox (tJ; 


k 
A+B —*> 2B+C+S (-AH, 4) 
intrinsic decomposition 


lot 2 but then this changes and HNF lot 3 becomes 
significantly the more unstable substance with 
higher decomposition rates. At temperatures above 
70 °C the difference diminishes because of the 
already very high decomposition rates. Because of 
the high HGR and high gas generation rates the 
probability of the glass ampoule breaking is high. 


Evaluation of heat generation data 


Equation (1) shows the connection between the 
quantities heat generation rate dQ,/dt, molar 
amount changing rate d4/dt and mass changing 


Q [/g] 


A(t,T) = A(t,T)/ A(0); 


A —15B+C+S (-AA,,) 


M,.(t,T)=M,(t,T)/M,(0) 


(4) 


autocatalytic decomposition 


rate dM,/dt of decomposing substance A. With 
Equation (2) the total amount of evolved heat 
Ox(t.) by the complete decomposition of A is 
given, using the molar heat of reaction (-AHR 4). 
The complete decomposition of A is achieved at 
infinite time according to a first order reaction. 
This is not a real world property and the time ¢, is 
introduced, at which the decomposition of A is no 
longer measurable. The reaction kinetic determined 
quantity O,(t.) is often not obtainable because the 
measurements may take too long a time or the 
suitable apparatus is not available. As reference 


comparison 
lot 2, full lines 
lot 3, broken lines , 


time [d] 


Figure 7. Comparison of HNF lot 2 and HNF lot 3. Heat generation (HG) as function of time at 


temperatures between 60 °C and 75 °C. 
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(Ae?) =-k,(T)-A(t,T)-k,(T)- A(T) BT) (5) 
[Seat-7)) = +k,(T)-(-AH, ,)- A(t,T) + ky(T)-(-AHy 9) A(t,T)- BOLT) 6) 
d 1 2 
[Seat 7D) = K(T)-(1- Og, (t,7)) +2 (T)-(1- Ox, (tT) -(F + Ou(t,7)) 7) 
(—AHR a) = (AHR ;) = (AHR >) (8) 
Ki (L) = ky (1) = ki (LT) = (ZT) ko (L) = ky (LT) = kx (1) =k, (T)- A(0) (9) 
O,(t,T) = OF, + Q,(te)-| 1 ky D) +k, (2) 
R(T) + ky (T)-exp((k\ (2) + (T))-t) 
(10) 
0,(,7)= OF, +0,(,)4 1 ky (T)+(F +1): (7) 
R(T) + (ky (T) + Fk (T))-exp (kx (7) + (F +1), (1))-1) 
k\(T) +k (1) 
t Ih, R(T) +k, (1) -exp((k\ (1) + kx) t) 
ky (1) +-kQ (LT 
wr (k\(T) +k, (7)) = 


R(T) +KA(D)-exp (kT) +2 (2))-t) 


(A, (7) +k, (T)) 


KR (T) + ky (T)-exp((k, (1) +i (T))-2) 


quantity the heat of explosion Qgx, is taken 
therefore, with the value from the thermodynamic 
calculation.' Equation (3) is the normalized form 
of equation (1) using the normalizations also given 
in equation (3). 


The data evaluation was done with reaction 
kinetic modelling described in detail already 
elsewhere.” Here the important relations are given 
in short. HNF decomposes autocatalytically. 
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The corresponding reaction scheme is shown in 
equation (4). Substance A decomposes into gases 
C, solids S and into an autocatalytically effective 
product B. In the second parallel reaction B reacts 
with A and accelerates its consumption as shown 
in equation (4). 


The intrinsic decomposition of A cannot be 
suppressed by stabilizers with their common 
function. To use stabilizers for this purpose would 
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require the 1:1 addition of a substance which is 
able to stabilize the electronic system of HNF by 
complex formation. With stabilizers of common 
function only the autocatalytic reaction can be 
influenced by removing B by chemically bonding 
it to the stabilizing substance. In equation (5) 
the intrinsic or inherent decomposition of A is 
included as a first order reaction. Reformulation to 
heat generation rates leads to equation (6).” 


Equation (6) is transformed into equation (7) used 
here by normalization with O,.¢= Q,(t.). B(0) must 
not be zero. This is considered in equation (7) by 
the factor F = B(0)/A(O). The effect of F # 0 is 
to emphasize the autocatalytic channel. This can 
change the course of the measured curve strongly.” 
In the following F is set to zero here. If the two 
heats of reaction in equation (4) are set equal then, 
with equations (8) and (9), equation (10) results. 
The equation used here for the modelling of the 
heat generation is this equation (10) with F = 0 as 
argued above. For completeness the corresponding 
equation is given with F 4 0 in equation (10). It 
is named model 'O: first order + autocatalytic'. To 
calculate the heat generation rate with the reaction 
rate constants obtained, equation (11) is used. 
Notice the equivalence of the rate constants with 
different indexing, see equation (9). 


To calculate the times ¢,p;,, to reach a certain 
energy loss EL by decomposition in substance A 
equation (12) is used in the case of autocatalytic 
decomposition. For completeness both versions 


are given, for F = 0 and F' # 0. The corresponding 
quantities are explained as follows. 


On ret = Qal(t.,T) is maximum energy of the sample 
(here Oa(t.) = Ogx.a = 5580 J g') 


Oa(t,o,a(7)) is amount of energy released by the 
sample A 


Os(tT) — QOsltox(T)) is momentary energy 
content of the sample A 


Von = (QalteT) — Onlhoa(D))) / Qa rer 


typ,a() is the time to reach yp, 4= 100% (1 —yo.a); 
BELA 40,4 

If only a first order decomposition reaction is 
active, equation (13) applies to calculate the times 
tyet,a- A smaller Ogx 4 value does not change the 
E,0; Values; it changes Zo;. 


Figure 8 shows an example of the modelling of 
heat generation with HNF lot 2. The obtained 
reaction rate constants are compiled in Tables 4 
and 5 for lots 2 and 3, together with the Arrhenius 
parameters. The corresponding Arrhenius plots 
can be seen in Figures 9 and 10 for HNF lot 2 and 
Figures 11 and 12 for HNF lot 3. The Arrhenius 
parameters are given for measurements in air only 
and for the combined measured data ‘air + argon’. 
Again the description of the data with the model 
is very good as can be seen from the correlation 
coefficients for the individual temperatures. 


| ky (L) +k (T))— ka (1) 
7 1 Yaa 
_ (= 2 “In 1 
YELA a (T) + ki (T) ky (T) 
(12) 
7 ky (1) +(F +1)-kq(1))- A(T) 
1 QA 
woes aan ; 
VELA kK (T)+(F+)-ky (7) kK, (T)+F ky, (T) 
1 1 
nak ) ki (T) {2 ( 3) 
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for first order decomposition only. 
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Table 4. HNF lot I and lot 2. Reaction rate constants from autocatalytic modelling of heat generation 
HG, with O,.¢= Ory = 5580 J g', air and argon atmosphere 


Temp./°C 
60 
65 
70 
75 
75 
65 
70 


ko. [1/d] 
1.4627 x 1044+1.8 x 107 
2.0876 x 1074+2.9 x 107 
5.2097 x 107+1.1 x 10° 
1.2042 x 107+2.0 x 10° 
1.0887 x 107 +3.3 x 10° 
1.9363 x 104+5.7 x 107 
6.3628 x 107+ 1.0 x 10° 


ko» [1/d] 


R2 


9.0094 x 107+ 6.3 x 10° 
1.7909 x 107! + 1.3 x 10% 
2.9501 x 10'+3.9 x 107 
7.0375 x 101+47.1 x 10% 
7.7276 x 10'41.2 x 10° 
1.7961 x 101 43.1 x 104 
2.8351 x 107+3.5 x 107+ 


Arrhenius parameters from measurements in air and argon 


E,;/kJ mot ! 
Ig(Z; [1/d]) 
R?, 


E,;/kJ mol! 
Ig(Z; [1/d]) 


146.2 + 14 132.445 
18.98 + 2.2 19.69 + 1.4 
0.978 0.988 
Arrhenius parameters from measurements in air only 
139.3 + 18 128.4+ 10 
17.94 + 2.8 19.08 + 1.6 
0.983 0.994 


R’, 


1 


Mass loss 


Measurement data 


Because of the limited amount of substance for 


0.99951 
0.99969 
0.99956 
0.99989 
0.99989 
0.99853 
0.99973 


Atmosphere 


air 
air 
air 
air 
air 
argon 
argon 


Lot 


NU NF NY NN NY 


HNF lot | only mass loss measurements could be 
performed in the necessary extension. Figures 13 
to 17 show the mass loss as function of time for 


all three lots. Figure 13 gives an overview of the 


Table 5. HNF lot 3. Reaction rate constants from autocatalytic modelling of heat generation HG, with 
Orep= Ory = 5580 J g’, air and argon atmosphere 


Temp./°C 


60 
60 
65 
65 
65 
70 
70 
70 
75 
75 


E,i/kJ mol | 
Ig(Z; [1/d]) 
R?, 


E,;/kJ mot ! 
1g(Z; [1/d]) 


kg, [1/d] kao [1/d] 
1.2965 107+4.4x107 1.1702 x 10° 
1.0634 x 10%+41.4x107 1.2863 x 10° 
2.4849 x 107410107 2.5756 x 10° 
1.6623 x 107445107 2.5640 x 10° 
1.5178 x 10%+4.2x107 2.3540 x 10° 
4.2145= 107+5.7x107 3.9323 x 10° 
3.7576 x 107418 10° 4.5477 x 10° 
3.3413 x 10741.2x10° 4.6363 x 10° 
8.6541 x 10442.2x 10° 8.0396 x 10° 
1.0330 x 10744.2x 10° 7.4789 x 10° 
Arrhenius parameters from measurements in air and argon 
134.2 + 13 117.144 
17.05 + 1.9 17.46 + 0.6 
0.967 0.995 
Arrhenius parameters from measurements in air only 
132.4+7 116.6+5 
16.82 + 1.1 17.38 + 0.8 
0.994 0.997 


R’, 
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b2.1x 104 
t 7.9 E-05 

t 4.6 x 10% 
+ 3.1 x 107 
3.3 x 10% 
b2.6 x 104 
+ 92x10" 
t 7.0 x 107 
t1.1x10° 
E1.7x 103 


R? 
0.99724 
0.99478 
0.99835 
0.99925 
0.99909 
0.99987 
0.99873 
0.99928 
0.99976 
0.99936 


Atmosphere 


air 
air 
air 
argon 
argon 
air 
argon 
argon 
air 
air 
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Q [J/g] 


HNF-S13, 1/2002 

75°C 

air 

autocatalytic modelling of heat generation 

Q(t) = OF +Qe*(1-(k, +k,)/(k,+k, *exp((k,+k,)*t))) 
OF =0 


Qe = 5580 J/g 
k, = 1.2042 E-3 + 2.0 E-6 1/d — — measurement 


k, = 7.0375 E-1 + 7.1 E-4 1/d model 
R’ = 0.99989 


time [d] 


Figure 8. HNF lot 2: Autocatalytic modelling of the heat generation at 75 °C. 


[Aa | = — Ky (T)- Mg, (7) — key (2): M 9, (7)-(F +1-M a, (t,7)) (14) 
kK (T)=k(T) and k2,(T)=k,(T)-A(0) and a (15) 
MOONE Dba eae ae (T))-1) 7 
ML(,7) = OF a, +100%- OED = OF , +100%-(1-M,(t,7)) (17) 
M(t,T)=1 ser M,,(t,T)) (18) 
ML(,T) = OF, +100%- ©. 1 ko (1) + an (P) (19) 


My | ke (7) + kin (T)-exp (kin. (7) + ku. ()) 1) 
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In(k, [1/d]) 


-9.5 


0.00285 0.00288 0.00291 


0.00294 


HNF-S13, 1/2002 

60°C - 75°C, only air 

modelling of heat generation 
intrinsic reaction rate constants 
Ea, = 139.3 + 18.2 kJ/mol 


Z, = 8.705 E+17 1/d 


In(Z,[1/d]) = 17.940 + 2.8 
R =-0.98336 


4/T [1/K] 


0.00297 0.00300 0.00303 


Figure 9. HNF lot 2: Arrhenius plot of the reaction rate constant of the intrinsic decomposition of HNF, 
obtained by autocatalytic modelling of the heat generation. 


ML data of lot 1 and Figure 14 details the data for 
the lower temperatures. The strong autocatalytic 
decomposition of HNF can be seen in Figure 15. 
HNF decomposes nearly completely into gaseous 
products at 80 °C. But at lower temperatures it was 
also found that HNF decomposes nearly completely 
into gaseous products. From Figure 15 it is clear 
that the course of the complete decomposition 


In(k, (1/d]) 


0.00285 


0.00288 0.00291 


0.00294 


of HNF is complicated. At least a two regime 
autocatalytic acceleration must be used to describe 
the complete decomposition. A self heating effect 
during the strong increase of mass loss may 
increase the decomposition rate additionally. 
Figure 16 shows the mass loss measurements up 
to 5% for HNF lot 2. The decomposition rate at 
75 °C is somewhat higher with lot 2 than with lot 


HNF-S13, 1/2002 

60°C - 75°C, only air 

modelling of heat generation 
autocatalytic reaction rate constants 
Ea, = 128.4 + 10.2 kJ/mol 

Z, = 1.208 E+19 1/d 


In(Z,[1/d]) = 19.082 + 1.6 
R =-0.98336 


0.00297 0.00300 0.00303 


Figure 10. HNF lot 2: Arrhenius plot of the reaction rate constant of the autocatalytic decomposition of 
HINF, obtained by autocatalytic modelling of the heat generation. 
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In(k,[1/d]) 
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0.00285 0.00288 0.00291 
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HNF-S19, 1/2002 

60°C - 75°C, only air 

modelling of heat generation 
intrinsic reaction rate constants 


Ea, = 132.4 + 7.1 kJ/mol 
Z, = 6.669 E+16 1/d 
In(Z,[4/d]) = 16.824 + 1.1 
R =-0.99427 


4/T [1/K] 


0.00297 0.00300 0.00303 


Figure 11. HNF lot 3: Arrhenius plot of the reaction rate constant of the intrinsic decomposition of HNF, 
obtained by autocatalytic modelling of the heat generation. 


1. Figure 17 shows the mass loss measurements up 
to 5% for HNF lot 3. Here also the decomposition 
rate is higher for lot 3 than for lot 1, especially 
after an initial period. There lot 3 seems somewhat 
more stable. 


Evaluation of mass loss data 


The data evaluation was done with reaction kinetic 


In(k, [1/d]) 
fick & 


0.0 


-1.0 


-1.5 


-2.5 


0.00285 0.00288 0.00291 


0.00294 


modelling described in detail already elsewhere? 
and in short below. Here also the reaction scheme 
equation (4) with the corresponding rate equation, 
equation (5), is used. Equation (5) is converted to 
the description with masses whereby M,,(1)=M,(t)/ 
M,(0). The result is equation (14) with the reaction 
rate constants and F’ given in equation (15). The 
integration of equation (14) yields equation (16), 


HNF-S19, 1/2002 

60°C - 75°C, only air 

modelling of heat generation 
autocatalytic reaction rate constants 


Ea, = 116.6 + 4.9 kJ/mol 
Z, = 2.398 E+17 1/d 
In(Z,[1/d]) = 17.380 + 0.76 
R =-0.99647 


4/T [1/K] 


0.00297 0.00300 0.00303 


Figure 12. HNF lot 3: Arrhenius plot of the reaction rate constant of the autocatalytic decomposition of 
HINF, obtained by autocatalytic modelling of the heat generation. 
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1 


5 (ba) +R) Ki 


ne Oo OE RD) by) G0) 
= M(,,, (T)) = _ ML 
MQ) =M,(t, (T))=1 apnee (21) 


which is used for the autocatalytic modelling of the 
mass loss data. Because HNF decomposes nearly 
without residue, the summarized molar mass mc 
of the gaseous species equals the molar mass m, 
of the decomposing substance A, which is here 
HNF. With freshly prepared substances one may 
assume F’= 0. Equation (17) stands for the actually 
measured mass loss and the normalized mass 
M(t,T) is given by equation (18), which shows the 
connection of mass loss M,,(t,7T) of investigated 
substance A with the measured quantity M,(t,7). 
Together with the assumptions about m;/m, and 
F equation (19) results, which was used for the 
kinetic description of the measured data. The 
quantity OFML represents an offset not caused by 
decomposition of substance A. Figure 18 shows 
an example of the data modelling. Equation (20) 


ML [%] 


is analogous to equation (12) to calculate the 
times ty, to reach preset mass losses ML in the 
autocatalytic decomposition case. Equation (21) 
gives the relation between degree of mass change 
Ym and mass loss ML and yyr. 


As said above HNF shows a_ complicated 
decomposition course. But for the assessment of 
the in-service time (often named ‘lifetime’) of 
HNF the decomposition up to 5% ML is enough, 
because the performance data are after 5% out 
of tolerance already. During this decomposition 
range the curve is homogeneous’ which means 
the reaction scheme of equation (4) is applicable, 
which is sustained by the good description of 
the data, see Figures 18a and 18b. The reaction 
rate constants of the intrinsic and autocatalytic 


HNF lot 1 
—* 80°C 
—B- 75°C 
—&- 70°C 
—K- 65°C 
—@- 60°C 
—&- 50°C 


Figure 13. HNF lot 1: Mass loss as function of time at the temperatures 50 °C, 60 °C, 65 °C, 70 °C, 


75 °C and 80 °C. 
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ML [%] HNF lot 1 


—*K- 65°C 
—*e- 60°C 
—48— 50°C 


0 50 100 150 200 250 300 350 


time [d] 


400 450 500 550 600 650 700 


Figure 14. HNF lot 1: Same as Figure 13, but coordinate scaling adapted to the curves at 50 °C, 60 °C 


and 65 °C. 


reaction are compiled in Table 8. There also the 
determined Arrhenius parameters of these two 
rate constants are given. The description of the 
data is always very good, recognizable with the 
high correlation coefficients. The Figures 19 and 
20 show the Arrhenius plots of the two reaction 
rate constants. In Table 9 the times to reach preset 
mass loss values are given for the measurement 
temperatures, calculated with the reaction rate 


100 
99 | ML [%] 


constants via equation (20). The times tyy,(7) 
can be determined by calculation of the two 
rate constants via their Arrhenius parameters. 


Discussion 
Comparison of the three HNF lots 


Table 10 shows all Arrhenius parameters obtained 
by reaction kinetic modelling. The activation 
energies from heat generation data are somewhat 


HMF lot 1 
80°C 


—* Mittel 
—* 80/1 
—*— 80/2 


time [d] 


0 2 4 6 8 10 12 14 16 18 20 


Figure 15. HNF lot 1: Mass loss at 80 °C up to complete decomposition. HNF forms only a very small 


residue in slow thermal decomposition. 


Page 78 


Journal of Pyrotechnics, Issue 26, Winter 2007 


lot 1, 75°C 


10 


HM lot 2 


12 14 16 


Figure 16. HNF lot 2: Mass loss at 65 °C, 70 °C and 75 °C, compared with lot I at 75 °C. 


smaller than from mass loss. The interpretation 
could be that with HNF the two types of probing 
the decomposition reactions do not go in parallel 
as found with other systems.’ The activation 
energies of HNF lot 2 and HNF lot 3 are on 
average somewhat different, but the values are 
within the error margins in congruence. In 
Tables 11 and 12 some characteristic data of the 
decomposition behaviour are listed, all at 75 °C. 


10 


The initial mass loss rates and energy loss rates 
differ from substance to substance. HNF lot 3 has 
the lowest initial rates. But in terms of mass loss 
it is the substance with the fastest decomposition 
rate, Table 13, followed by HNF lot 2 and then 
HNF lot 1. 


In Table 14 the times to reach energy loss (EL) 
data are presented. To reach the same conversion 
expressed as ML or as EL needs significantly 


time [d] 


12 14 16 18 20 


Figure 17. HNF lot3: Mass loss at 65 °C, 70 °C and 75 °C, compared with lot I at 75 °C. 
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ML [%] 


HNF S13 8/2000 

75°C 

autocatalytic modelling of mass loss 
OF = 0.07 + 0.03 % 

k, = 1.242 E-3+7.3 E-5 1/d 

k, = 2.924 E-1+9.5E-3 1/d 


R® = 0.9989 


time [d] 


Figure 18a. HNF Jot 1: Autocatalytic modelling of the mass loss at 75 °C. 


less time as EL. This is surprising but confirms 
the above made assumption that the two criteria, 
evolution of gases and net sum of the heats of 
reaction, do not go in parallel with HNF. To solve 
this discrepancy a third method must be employed. 
The best is to follow the decomposition behaviour 
by detailed analysis of decomposition gases and 
of undecomposed HNF. Figure 21 shows the 
times 1,1 to preset energy losses for lots 2 and 3. 


ML [%] 


HNF S13 8/2000 

60°C 

autocatalytic modelling of mass loss 
OF =0.12 + 0.02 % 

k, = 5.872 E-5 + 3.5 E-6 1/d 

k, = 2.137 E-2 + 6.4 E-4 1/d 


2 


R? = 0.995 


To reach 1% EL at 30 °C needs about 4.2 years 
for lot 3 and 7.4 years for lot 2, see Table 15. 
In Figure 22 the comparison of the two kinetic 
descriptions can be seen. It shows the times to 
1% EL with autocatalytic description and only 
first order description. The times of the first order 
description are longer by a temperature dependent 
factor between 3.1 and 4.5. 


Figure 18b. HNF lot 1: Autocatalytic modelling of the mass loss at 60 °C. 
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In(k, [1/d]) 


0.00280 0.00285 0.00290 


0.00295 


HNF-S13, 8/2000 

50°C - 80°C. 

modelling of mass loss 
intrinsic reaction rate constants 
Ea, = 166.2 + 9.9 kJ/mol 

Z, = 8.229 E+21 1/d 


In(Z, [1/d])= 21.915 + 1.52 
R =-0.98336 


0.00300 0.00305 0.00310 0.00315 


Figure 19. HNF lot 1: Arrhenius plot of the reaction rate constant of the intrinsic decomposition of HNF, 


obtained by autocatalytic modelling of the mass loss. 


Comparison of HNF with other energetic 
material 


In order to assess HNF as a possible ingredient 
in rocket propellants for military use, it is useful 
to compare it with other ingredients. Figure 23 
compares the mass loss behaviour of HNF with that 
of unstabilized nitrocellulose (NC). This standard 


0 Ink, [1/d]) 


-7 
0.00280 0.00285 0.00290 


0.00295 


material with nitrogen content of 13.15% shows 
better stability than HNF. The adiabatic self heat 
rate, Figure 24, provides a similar insight. HNF is, 
in the series HNF < ADN < NC < CL20 < RDX < 
HMxX, the most unstable substance. In comparison 
with ADN at lower temperatures HNF is less 
stable than ADN, see Figure 25, which compares 
the mass loss curves of HNF lot 1 and ADN ICT 


HNF-S13, 8/2000 

50°C - 80°C 

modelling of mass loss 
autocatalytic reaction rate constants 
Ea, = 158.6 + 5.3 kJ/mol 

Z, = 1.589 E+23 1/d 


In(Z, [1/d]) = 23.201 + 0.81 


R =-0.9978 


60°C 


4/T [1/K] 


0.00300 0.00305 0.00310 0.00315 


Figure 20. HNF lot 1: Arrhenius plot of the reaction rate constant of the autocatalytic decomposition of 
HINF, obtained by autocatalytic modelling of the mass loss. 
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10000 + ; 27.379 
F a oy ss lifetime from heat generation HG time [a] to 
br times to energy loss EL % EL 
[ comparison HNF lot 2 and 3 
1000 + | 2.7379 
t lot 2, 5% 
r —Lot 2, 1% EL 
i —- Lot 3, 1% EL lot 3, 1% 
100 + — Lot 2, 5% EL 4 0.2738 
c —- Lot 3, 5% EL 
lot 3, 5% 
10+ + 0.0274 
T PC] 
ri \ 1 L L 1 1 L L n 1 1 n 0.0027 


Figure 21. Times tyEL to reach preset energy losses EL determined from the measured heat generation 
data of HNF lot 2 and 3. 


at 70 °C and 75 °C. It is to be noticed that the 
so-called anomalous decomposition behaviour 
of ADN at around 60 °C only appears at special 
conditions of ADN purity. In particular, the water 
content seems to have an important effect.'? In 


100000 


HNF Lot 2 


10000 


80.2a/17.9a 
1000 


first order + autocatalytic 
100 calculation 


——1% EL, autocatalytic 
10 —1% EL, only 1st Order 
——time factor 


most cases these conditions are not encountered 
with ADN in practice and in use one can adjust the 
water content in such a way that the anomalous 
decomposition of ADN will not be effective. 


time factor [-] 


factor in time 
‘first order' / ‘first order + 


only first order 


9.2d / 3d 


Figure 22. Comparison of times t,, to reach preset energy loss EL determined with the autocatalytic 
model and with the reaction of first order only. The times determined by the last description are longer by 


a factor of about 3.1 to 4.5, depending temperature. 
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Table 8. HNF lot 1: Reaction rate constants obtained from autocatalytic modelling of mass loss up to 5%. 


The offset found equals in average nearly to the water content found at ICT, namely 0.18 mass%. 


Temp./°C 
80 
75 
70 
65 
60 
50 


kya (1/d] 

2.3120 x 10 34 
1.2420 x 10°34 
3.0670 x 10 *4 
1.6060 x 10+ 
5.8720 x 10° 4 


b2.8x 104 
E7.3x 10° 
E11 10° 
t 6.8 x 10° 
t3.5x 10° 


1.3778 x 10° 4 


Arrhenius parameters 


E,;/kJ mol! 
Ig(Z; [1/d]) 
R’; 


166.2 + 10 
21.92 + 1.5 
0.993 


E75 * 107 


kyn2 U/d] 


4.6470 x 10°14 
2.9240 x 10°74 
1.3300 x 1014 
5.1900 x 107 4 
2.1370 x 10° 4 
3.5627 x 1034 


158.645 
23.20 + 0.8 
0.998 


Offset [%] R 
t3.4x 107 0.19 + 0.07 0.997 
t9.5 x 103 0.07 + 0.03 0.999 
t2.x 103 0.18 + 0.03 0.999 
H1.1x 10° 0.11 +0.02 0.998 
16.4 x 104 0.12 + 0.02 0.995 
t1.2x 104 0.17 + 0.02 0.992 


Table 9. Times in days to reach preset ML at different temperatures for HNF lot 1, calculated with the 
corresponding reaction rate constants obtained by autocatalytic modelling, using equ\tion (19) and equ\ 
tion (20). The offset, caused by residual water not included. 


ML [%] 


nA WN 


50 °C [d] 
360 
514 
615 
751 


60 °C [d] 


72.1 
99.6 
117.1 
140.3 


Table 10. Arrhenius parameters compiled for overview. 


HNF lot Atmosphere Type 


lot 1 air 
lot 2 


lot 3 


air 
air 


ML 
HG 
HG 


Temp. 


range/°C_ kJ mol! 


50-80 
60-75 
60-75 


65°C [d] 70°C[d} —- 75°C [d] 80°C [d] 
27.9 12.65 4.15 2.38 
39.0 17.16 6.00 3.50 
46.1 20.03 7.21 4.24 
55.6 23.80 8.85 5.26 
E E 
at! Ig(Zi (dl) Ry Ag Ed)? 
166.2410 2192415 0.993 158645 2320408 0.998 
139.3418 1794428 0.983 1284410 19.08+1.6 0.994 
132.447  16.824£1.1 0.994 116.645 17.38+0.8 0.997 


Table 11. Characteristic decomposition data at 75 °C from mass loss. MLR means mass loss rate. 


HNF lot [o/dj at 75 °C [I/d] 
HNF lot! 0.124 1.242 x 103 
HNF lot2 0.081 8.115 x 10-4 
HNF lot 3 0.072 7.193 x 104 


Comparison of obtained kinetic data with 
literature data 


HNF was investigated years ago with regard 
to thermal stability and decomposition. Such 
literature will be taken here which deals with the 
determination of decomposition parameters. An 
early work is that of Koroban eral.'' He investigated 
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Initial MLR at 75°C Intrinsic reaction rate constant ky! 


Autocatalytic reaction rate constant ky)? 
at 75 °C [1/d] 

2.924 x 10°! 

4.914x 10! 

6.418 x 10"! 


the thermal decomposition of HNF between 70 
and 100 °C by gas generation and has determined 
kinetic parameters. Also Brill'* has investigated 
HNF with the fast heating method (T-jump) 
together with FTIR analyses of decomposition 
gases in the range 130 to 400 °C. Further studies 
containing decomposition parameters are from 
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Table 12. Characteristic decomposition data at 75 °C from heat generation, measured in air. HGR means 
heat generation rate, ELR means energy loss rate determined with Q,.¢= Ory = 5580 J g'. 


Initial HGR and ELR at 75 °C 


Intrinsic reaction rate 


Autocatalytic reaction rate 


FINE Tot constant kp! at 75 °C [1/d] constant kp” at 75 °C [1/d] 
[MW g"] [%/d] - 7 

HNF lot 1 10:32 0.109 1.089 x 10° 7.728 * 10" 

HNF lot2 74.08 0.115 1.147 x 103 7.383 x 10" 

HNF lot 3 61.30 0.095 9.492 x 10+ 7.759 x 10! 


Table 13. Times in days to reach given mass loss ML values at 75 °C for the three HNF lots, calculated 
with reaction rate constants obtained by autocatalytic modelling. The offset is excluded. 


ML [%] Lot 1 [d] 
1 4.15 
6.00 

3 7.21 

5 8.85 

10 11.26 


TNO*», from FOI’? and from Sinditskii et al. by 
flame front measurements. !> 


Figure 26 compares the ML data of this work with 
the gas generation data of Koroban.'! The curves 
look very similar. Also with gas generation the 
clear autocatalytic behaviour was observed. In 
Figure 27 the evaluation of Koroban’s data can be 
seen up to a conversion of 18%, whereby the total 


HNF, 75°C 


HNF, 7Q°C 


0 5 10 15 20 


Lot 2 [d] Lot 3 [d] 
3.99 3.59 
5.27 4.60 
6.06 509 
7.10 6.03 
8.59 7.17 


decomposition gas amount GG(¢,) was obtained 
from Figure 26, right side. The description 
of the GG data is very good. In Figure 28 the 
measurements of heat generation rate on HNF 
at 60 °C are compared. The left side shows data 
from this work, the right side from TNO," using an 
HNF from APP. That lot shows higher HGR than 
the lots used in this work. The last comparison 
between measurements of the two laboratories 


HMF lot 1 
compared with NC 


—-75°C 
—*-—70°C 
—EF 50°C 
—G-NC, 75°C 
—H-NC, 70°C 


I, — NC, 70°C 


HNF, 50°C 


25 30 35 40 45 50 


Figure 23. Comparison between nitrocellulose (N content = 13.15%) and HNF lot 1 by mass loss 


measurements. 
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Table 14. Times in days to reach given energy loss EL values at 75 °C for the three HNF lots, calculated 
with reaction rate constants obtained by autocatalytic modelling. Oye = Ory = 5580 J g '. Measurements 


in air. 
EL [%] Lot 1 [d] 
1 2.72 
2 3.54 
3 4.05 
5 4.71 
10 5.66 


is made in Figure 29. Here also heat generation 
rates are compared obtained in this work and at 
FOI, Sweden.’” The data are comparable; the HNF 
from FOI seems somewhat more stable than the 
lots used here. 


Table 16 lists all Arrhenius parameters determined 
in this work and found in the cited literature. 
The data*'* are not in the range found from the 
other measurements. The parameters from gas 
generation'' fit with data obtained from mass loss 
in the temperature range 50 to 80 °C. In the second 
last row of Table 16 autocatalytic description 
parameters are given, calculated from Figure 3, 
C,*k, inref. 11, as was done for the initial first order 
decomposition. In the last row Koroban's data are 
given as cited in his paper’! for the autocatalytic 


h [°C/min] 


HNF lot 1 


Lot 2 [d] Lot 3 [d] 
2.73 2.86 
3.58 3.70 
4.11 4.21 
4.81 4.87 
5.80 5.82 


part. But for the first order or initial part this new 
set of £,, and Z, was taken from Koroban’s notes 
by Sinditskii and communicated to the author of 
this paper. The evaluated heat generation data of 
four lots are in agreement with each other within 
the limits of accuracy. 


The parameters obtained from flame front 
measurements by Sinditskii'® are a special 
case. The activation energy given in the row 
with E, = 140.7 kJ mol! was calculated using 
equation (22) with data on TS and 7,(TS) obtained 
from Sinditski1. Equation (22) is an approximation 
for equation (23), which gives the results in the 
row below with E, = 135.2 kJ mol '. These results, 
extractable from Sinditskii’s data,'® fit well with 
the determined data from heat generation. In 


adiabatic self heating 


Figure 24. Adiabatic self heat rate h as function of the adiabatically reached temperature of HNF, ADN, 
NC and some high explosives. Measurements performed with an ARC™. Sample mass between 200 mg 


and 300 mg. 
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ADN, 75°C 


HMF lot 1 
compared with ADN ICT 


HNF, 70°C 


ADN, 70°C 


Figure 25. Comparison by mass loss measurements of HNF lot 1 with un-stabilized ADN manufactured at 


ICT: 


equation (23) a lot of material properties are needed 
and data on the heat of reaction Qp are often not 
easy to get. Information about the derivation of 
equation (23) can be found for example around 
equation (11).'” 


Certainly the activation energy £, is oftena complex 
quantity, especially the one determined from HGR 
or HG. But when the measurement method and 
the evaluation procedure have enough ‘resolution’ 
to separate the individual kinetic processes one 
can identify the activation energy of the limiting 
decomposition step of the material. Here one 
could assume that in the hot surface of the burning 


strand the activation energy found is assigned to 
the gasification process. This means one has with 
such methods high material conversion in the 
solid layer. The resulting description by Arrhenius 
parameters gives for E, a mix of intrinsic first 
order reaction and the autocatalytic reaction. Then 
the apparent Arrhenius parameters derived from 
formal first order description as made by Brill and 
Sinditskii are composed of those of the individual 
reaction parameters. A sort of average is obtained. 
This can be seen if one uses the reciprocal times 
to reach a certain conversion (calculated with 
equation (12) for EL and equation (20) for ML) as 


Table 15. Times t,p, to reach preset energy losses EL for HNF lot 2 and 3. 


1% EL 3% EL 5% EL 
TC Lot 2 Lot 3 Lot 2 Lot 3 Lot 2 Lot 3 
[d] [a] [d] [a] [d] [a] [d] [a] [d] [a] [d] [a] 
20 44.4 21.9 62.8 29.4 71.7 32.9 
25 17.9 9.51 25.4 12.8 29.1 14.4 
30 7.40 4.24 10.6 5.76 12.2 6.50 
40 1.38 0.91 2.00 1.25 2.31 1.42 
50 0.28 0.21 0.42 0.30 0.48 0.34 
60 23.6 19.8 35.2 28.3 40.9 32.4 
70 5.82 5.45 8.8 7.91 10.29 9.12 
75 2.98 2.94 4.53 4,3 5.31 4.97 
80 1.55 1.61 2.38 2.37 2.79 275 
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ML [%] 


HNF lot 1 
80°C 
—* Mittel 


—#— 80/1 
—*— 80/2 


time [d] 


10 12 14 16 18 


mil/g gas generation 


of HNF at 
100°C 


time [min] 


600 800 


200 400 


Figure 26. Comparison of experimental data obtained for mass loss at 80 °C, this work, and on the right 


gas generation at 100 °C, figure taken from Koroban." 


rate constants and determines from them the formal 
first order Arrhenius parameters. The apparent EF, 
has then values between £,, and E,,, the actual 
value depending on conversion. This is shown with 
data in Table 17 and their graphical representation 
in Figure 30. If a broader temperature range is 
used for the calculation of ¢, then the process with 
greater activation energy is somewhat pronounced 
and the determined formal first order E,,, values 
are a bit greater, as can be seen in Table 17. 


Aninstructive comparison of Arrhenius parameters 
is achieved by plotting the corresponding 


GG [ml/g] 


HNF 
gas generation at 100°C 


apparatur loading density 0.002 g/ml 


GF =3 ml/g 

GGe = 680 ml/g 

k1 = 1.8804 E-5 + 2.0 E-6 1/min 
k2 = 1.1461 E-2 + 2.9 E-4 1/min 
S? = 7.12737 (ml/g) 

R? = 0.99905 


0 100 


reaction rate constants together as in Figures 31 
to 33 for the first order and the autocatalytic 
rate constants. From Figure 31 one recognizes 
that the Arrhenius parameters of Brill and de 
Klerk deviate significantly from the others. The 
Arrhenius parameters of Sinditskii correspond 
well with those obtained from heat generation 
data, but the correlation with mass loss and gas 
generation derived values is still acceptable. 
Figure 32 shows only the first order rate constants 
from autocatalytic descriptions. Gas generation 
and mass loss correspond very well and also the 


© measurement 
model 


time [min] 


300 400 500 


Figure 27. Modelling of the data" with the autocatalytic model. 
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HNF lot 3 


60°C, air 


time [d] 


12 15 18 21 24 27 30 


KLERK et al.: HYDRAZINIUM NITROFORMATE (ANF) 
JTAC V64, 973 (2001) 


dQ/dt [uW/g] 


time [h] 


180 


Figure 28. Comparison of heat generation rates measured with HNF at 60 °C, on the left this work, and 


on the right figure taken from W. de Klerk.* 


HG derived data are in good agreement. Only 
FOI-E8 deviates somewhat because of a somewhat 
small pre-exponential factor. In Figure 33 the rate 
constants k, from the autocatalytic part can be seen. 
Again gas generation and mass loss correspond 
well and the group of HG derived values also. For 
both cases, k, and ky, there is a certain difference 
between the group of HG data and the data group 
from mass loss and gas generation data. The reason 
may be the different loading densities applied. 
Table 18 lists the situation for the measurements 
of this work. With mass loss the loading density 


HNF lot 3 


12 15 18 21 


was relatively small with about 0.08 g cm? for 
the most measurements. Additionally one has 
some escaping of decomposition gases due to the 
nature of the method. In HGR measurements one 
has sealed ampoules and the average of loading 
density was about 0.43 g cm? and for the 60 °C 
measurements some 0.68 g cm. The loading 
density during ML measurements is nominally 
only 18% of that of HG measurements, and it may 
be accounted smaller because of the not really 
sealed ML vials. Koroban'’ reports a remarkable 
dependence of HNF decomposition behaviour 


2000 


HNF-E8 

1500 + 80 75 
= 1000 + 
~ 500 + 
i.) 
=] 
2 70 
= 65 

) 55 60 
0 2 4 6 8 
Time (days) 


Figure 29. Comparison of heat generation rates measured with HNF at several temperatures, on the left 
this work, and on the right measurements from FOI, Grindsjén Research Center, Tumba, Sweden,’” used 


with permission. 
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Eas, 
1, (Ts) = Z,, (7s) Xp ea = ,/F (Tg): Zo, exp(—Edg, / RTs) = J F (Ty) kg, (Ty) (22) 
Ss 


2 2 
Z,°€xp(—Edg, / RT.)-2-p-A-Q,- aa 2pd-O.: Ne sl Pipa 
r.)= Edy, _ Edy, i Edy, \(23) 
os pc, (L.-T +L, /c,) p-0,-(L,—Th+Ly/¢,) 2-RT. 
lp linear burning rate of strand 
kor reaction rate constant for heat release in solid phase at combustion 
Lib pre-exponential factor for burning rate (with weak dependence on T;) 
E,or activation energy for decomposition, heat release in solid phase at combustion 
Zor pre-exponential factor for decomposition, heat release in solid phase at comb. 
F proportionality factor, depended on material properties and weakly on Tg; 
Ts surface temperature of strand (can be measured by micro thermocouples) 
Ty strand temperature not heated 
QO.” heat of decomposition reaction 
p mass density 
A heat conductivity 
Cp specific heat 
Li heat of melting 
R general gas constant 


20 
Eaty [kJ/mol HNF lot 2 Ig(Zty [1/d]) 
F 19.6 
conversion followed by energy loss 
apparent Arrhenius parameters Eat and Zty —?— Eaty-80°C 
as function of conversion -- K - - Eaty-500°C 
evaluation with two temperature range —--—- Ea of first order reaction 19.2 
20°C to 80°C and 20°C to 500°C ——-— Ea of autocatalytic reaction 
P< —’— Ig(Zty)-80°C 
oe ee - +X - -Ig(Zty)-500°C 18.8 
18.4 
18 
conversion [%] 
17.6 
0 10 20 30 40 50 60 70 80 90 100 


Figure 30. Apparent Arrhenius parameters are obtained if one does not resolve the underlying reactions. 
Here with the real Arrhenius parameters of autocatalytic description the times t,g;, have been calculated 
as function of temperature up to given conversions. The temperature dependence of these times can 
formally be described with on set of Arrhenius parameters E,,, and Z,. These two parameters are a mix of 
the underlying parameters obtained by autocatalytic description. 
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Table 16. Compilation of Arrhenius data obtained for HNF from this work (APP type S) and from 
literature. 


Temp. Description with model ‘Q: first order + autocatalytic’ 
Lot Data type range 


[°C] Ey/kI mol 1g(Z, (vd), EMO! toc, [ad] Rs 

HNF 1 Mass loss 50-80 166.2+10 21.9241.5 0.993 158.645 23.20+0.8 0.998 

HNF2 Heat 60-75 139.3418 17.9442.8 0.983 1284410 19.0841. 0.994 
generation 

HNF3 Heat 60-75 132.447 16.82+1.1 0.994 116.645 17.38+0.8 0.997 
generation 


Description with kinetics of ‘first order’ or of single step 
Mass loss, 


TGA 100-150 226.5+17 26.37+2.3 (0.97)? de Klerk et al., 2003° 
Gas Koroban et al., 1979! 
aceaten. 10200 176.6 + 21 23.50 3.1 0.972 Ig(Z, [ml g | d')= 26.33 + 3.1 for 

a aes initial decomposition only 

-jJump, 
using 7 is i Williams, Brill, 1995" 
induction Pe we Pale pre-exp. factor estimated 
time 


Sinditskii et al., 2002!° 
1g(Z,, [mm s ']) =5.95+0.18 


Flame front 292-477 135.2+6.8 17.59+0.47 0.956 — Sinditskii, 2006,'° equation (23) 


Flame front 292-477 140.7+4.6 17.93+0.41 0.957 


Following results based on personal communication.'> HNF samples from APP, type C and E. Measurements 
from ref. 12. Data description with model ‘Q: first order + autocatalytic’ 


Etype Heat. 60-80 137.3 + 6.6 17.1141.0 0.993 117.7410.6 17.61+1.6 0.976 
generation 
Ctype Heat 60-80 134.0+ 9.8 17.0741.5 0.979 118.2410 17.5641.6 0.971 


generation ee 
Complete data from Koroban et al.,'' recalculated from data given in Figure 3 in ref. 11, k,/670 [ml g '] and C)*k, 


Gas 
eneration 


70-100 176.6 +21 23.50+43.1 0.972 135.0406 19.6440.1 0.9999 
Data from Koroban et al.'' and by personnel communication’® 


Gas = n 7 - 


Table 17. Apparent Arrhenius parameters obtained by calculation of times typ, as function of temperature 
to reach given conversions (energy losses) for HNF lot 2. First order E,, = 139.3 kJ mol ', autocatalytic 
Ey = 128.4 kJ mol’. 


Conversion by energy loss _29 to 500 °C 20 to 80 °C 
EL [%] Ey [kJ/mol] Ig(Zy [1/d]) Eqy [kJ mol] Ig(Zy [1/d]) 
134.4 19.709 132.77 19.446 
5 132.2 19.117 131.23 18.963 
10 131.5 18.924 130.79 18.815 
20 131.0 18.774 130.44 18.688 
40 130.5 18.626 130.14 18.564 
60 130.3 18.529 129.94 18.480 
80 123.0 18.432 129.75 18.393 
90 129.8 18.364 129.64 18.332 
95 129.7 18.309 129.54 18.282 
99 129.5 18.207 129.37 18.186 
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Ig(ky [1/d]) 


HNF 
first order reaction rate constants k, 


ICT-ML, Lot1, extrapolated 


Sinditskii 


ICT-HG, Lot2, extrapolated 


——ICT-ML, L1 

—— ICT-HG, L3 
Sinditskii 

—-—-Koroban Sindit Brill 


de Klerk 


—— FOI-HG, E8 ~——~FOI-HG, C9 
ICT-ML-100 


ICT-L2-100 


0.001 0.0015 0.002 


ICT-HG, L2 


Koroban Fig. 3 


Koroban, GG 


ICT-HG, L2 


ICT-ML, L1 


4/T [4/K] 


0.0025 0.003 0.0035 


Figure 31. Comparison of reaction rate constants k, of first order reaction in autocatalytic descriptions 
and apparent first order rate constants from other descriptions. 


on loading density. The differences found in the 
Arrhenius parameters of this work could originate 
in such a dependence on loading density. But also 
the probing itself of the two methods could result 
in differences. The two methods gas evolution and 
release of heat may see a different weighting of the 
ongoing complex decomposition processes and or 
may recognize different rate determining steps. 


Conclusion 


The thermal-chemical stability of the energetic 
substance HNF is such that it cannot be seen as 


non-handling substance. At room temperature it 
can be stored for several years as pure material. 
Substances like unstabilized NC and ADN have 
better stability behaviourthan HNF. The determined 
relatively small values of activation energies, 
from mass loss (lot 1) data 166 kJ mol! and with 
heat generation between 132 and 140 kJ mol! for 
lots 2 and 3, each for the intrinsic decomposition 
reaction, reflect the reduced thermal stability of 
HNF. HNF decomposes autocatalytically. The 
activation energies for this part are smaller than 
for the intrinsic decomposition. A comparative 


Table 18. Loading densities (LD) of the ML and HGR measurements. Mass loss was determined in vials 
with volume of 12.6 ml, closed by ground stoppers. The HGR were determined in glass ampoules closed 
with a rubber sealing, volume 3 ml. E means amount weighed in, free volume in vials and ampoules 
calculated for the situation at begin of the measurements with HNF density of 1.91 g cm®. 


Temp./ ML, Lot 1 HGR, Lot 2 HGR, Lot 3 

bi Elg : eos ae Els gem? yee Els LD/ g cm’ val fom? 
80 0.5 0.040 12.338 
75 1.0 0.079 12.076 1.29 0.430 2.325 1.30 0.433 2.319 
70 1.0 0.079 12.076 1.30 0.433 2.319 1.25 0.417 2.346 
65 1.0 0.079 12.076 1.28 0.427 2.330 1.28 0.427 2.330 
60 1.0 0.079 12.076 1.93 0.643 1.990 2.14 0.713 1.880 
50 1.0 0.079 12.076 
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Ig(k1 [1/d]) RE 


first order reaction rate constants k; 
Koroban, Sinditskii 


Koroban ~N ICT-HG, L3 


ICT-ML, L1 ICT-HG, L2 

ICT-HG, L2 

ICT-HG, L3 FOI-HG, E8 

Koroban Fig. 3 FOI-HG. C9 ICT-ML, L1 
— = Koroban Sindit : 
—— FOI-HG, E8 
== =FOI-HG, C9 


4/T [1/K] 


0.0026 0.0027 0.0028 0.0029 0.003 0.0031 0.0032 


Figure 32. First order reaction rate constants k, in autocatalytic descriptions only. 


discussion of the Arrhenius parameters found here evaluations based on heat generation. Nevertheless 

and those which can be found in the literature the stability of HNF is such that it may be stored 

shows a good agreement between the methods at 25 °C for several years until reaching a mass 

of the same probing type: this means agreement loss of 1%. The autocatalytic decomposition of 

between mass loss and gas generation and between HNF is not delayed as with nitrocellulose, it starts 
1.5 


Ig(ka [1/d]) HNFE 


autocatalytic reaction rate constants k2 


0.5 
Koroban, Fig. 3 ICT-HG, L3 
0 
-0.5 
FOI-HG, C9 
a ICT-ML,L1  Koroban, Sinditskii 
He ICT-HG, L2 a ICT-HG, L2 
——— ICT-HG, L3 
-2 Koroban Fig. 3 ICT-ML, L1 
== = Koroban Sindit 
“2.5 = —_FOLHG, E8 
3 ———~ FOI-HG, C9 
4/7 [1/K] 
-3.5 
0.0026 0.0027 0.0028 0.0029 0.003 0.0031 0.0032 


Figure 33. Comparison of reaction rate constants k, of autocatalytic reactions in autocatalytic 
descriptions. 
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from the beginning on. An additional source of 
risk is the high heat generation rate HNF is able to 
show. This limits non-stabilized HNF to handling 
only in containers of small geometric dimensions. 
The use in military ammunition with the demand 
to handle it up to 71 °C (highest environmental 
temperature according to STANAG 2895) does 
not seem advisable. The high heat generation rates 
during HNF decomposition would be a source 
of a non-manageable danger. Further to this the 
possibility of stabilizing HNF seems questionable. 
HNF decomposes via very reactive species, and 
probably they cannot be trapped fast enough by 
stabilizers added in the usual amounts. 
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HNF 


ADN 
NC 
RDX 


CL20 


HMxX 


AP 
AN 
ARC™ 


TAM™ 


List of abbreviations 


hydrazinium nitroformate (older name 
also hydrazine trinitromethane) 


ammonium dinitramide 

nitrocellulose 
1,3,5-trinitro-1,3,5-triaaza-cyclohexane 
(Hexogen, Research Development 
eXplosive or Royal Demolition 
eXplosive ) 


China Lake explosive no. 20, 
laboratory short name of HNIW; 
2,4,6,8,10,12-hexanitro-2,4,6,8,10,12- 
hexaaza-tetracyclo-[5.5.0.0°?.07!"]- 
dodecane 


octogen, cyclo-tetramethylene- 
tetranitramine, or 1,3,5,7-tetranitro- 
1,3,5,7-tetraaza-cyclooctane(High 
Melting Explosive or Her Majesty’s 
Explosive) 


ammonium perchlorate 
ammonium nitrate 


Accelerating Rate Calorimeter, 
determines adiabatic self heat rate 


Thermal Activity Monitor, 
microcalorimeter, determines heat 
generation rate 


mass loss 

energy loss 

heat generation rate (dQ/d?) 

heat generation (Q) 

gas generation 

vacuum stability test 

thermo gravimetric analysis 
activation energy from mass loss data 


pre-exponential factor from mass loss 
data; ky, = Zr *exp(—Egmr/RT) 
activation energy from heat generation 
or HGR data 

pre-exponential factor from HG or 
HGR data; kg = Zo *exp(—E,,0/RT) 
time to reach a certain value of mass 
loss 


time to reach a certain value of energy 
loss 


heat of explosion 
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Caution 


The experimentation with, and the use of, pyrotechnic materials can be dangerous and may require licences 
or permits in certain countries; it is felt to be important for the reader to be duly cautioned. Without the 
proper training and experience no one should ever experiment with or use pyrotechnic materials. Also, the 
amount of information presented in this Journal is not a substitute for necessary training and experience, 
nor does it remove the relevant application of national or local laws and regulations. 


A major effort has been undertaken to review all articles for correctness. However it is possible that 
errors remain. It is the responsibility of the reader to verify any information herein before applying that 
information in situations where death, injury or property damage could result. 
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Information for Readers 
Editorial Policy 


Articles accepted for publication in the Journal of Pyrotechnics can be on any technical subject in 
pyrotechnics. However, a strong preference will be given to articles reporting on research (conducted by 
professional or serious individual experimenters) and to review articles (either at an advanced or tutorial 
level). Both long and short articles will be gladly accepted. Also, responsible letters commenting on past 
Journal articles will be published along with responses by the authors. 


Back Issues 


Back issues of the Journal will be kept in print permanently as reference material. In addition all past 
articles and issues are available at the Journal of Pyrotechnics - http://archive.jpyro.com. 


New Publication Procedures 


From the next issue, Issue 27, the JPyro Board has decided that the Journal will be first published in 
electronic form and that the hard copy will be published annually. This has been driven by three main 
factors - the needs of authors to have their papers published in a reasonable timescale, the success of the 
JPyro archive, and the costs and time required to produce hard copies. In addition the modern publishing 
and research worlds increasing rely on early electronic publication to provide information to others in 
the most timely manner and, as appropriate, to be able to claim Intellectual Property rights on material 
published (which are taken from the date of publication). 


In future articles accepted for publication in the Journal will be published, in the first instance, electronically 
on the JPyro archive website (http://archive.jpyro.com ) within a target timescale of three weeks from 
acceptance (i.e. when all editing is complete) and eight weeks from original submission. All articles will be 
dated and published in acceptance order. The timescales rely, of course, on referees and authors responding 
in a reasonable timescale to any queries raised. 


In this way no article will be held up awaiting the collation of all the other material necessary for economic 
and efficient production of the hard copy. Instead, a hard copy of the Journal will be available to those who 
require it, in most cases annually - incorporating all the articles published within the previous year - for an 
additional charge. In addition, the JPyro archive will continue to offer online purchase of single articles 
from the Journal and other publications, or time-limited access to the entire archive. 


If you have any questions regarding the new approach please do contact any of the production team 


Publication Frequency 


Articles in the Journal of Pyrotechnics will appear first electronically at http://archive.jpyro.com and will 
subsequently be produced in hard copy annually in the early spring of the year following publication. 


Subscriptions 
The following subscription types are available 


¢ Electronic only - access to the electronic versions of articles published during the subscription period 
indefinitely 


¢ Electronic and Hard copy - as above + hard copy when published 
e Single article access - indefinite access to the specified article 


¢ Access to the entire archive - annual access to the full archive 
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